-
A Bayesian method for point source polarization estimation
Authors:
D. Herranz,
F. Argüeso,
L. Toffolatti,
A. Manjón-García,
M. López-Caniego
Abstract:
The estimation of the polarization $P$ of extragalactic compact sources in Cosmic Microwave Background images is a very important task in order to clean these images for cosmological purposes -- as, for example, to constrain the tensor-to-scalar ratio of primordial fluctuations during inflation -- and also to obtain relevant astrophysical information about the compact sources themselves in a frequ…
▽ More
The estimation of the polarization $P$ of extragalactic compact sources in Cosmic Microwave Background images is a very important task in order to clean these images for cosmological purposes -- as, for example, to constrain the tensor-to-scalar ratio of primordial fluctuations during inflation -- and also to obtain relevant astrophysical information about the compact sources themselves in a frequency range, $ν\sim 10$--$200$ GHz, where observations have only very recently started to be available. In this paper we propose a Bayesian maximum a posteriori (MAP) approach estimation scheme which incorporates prior information about the distribution of the polarization fraction of extragalactic compact sources between 1 and 100 GHz. We apply this Bayesian scheme to white noise simulations and to more realistic simulations that include CMB intensity, Galactic foregrounds and instrumental noise with the characteristics of the QUIJOTE experiment Wide Survey at 11 GHz. Using these simulations, we also compare our Bayesian method with the frequentist Filtered Fusion method that has been already used in WMAP data and in the \emph{Planck} mission. We find that the Bayesian method allows us to decrease the threshold for a feasible estimation of $P$ to levels below $\sim 100$ mJy (as compared to $\sim 500$ mJy that was the equivalent threshold for the frequentist Filtered Fusion). We compare the bias introduced by the Bayesian method and find it to be small in absolute terms. Finally, we test the robustness of the Bayesian estimator against uncertainties in the prior and in the flux density of the sources. We find that the Bayesian estimator is robust against moderate changes in the parameters of the prior and almost insensitive to realistic errors in the estimated photometry of the sources.
△ Less
Submitted 5 February, 2021;
originally announced February 2021.
-
Cosmological constraints with the sub-millimetre galaxies Magnification Bias after large scale bias corrections
Authors:
J. Gonzalez-Nuevo,
M. M. Cueli,
L. Bonavera,
A. Lapi,
M. Migliaccio,
F. Argüeso,
L. Toffolatti
Abstract:
The study of the magnification bias produced on high-redshift sub-millimetre galaxies by foreground galaxies through the analysis of the cross-correlation function was recently demonstrated as an interesting independent alternative to the weak-lensing shear as a cosmological probe. In the case of the proposed observable, most of the cosmological constraints mainly depend on the largest angular sep…
▽ More
The study of the magnification bias produced on high-redshift sub-millimetre galaxies by foreground galaxies through the analysis of the cross-correlation function was recently demonstrated as an interesting independent alternative to the weak-lensing shear as a cosmological probe. In the case of the proposed observable, most of the cosmological constraints mainly depend on the largest angular separation measurements. Therefore, we aim to study and correct the main large-scale biases that affect foreground and background galaxy samples to produce a robust estimation of the cross-correlation function. Then we analyse the corrected signal to derive updated cosmological constraintsWe measured the large-scale, bias-corrected cross-correlation functions using a background sample of H-ATLAS galaxies with photometric redshifts > 1.2 and two different foreground samples (GAMA galaxies with spectroscopic redshifts or SDSS galaxies with photometric ones, both in the range 0.2 < z < 0.8). These measurements are modelled using the traditional halo model description that depends on both halo occupation distribution and cosmological parameters. We then estimated these parameters by performing a Markov chain Monte Carlo under multiple scenarios to study the performance of this observable and how to improve its results. After the large-scale bias corrections, we obtain only minor improvements with respect to the previous magnification bias results, mainly confirming their conclusions: a lower bound on $Ω_m > 0.22$ at $95\%$ C.L. and an upper bound $σ_8 < 0.97$ at $95\%$ C.L. (results from the $z_{spec}$ sample). However, by combining both foreground samples into a simplified tomographic analysis, we were able to obtain interesting constraints on the $Ω_m$-$σ_8$ plane as follows: $Ω_m= 0.50_{- 0.20}^{+ 0.14}$ and $σ_8= 0.75_{- 0.10}^{+ 0.07}$ at 68\% CL.
△ Less
Submitted 23 December, 2020; v1 submitted 29 July, 2020;
originally announced July 2020.
-
Confusion noise due to clustered extragalactic point sources. Application of logarithmic cumulants for parameter estimation
Authors:
Francisco Argüeso,
Diego Herranz,
Luigi Toffolatti,
Joaquín González-Nuevo
Abstract:
The calculation of the characteristic function of the signal fluctuations due to clustered astrophysical sources is performed in this paper. For the typical case of power-law differential number counts and two-point angular correlation function, we present an extension of Zolotarev's theorem that allows us to compute the cumulants of the logarithm of the absolute value of the intensity. As a test,…
▽ More
The calculation of the characteristic function of the signal fluctuations due to clustered astrophysical sources is performed in this paper. For the typical case of power-law differential number counts and two-point angular correlation function, we present an extension of Zolotarev's theorem that allows us to compute the cumulants of the logarithm of the absolute value of the intensity. As a test, simulations based on recent observations of radio galaxies are then carried out, showing that these cumulants can be very useful for determining the fundamental parameters defining the number counts and the correlation. If the angular correlation scale of the observed source population is known, the method presented here is able to obtain estimators of the amplitude and slope of the power-law number counts with mean absolute errors that are one order of magnitude better than previous techniques, that did not take into account the correlation. Even if the scale of correlation is not well known, the method is able to estimate it and still performs much better than if the effect of correlations is not considered.
△ Less
Submitted 4 June, 2019;
originally announced June 2019.
-
SHALOS: Statistical Herschel-ATLAS Lensed Objects Selection
Authors:
J. González-Nuevo,
S. L. Suárez Gómez,
L. Bonavera,
F. Sánchez-Lasheras,
F. Argüeso,
L. Toffolatti,
D. Herranz,
C. González-Gutiérrez,
F. García Riesgo,
F. J. de Cos Juez
Abstract:
The statistical analysis of large sample of strong lensing events can be a powerful tool to extract astrophysical and/or cosmological valuable information. However, the number of such events is still relatively low, mostly because of the lengthily observational validation process on individual events. In this work we propose a new methodology with a statistical selection approach in order to incre…
▽ More
The statistical analysis of large sample of strong lensing events can be a powerful tool to extract astrophysical and/or cosmological valuable information. However, the number of such events is still relatively low, mostly because of the lengthily observational validation process on individual events. In this work we propose a new methodology with a statistical selection approach in order to increase by a factor of $\sim 5$ the number of such events. Although the methodology can be applied to address several selection problems, it has particular benefits in the case of the identification of strongly lensed galaxies: objectivity, minimal initial constrains in the main parameter space, preservation of the statistical properties. The proposed methodology is based on the Bhattacharyya distance as a measure of the similarity between probability distributions of properties of two different cross-matched galaxies. The particular implementation for the aim of this work is called SHALOS and it combines the information of four different properties of the pair of galaxies: angular separation, luminosity percentile, redshift and optical/sub-mm flux density ratio. The SHALOS method provided a ranked list of strongly lensed galaxies. The number of candidates for the final associated probability, $P_{tot}>0.7$, is 447 with an estimated mean amplification factor of 3.12 for an halo with a typical cluster mass. Additional statistical properties of the SHALOS candidates, as the correlation function or the source number counts, are in agreement with previous results indicating the statistical lensing nature of the selected sample.
△ Less
Submitted 15 March, 2019;
originally announced March 2019.
-
Planck 2018 results. II. Low Frequency Instrument data processing
Authors:
Planck Collaboration,
Y. Akrami,
F. Argüeso,
M. Ashdown,
J. Aumont,
C. Baccigalupi,
M. Ballardini,
A. J. Banday,
R. B. Barreiro,
N. Bartolo,
S. Basak,
K. Benabed,
J. -P. Bernard,
M. Bersanelli,
P. Bielewicz,
L. Bonavera,
J. R. Bond,
J. Borrill,
F. R. Bouchet,
F. Boulanger,
M. Bucher,
C. Burigana,
R. C. Butler,
E. Calabrese,
J. -F. Cardoso
, et al. (126 additional authors not shown)
Abstract:
We present a final description of the data-processing pipeline for the Planck, Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters.…
▽ More
We present a final description of the data-processing pipeline for the Planck, Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters. We provide a brief pedagogical introduction to the complete pipeline, as well as a detailed description of the important changes implemented. Self-consistency of the pipeline is demonstrated using dedicated simulations and null tests. We present the final version of the LFI full sky maps at 30, 44, and 70 GHz, both in temperature and polarization, together with a refined estimate of the Solar dipole and a final assessment of the main LFI instrumental parameters.
△ Less
Submitted 11 September, 2018; v1 submitted 17 July, 2018;
originally announced July 2018.
-
Planck intermediate results. LIV. The Planck Multi-frequency Catalogue of Non-thermal Sources
Authors:
Planck Collaboration,
Y. Akrami,
F. Argüeso,
M. Ashdown,
J. Aumont,
C. Baccigalupi,
M. Ballardini,
A. J. Banday,
R. B. Barreiro,
N. Bartolo,
S. Basak,
K. Benabed,
J. -P. Bernard,
M. Bersanelli,
P. Bielewicz,
L. Bonavera,
J. R. Bond,
J. Borrill,
F. R. Bouchet,
C. Burigana,
R. C. Butler,
E. Calabrese,
J. Carron,
H. C. Chiang,
C. Combet
, et al. (116 additional authors not shown)
Abstract:
This paper presents the Planck Multi-frequency Catalogue of Non-thermal (i.e. synchrotron-dominated) Sources (PCNT) observed between 30 and 857 GHz by the ESA Planck mission. This catalogue was constructed by selecting objects detected in the full mission all-sky temperature maps at 30 and 143 GHz, with a signal-to-noise ratio (S/N)>3 in at least one of the two channels after filtering with a part…
▽ More
This paper presents the Planck Multi-frequency Catalogue of Non-thermal (i.e. synchrotron-dominated) Sources (PCNT) observed between 30 and 857 GHz by the ESA Planck mission. This catalogue was constructed by selecting objects detected in the full mission all-sky temperature maps at 30 and 143 GHz, with a signal-to-noise ratio (S/N)>3 in at least one of the two channels after filtering with a particular Mexican hat wavelet. As a result, 29400 source candidates were selected. Then, a multi-frequency analysis was performed using the Matrix Filters methodology at the position of these objects, and flux densities and errors were calculated for all of them in the nine Planck channels. The present catalogue is the first unbiased, full-sky catalogue of synchrotron-dominated sources published at millimetre and submillimetre wavelengths and constitutes a powerful database for statistical studies of non-thermal extragalactic sources, whose emission is dominated by the central active galactic nucleus. Together with the full multi-frequency catalogue, we also define the Bright Planck Multi-frequency Catalogue of Non-thermal Sources PCNTb, where only those objects with a S/N>4 at both 30 and 143 GHz were selected. In this catalogue 1146 compact sources are detected outside the adopted Planck GAL070 mask; thus, these sources constitute a highly reliable sample of extragalactic radio sources. We also flag the high-significance subsample PCNThs, a subset of 151 sources that are detected with S/N>4 in all nine Planck channels, 75 of which are found outside the Planck mask adopted here. The remaining 76 sources inside the Galactic mask are very likely Galactic objects.
△ Less
Submitted 11 September, 2018; v1 submitted 23 February, 2018;
originally announced February 2018.
-
Statistics of the fractional polarisation of extragalactic dusty sources in Planck HFI maps
Authors:
L. Bonavera,
J. Gonzalez-Nuevo,
B. De Marco,
F. Argüeso,
L. Toffolatti
Abstract:
We estimate the average fractional polarisation at 143, 217 and 353 GHz of a sample of 4697 extragalactic dusty sources by applying stacking technique. The sample is selected from the second version of the Planck Catalogue of Compact Sources at 857 GHz, avoiding the region inside the Planck Galactic mask (fsky ~ 60 per cent). We recover values for the mean fractional polarisation at 217 and 353 GH…
▽ More
We estimate the average fractional polarisation at 143, 217 and 353 GHz of a sample of 4697 extragalactic dusty sources by applying stacking technique. The sample is selected from the second version of the Planck Catalogue of Compact Sources at 857 GHz, avoiding the region inside the Planck Galactic mask (fsky ~ 60 per cent). We recover values for the mean fractional polarisation at 217 and 353 GHz of (3.10 \pm 0.75) per cent and (3.65 \pm 0.66) per cent, respectively, whereas at 143 GHz we give a tentative value of (3.52 \pm 2.48) per cent. We discuss the possible origin of the measured polarisation, comparing our new estimates with those previously obtained from a sample of radio sources. We test different distribution functions and we conclude that the fractional polarisation of dusty sources is well described by a log-normal distribution, as determined in the radio band studies. For this distribution we estimate μ_{217GHz} = 0.3 \pm 0.5 (that would correspond to a median fractional polarisation of Π_{med} = (1.3 \pm 0.7) per cent) and μ_{353GHz} = 0.7 \pm 0.4 (Π_{med} = (2.0 \pm 0.8) per cent), σ_{217GHz} = 1.3 \pm 0.2 and σ_{353GHz} = 1.1 \pm 0.2. With these values we estimate the source number counts in polarisation and the contribution given by these sources to the CMB B-mode angular power spectrum at 217, 353, 600 and 800 GHz. We conclude that extragalactic dusty sources might be an important contaminant for the primordial B-mode at frequencies > 217 GHz.
△ Less
Submitted 29 May, 2017;
originally announced May 2017.
-
Statistics of the fractional polarisation of compact radio sources in Planck maps
Authors:
Laura Bonavera,
Joaquin González-Nuevo,
Francisco Argüeso,
Luigi Toffolatti
Abstract:
In this work we apply the stacking technique to estimate the average fractional polarisation from 30 to 353 GHz of a primary sample of 1560 compact sources - essentially all radio sources - detected in the 30 GHz Planck all-sky map and listed in the second version of the Planck Catalogue of Compact Sources (PCCS2). We divide our primary sample in two subsamples according to whether the sources lay…
▽ More
In this work we apply the stacking technique to estimate the average fractional polarisation from 30 to 353 GHz of a primary sample of 1560 compact sources - essentially all radio sources - detected in the 30 GHz Planck all-sky map and listed in the second version of the Planck Catalogue of Compact Sources (PCCS2). We divide our primary sample in two subsamples according to whether the sources lay (679 sources) or not (881 sources) inside the sky region defined by the Planck Galactic mask (fsky ~ 60 per cent) and the area around the Magellanic Clouds. We find that the average fractional polarisation of compact sources is approximately constant (with frequency) in both samples (with a weighted mean over all the channels of 3.08 per cent outside and 3.54 per cent inside the Planck mask). In the sky region outside the adopted mask, we also estimate the μ and σ parameters for the log-normal distribution of the fractional polarisation, finding a weighted mean value over all the Planck frequency range of 1.0 for σ and 0.7 for μ (that would imply a weighted mean value for the median fractional polarisation of 1.9 per cent).
△ Less
Submitted 26 May, 2017; v1 submitted 29 March, 2017;
originally announced March 2017.
-
Planck 2015 results. XXVI. The Second Planck Catalogue of Compact Sources
Authors:
Planck Collaboration,
P. A. R. Ade,
N. Aghanim,
F. Argüeso,
M. Arnaud,
M. Ashdown,
J. Aumont,
C. Baccigalupi,
A. J. Banday,
R. B. Barreiro,
N. Bartolo,
E. Battaner,
C. Beichman,
K. Benabed,
A. Benoît,
A. Benoit-Lévy,
J. -P. Bernard,
M. Bersanelli,
P. Bielewicz,
J. J. Bock,
H. Böhringer,
A. Bonaldi,
L. Bonavera,
J. R. Bond,
J. Borrill
, et al. (218 additional authors not shown)
Abstract:
The Second Planck Catalogue of Compact Sources is a catalogue of sources detected in single-frequency maps from the full duration of the Planck mission and supersedes previous versions of the Planck compact source catalogues. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. Compact sources detected in the lower frequency channels are assigned to the PC…
▽ More
The Second Planck Catalogue of Compact Sources is a catalogue of sources detected in single-frequency maps from the full duration of the Planck mission and supersedes previous versions of the Planck compact source catalogues. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. Compact sources detected in the lower frequency channels are assigned to the PCCS2, while at higher frequencies they are assigned to one of two sub-catalogues, the PCCS2 or PCCS2E, depending on their location on the sky. The first of these catalogues covers most of the sky and allows the user to produce subsamples at higher reliabilities than the target 80% integral reliability of the catalogue. The PCCS2E contains sources detected in sky regions where the diffuse emission makes it difficult to quantify the reliability of the detections. Both the PCCS2 and PCCS2E include polarization measurements, in the form of polarized flux densities, or upper limits, and orientation angles for all seven polarization-sensitive Planck channels. The improved data-processing of the full-mission maps and their reduced noise levels allow us to increase the number of objects in the catalogue, improving its completeness for the target 80 % reliability as compared with the previous versions, the PCCS and ERCSC catalogues.
△ Less
Submitted 15 February, 2016; v1 submitted 8 July, 2015;
originally announced July 2015.
-
Planck 2013 results. XXVIII. The Planck Catalogue of Compact Sources
Authors:
Planck Collaboration,
P. A. R. Ade,
N. Aghanim,
F. Argüeso,
C. Armitage-Caplan,
M. Arnaud,
M. Ashdown,
F. Atrio-Barandela,
J. Aumont,
C. Baccigalupi,
A. J. Banday,
R. B. Barreiro,
J. G. Bartlett,
E. Battaner,
K. Benabed,
A. Benoît,
A. Benoit-Lévy,
J. -P. Bernard,
M. Bersanelli,
P. Bielewicz,
J. Bobin,
J. J. Bock,
A. Bonaldi,
L. Bonavera,
J. R. Bond
, et al. (218 additional authors not shown)
Abstract:
The Planck Catalogue of Compact Sources (PCCS) is the catalogue of sources detected in the first 15 months of Planck operations, the "nominal" mission. It consists of nine single-frequency catalogues of compact sources, both Galactic and extragalactic, detected over the entire sky. The PCCS covers the frequency range 30--857\,GHz with higher sensitivity (it is 90% complete at 180 mJy in the best c…
▽ More
The Planck Catalogue of Compact Sources (PCCS) is the catalogue of sources detected in the first 15 months of Planck operations, the "nominal" mission. It consists of nine single-frequency catalogues of compact sources, both Galactic and extragalactic, detected over the entire sky. The PCCS covers the frequency range 30--857\,GHz with higher sensitivity (it is 90% complete at 180 mJy in the best channel) and better angular resolution (from ~33' to ~5') than previous all-sky surveys in this frequency band. By construction its reliability is >80% and more than 65% of the sources have been detected at least in two contiguous Planck channels. In this paper we present the construction and validation of the PCCS, its contents and its statistical characterization.
△ Less
Submitted 25 November, 2013; v1 submitted 20 March, 2013;
originally announced March 2013.
-
Extragalactic Compact Sources in the Planck Sky and their Cosmological Implications
Authors:
Luigi Toffolatti,
Carlo Burigana,
Francisco Argueso,
Jose M. Diego
Abstract:
The Planck satellite has proved to be a very successful mission, which has been operating flawlessly for more than 36 months. Its main purpose is to map the anisotropies of the CMB radiation at 9 frequencies, between 30 and 857 GHz, with unprecedented resolution and sensitivity. After a description of the properties and cosmological aims of Planck, we review the results obtained during the first 1…
▽ More
The Planck satellite has proved to be a very successful mission, which has been operating flawlessly for more than 36 months. Its main purpose is to map the anisotropies of the CMB radiation at 9 frequencies, between 30 and 857 GHz, with unprecedented resolution and sensitivity. After a description of the properties and cosmological aims of Planck, we review the results obtained during the first 1.6 full sky surveys relative to compact extragalactic sources and their implications. The most recent and efficient methods and filters for compact source detection in the presence of Gaussian noise and CMB anisotropies are discussed. The first surveys of Planck have enabled the detection of about 200 galaxy clusters, by the SZ effect, and from several hundreds to many thousands of extragalactic point sources (EPS), presented in the Planck ERCSC. We review the results on the SZ effect in galaxy clusters. We discuss the possibilities of future Planck data to study clusters through the SZ effect and the interesting possibilities of combining these data with X ray observations. These EPS can be classified into two main source classes: radio sources, at intermediate to high z; far IR sources, typically at very z. Planck number counts of EPS at LFI frequencies are in agreement with WMAP counts. A clear steepening of the mean spectral index of bright radio sources is observed above 70 GHz, that can be interpreted in terms of a break frequency in their spectra. For nearby dusty galaxies, Planck observations find evidence of colder dust, with T < 20 K, than has previously been found. More recently, the number counts of bright local dusty galaxies and of synchrotron sources have been measured by Planck at submm wavelengths. Planck has provided interesting results on the angular distribution of CIB anisotropies, allowing us to put new constraints on clustering properties of dusty galaxies at high z.
△ Less
Submitted 14 February, 2013;
originally announced February 2013.
-
Forecasts on the contamination induced by unresolved point sources in primordial non-Gaussianity beyond Planck
Authors:
A. Curto,
M. Tucci,
J. Gonzalez-Nuevo,
L. Toffolatti,
E. Martinez-Gonzalez,
F. Argueso,
A. Lapi,
M. Lopez-Caniego
Abstract:
In this paper we present forecasts of the contamination on different shapes of the primordial non-Gaussianity fnl parameter -- detectable on future Cosmic Microwave Background (CMB) high--resolution anisotropy maps -- produced by unresolved extragalactic point sources at frequencies of cosmological interest (45--375 GHz). We consider two scenarios: an ideal (noiseless) mission and a possible futur…
▽ More
In this paper we present forecasts of the contamination on different shapes of the primordial non-Gaussianity fnl parameter -- detectable on future Cosmic Microwave Background (CMB) high--resolution anisotropy maps -- produced by unresolved extragalactic point sources at frequencies of cosmological interest (45--375 GHz). We consider two scenarios: an ideal (noiseless) mission and a possible future space-borne satellite, with instrumental characteristics similar to the ones proposed for the Cosmic Origins Explorer (COrE). The local, equilateral, orthogonal and flat shapes are considered in both temperature (intensity) and polarized emission data. The angular power spectrum and bispectrum of extragalactic point sources are estimated by state-of-the-art models of source number counts. The impact of all the most relevant (far--IR and radio selected) source populations on these shapes at COrE frequencies is studied. The results of this analysis show that unresolved extragalactic point sources should not induce a very relevant non-Gaussian signal in the frequency range 100--200 GHz, thus not preventing a correct estimate of the CMB primordial fnl parameter. Polarization information allows one to significantly reduce the error--bars in the fnl parameter and the bias induced by unresolved sources and, hence, to widen the range of frequencies for fnl studies. On the contrary, at nu < 100 GHz or nu > 225 GHz, important non-Gaussian deviations in CMB anisotropy maps are expected due to unresolved extragalactic sources.
△ Less
Submitted 20 March, 2013; v1 submitted 8 January, 2013;
originally announced January 2013.
-
Planck intermediate results. VII. Statistical properties of infrared and radio extragalactic sources from the Planck Early Release Compact Source Catalogue at frequencies between 100 and 857 GHz
Authors:
Planck Collaboration,
P. A. R. Ade,
N. Aghanim,
F. Argüeso,
M. Arnaud,
M. Ashdown,
F. Atrio-Barandela,
J. Aumont,
C. Baccigalupi,
A. Balbi,
A. J. Banday,
R. B. Barreiro,
E. Battaner,
K. Benabed,
A. Benoît,
J. -P. Bernard,
M. Bersanelli,
M. Bethermin,
R. Bhatia,
A. Bonaldi,
J. R. Bond,
J. Borrill,
F. R. Bouchet,
C. Burigana,
P. Cabella
, et al. (166 additional authors not shown)
Abstract:
(abridged for arXiv) We make use of the Planck all-sky survey to derive number counts and spectral indices of extragalactic sources -- infrared and radio sources -- from the Planck Early Catalogue (ERCSC) at 100 to 857GHz. Our sample contains, after the 80% completeness cut, between 122 and 452 and sources, with flux densities above 0.3 and 1.9Jy at 100 and 857GHz, over about 31 to 40% of the sky.…
▽ More
(abridged for arXiv) We make use of the Planck all-sky survey to derive number counts and spectral indices of extragalactic sources -- infrared and radio sources -- from the Planck Early Catalogue (ERCSC) at 100 to 857GHz. Our sample contains, after the 80% completeness cut, between 122 and 452 and sources, with flux densities above 0.3 and 1.9Jy at 100 and 857GHz, over about 31 to 40% of the sky. Using Planck HFI, all the sources have been classified as either dust-dominated or synchrotron-dominated on the basis of their spectral energy distributions (SED). Our sample is thus complete, flux-limited and color-selected to differentiate between the two populations. We find an approximately equal number of synchrotron and dusty sources between 217 and 353GHz; at 353GHz or higher (or 217GHz and lower) frequencies, the number is dominated by dusty (synchrotron) sources, as expected. For most of the sources, the spectral indices are also derived. We provide for the first time counts of bright sources from 353 to 857GHz and the contributions from dusty and synchrotron sources at all HFI frequencies in the key spectral range where these spectra are crossing. The observed counts are in the Euclidean regime. The number counts are compared to previously published data (earlier Planck, Herschel, BLAST, SCUBA, LABOCA, SPT, and ACT) and models taking into account both radio or infrared galaxies. We derive the multi-frequency Euclidean level and compare it to WMAP, Spitzer and IRAS results. The submillimetre number counts are not well reproduced by current evolution models of dusty galaxies, whereas the millimetre part appears reasonably well fitted by the most recent model for synchrotron-dominated sources. Finally we provide estimates of the local luminosity density of dusty galaxies, providing the first such measurements at 545 and 857GHz.
△ Less
Submitted 22 January, 2013; v1 submitted 19 July, 2012;
originally announced July 2012.
-
Compact source detection in multi-channel microwave surveys: from SZ clusters to polarized sources
Authors:
Diego Herranz,
Francisco Argüeso,
Pedro Carvalho
Abstract:
In this paper we describe the state-of-the art status of multi-frequency detection techniques for compact sources in microwave astronomy. From the simplest cases where the spectral behaviour is well-known (i.e. thermal SZ clusters) to the more complex cases where there is little a priori information (i.e. polarized radio sources) we will review the main advances and the most recent results in the…
▽ More
In this paper we describe the state-of-the art status of multi-frequency detection techniques for compact sources in microwave astronomy. From the simplest cases where the spectral behaviour is well-known (i.e. thermal SZ clusters) to the more complex cases where there is little a priori information (i.e. polarized radio sources) we will review the main advances and the most recent results in the detection problem.
△ Less
Submitted 17 April, 2012;
originally announced April 2012.
-
Planck Early Results: Statistical properties of extragalactic radio sources in the Planck Early Release Compact Source Catalogue
Authors:
Planck Collaboration,
P. A. R. Ade,
N. Aghanim,
F. Argüeso,
M. Arnaud,
M. Ashdown,
J. Aumont,
C. Baccigalupi,
A. Balbi,
A. J. Banday,
R. B. Barreiro,
J. G. Bartlett,
E. Battaner,
K. Benabed,
J. -P. Bernard,
M. Bersanelli,
R. Bhatia,
A. Bonaldi,
L. Bonavera,
J. R. Bond,
J. Borrill,
F. R. Bouchet,
M. Bucher,
C. Burigana,
P. Cabella
, et al. (179 additional authors not shown)
Abstract:
The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts (dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in very good agreement with e…
▽ More
The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts (dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in very good agreement with estimates based on WMAP data, being somewhat deeper at 30 and 70 GHz, and somewhat shallower at 44 GHz. Planck's source counts at 143 and 217 GHz join smoothly with the fainter ones provided by the SPT and ACT surveys over small fractions of the sky. An analysis of source spectra, exploiting Planck's uniquely broad spectral coverage, finds clear evidence of a steepening of the mean spectral index above about 70 GHz. This implies that, at these frequencies, the contamination of the CMB power spectrum by radio sources below the detection limit is significantly lower than previously estimated.
△ Less
Submitted 28 February, 2011; v1 submitted 11 January, 2011;
originally announced January 2011.
-
A Bayesian technique for the detection of point sources in CMB maps
Authors:
F. Argueso,
E. Salerno,
D. Herranz,
J. L. Sanz,
E. E. Kuruoglu,
K. Kayabol
Abstract:
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 numb…
▽ More
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.
△ Less
Submitted 7 January, 2011;
originally announced January 2011.
-
Filter Design for the Detection/Estimation of the Modulus of a Vector. Application to Polarization Data
Authors:
F. Argueso,
J. L. Sanz,
D. Herranz
Abstract:
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 devel…
▽ More
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.
△ Less
Submitted 4 January, 2011;
originally announced January 2011.
-
Polarization of the WMAP Point Sources
Authors:
M. Lopez-Caniego,
M. Massardi,
J. Gonzalez-Nuevo,
L. Lanz,
D. Herranz,
G. De Zotti,
J. L. Sanz,
F. Argueso
Abstract:
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…
▽ More
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.
△ Less
Submitted 24 September, 2009;
originally announced September 2009.
-
Detection/estimation of the modulus of a vector. Application to point source detection in polarization data
Authors:
F. Argueso,
J. L. Sanz,
D. Herranz,
M. Lopez-Caniego,
J. Gonzalez-Nuevo
Abstract:
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, o…
▽ More
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.
△ Less
Submitted 4 June, 2009;
originally announced June 2009.
-
Statistical properties of extragalactic sources in the New Extragalactic WMAP Point Source (NEWPS) catalogue
Authors:
J. González-Nuevo,
M. Massardi,
F. Argüeso,
D. Herranz,
L. Toffolatti,
J. L. Sanz,
M. López-Caniego,
G. De Zotti
Abstract:
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ópez-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 yi…
▽ More
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ópez-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.
△ Less
Submitted 16 November, 2007;
originally announced November 2007.
-
Non-blind catalogue of extragalactic point sources from the Wilkinson Microwave Anisotropy Probe (WMAP) first 3--year survey data
Authors:
M. Lopez-Caniego,
J. Gonzalez-Nuevo,
D. Herranz,
M. Massardi,
J. L. Sanz,
G. De Zotti,
L. Toffolatti,
F. Argueso
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 January, 2007;
originally announced January 2007.
-
Wavelets on the sphere. Application to the detection problem
Authors:
J. L. Sanz,
D. Herranz,
M. Lopez-Caniego,
F. Argueso
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 September, 2006;
originally announced September 2006.
-
Statistical analysis of undetected point sources in cosmic microwave background maps
Authors:
F. Argueso,
J. L. Sanz,
R. B. Barreiro,
D. Herranz,
J. Gonzalez-Nuevo
Abstract:
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 analys…
▽ More
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.
△ Less
Submitted 13 September, 2006;
originally announced September 2006.
-
Comparison of filters for the detection of point sources in Planck simulations
Authors:
M. Lopez-Caniego,
D. Herranz,
J. Gonzalez-Nuevo,
J. L. Sanz,
R. B. Barreiro,
P. Vielva,
F. Argueso,
L. Toffolatti
Abstract:
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 sourc…
▽ More
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.
△ Less
Submitted 8 June, 2006;
originally announced June 2006.
-
The Mexican Hat Wavelet Family. Application to point source detection in CMB maps
Authors:
J. Gonzalez-Nuevo,
F. Argueso,
M. Lopez-Caniego,
L. Toffolatti,
J. L. Sanz,
P. Vielva,
D. Herranz
Abstract:
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…
▽ More
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.
△ Less
Submitted 21 April, 2006; v1 submitted 18 April, 2006;
originally announced April 2006.
-
Extragalactic source contributions to arcminute-scale Cosmic Microwave Background anisotropies
Authors:
L. Toffolatti,
M. Negrello,
J. González-Nuevo,
G. De Zotti,
L. Silva,
G. L. Granato,
F. Argüeso
Abstract:
The possible contributions of the various classes of extragalactic sources (including, in addition to the canonical radio sources, GHz Peaked Spectrum sources, advection-dominated sources, starburst galaxies, high-redshift proto-spheroidal galaxies) to the arcminute scale fluctuations measured by the CBI, BIMA, and ACBAR experiments are discussed. At 30 GHz, fluctuations due to radio sources und…
▽ More
The possible contributions of the various classes of extragalactic sources (including, in addition to the canonical radio sources, GHz Peaked Spectrum sources, advection-dominated sources, starburst galaxies, high-redshift proto-spheroidal galaxies) to the arcminute scale fluctuations measured by the CBI, BIMA, and ACBAR experiments are discussed. At 30 GHz, fluctuations due to radio sources undetected by ancillary low-frequency surveys may be higher than estimated by the CBI and BIMA groups. High-redshift dusty galaxies, whose fluctuations may be strongly enhanced by the effect of clustering, could contribute to the BIMA excess signal, and dominate at 150 GHz (the ACBAR frequency). Moreover, in the present data situation, the dust emission of these high-redshift sources set an unavoidable limit to the detection of primordial CMB anisotropies at high multipoles, even at frequencies as low as $\simeq 30$ GHz. It is concluded that the possibility that the excess power at high multipoles is dominated by unsubtracted extragalactic sources cannot be ruled out. On the other hand, there is room for a contribution from the Sunyaev-Zeldovich effect within clusters of galaxies, with a density fluctuation amplitude parameter $σ_8$ consistent with the values preferred by current data.
△ Less
Submitted 21 April, 2005; v1 submitted 25 October, 2004;
originally announced October 2004.
-
Predictions on the angular power spectrum of clustered extragalactic point sources at CMB frequencies from flat and all--sky 2D-simulations
Authors:
J. González-Nuevo,
L. Toffolatti,
F. Argüeso
Abstract:
We present predictions on the angular power spectrum of CMB fluctuations due to extragalactic point sources (EPS) by using a method for simulating realistic 2D distributions of clustered EPS. Both radio and far--IR selected source populations are taken into account. To analyze different clustering scenarios, we exploit angular power spectra of EPS, $P(k)$, estimated either by data coming from cu…
▽ More
We present predictions on the angular power spectrum of CMB fluctuations due to extragalactic point sources (EPS) by using a method for simulating realistic 2D distributions of clustered EPS. Both radio and far--IR selected source populations are taken into account. To analyze different clustering scenarios, we exploit angular power spectra of EPS, $P(k)$, estimated either by data coming from currently available surveys or by means of theoretical predictions. By adopting the source number counts predicted by the Toffolatti et al. (1998) evolution model -- capable of accounting well for the available data at radio cm wavelengths -- we are able to reproduce current data on the two--point angular correlation functions, $w(θ)$, of radio sources. We can confirm that the detection of primordial CMB anisotropies is not hampered by undetected clustered sources at frequencies $\leq 150-200$ GHz. On the other hand, our current findings show that at higher frequencies the clustering signal could severely reduce the detectability of intrinsic CMB anisotropies, thus confirming previous theoretical predictions. We also show that unsubtracted EPS can account for the excess signal at high multipoles detected by recent CMB anisotropy experiments. Moreover, the additional power due to the clustering of sources gives rise to a small but not negligible contribution to the same excess signal. As a final result, we also present an example of a currently feasible {\it realistic map} of EPS at 70 GHz, by taking into account data on bright detected sources as well as the previously quoted model for number counts.
△ Less
Submitted 12 November, 2004; v1 submitted 27 May, 2004;
originally announced May 2004.
-
Tests of Gaussianity
Authors:
A. M. Aliaga,
E. Martinez-Gonzalez,
L. Cayon,
F. Argueso,
J. L. Sanz,
R. B. Barreiro,
J. E. Gallegos
Abstract:
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 m…
▽ More
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.
△ Less
Submitted 24 October, 2003;
originally announced October 2003.
-
Goodness-of-fit tests of Gaussianity: constraints on the cumulants of the MAXIMA data
Authors:
A. M. Aliaga,
E. Martinez-Gonzalez,
L. Cayon,
F. Argueso,
J. L. Sanz,
R. B. Barreiro
Abstract:
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 R…
▽ More
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.
△ Less
Submitted 22 September, 2003;
originally announced September 2003.
-
Contributions of point extragalactic sources to the Cosmic Microwave Background bispectrum
Authors:
F. Argüeso,
J. González-Nuevo,
L. Toffolatti
Abstract:
All the analyses of Cosmic Microwave Background (CMB) temperature maps up--to--date show that CMB anisotropies follow a Gaussian distribution. On the other hand, astrophysical foregrounds which hamper the detection of the CMB angular power spectrum, are not Gaussian distributed on the sky. Therefore, they should give a sizeable contribution to the CMB bispectrum. In fact, the first year data of…
▽ More
All the analyses of Cosmic Microwave Background (CMB) temperature maps up--to--date show that CMB anisotropies follow a Gaussian distribution. On the other hand, astrophysical foregrounds which hamper the detection of the CMB angular power spectrum, are not Gaussian distributed on the sky. Therefore, they should give a sizeable contribution to the CMB bispectrum. In fact, the first year data of the Wilkinson Microwave Anisotropy Probe (WMAP) mission have allowed the {\it first} detection of the extragalactic source contribution to the CMB bispectrum at 41 GHz and, at the same time, much tighter limits than before to non--Gaussian primordial fluctuations. In view of the above and for achieving higher precision in current and future CMB measurements of non--Gaussianity, in this paper we discuss a comprehensive assessment of the bispectrum due to either uncorrelated and clustered extragalactic point sources in the whole frequency interval around the CMB intensity peak. Our calculations, based on current cosmological evolution models for sources, show that the reduced angular bispectrum due to point sources, $b_{ps}$, should be detectable in all WMAP and Planck frequency channels. We also find agreement with the results on $b_{ps}$ at 41 GHz coming from the analysis of the first year WMAP data. Moreover, by comparing $b_{ps}$ with the primordial reduced CMB bispectrum, we find that only the peak value of the primordial bispectrum (which appears at $l\simeq 200$) results greater than $b_{ps}$ in a frequency window around the intensity peak of the CMB. The amplitude of this window basically depends on the capability of the source detection algorithms (i.e., on the achievable flux detection limit, $S_{lim}$, for sources).
△ Less
Submitted 8 July, 2003;
originally announced July 2003.
-
Goodness-of-Fit Tests to study the Gaussianity of the MAXIMA data
Authors:
L. Cayon,
F. Argueso,
E. Martinez-Gonzalez,
J. L. Sanz
Abstract:
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 Goodn…
▽ More
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.
△ Less
Submitted 9 June, 2003;
originally announced June 2003.
-
COBE-DMR constraints on the nonlinear coupling parameter: a wavelet based method
Authors:
L. Cayon,
E. Martinez-Gonzalez,
F. Argueso,
A. J. Banday,
K. M. Gorski
Abstract:
Nonlinearity introduced in slow-roll inflation will produce weakly non-Gaussian CMB temperature fluctuations. We have simulated non-Gaussian large scale CMB maps (including COBE-DMR constraints) introducing an additional quadratic term in the gravitational potential.The amount of nonlinearity being controlled by the so called nonlinear coupling parameter f_{nl}. An analysis based on the Spherica…
▽ More
Nonlinearity introduced in slow-roll inflation will produce weakly non-Gaussian CMB temperature fluctuations. We have simulated non-Gaussian large scale CMB maps (including COBE-DMR constraints) introducing an additional quadratic term in the gravitational potential.The amount of nonlinearity being controlled by the so called nonlinear coupling parameter f_{nl}. An analysis based on the Spherical Mexican Hat wavelet was applied to these and to the COBE-DMR maps. Skewness values obtained at several scales were combined into a Fisher discriminant. Comparison of the Fisher discriminant distributions obtained for different nonlinear coupling parameters with the COBE-DMR values, sets a constraint of abs(f_{nl})<1100 at the 68% confidence level. This new constraint being tighter than the one previously obtained by using the bispectrum by Komatsu et al. (2002).
△ Less
Submitted 18 November, 2002;
originally announced November 2002.
-
The performance of spherical wavelets to detect non-Gaussianity in the CMB sky
Authors:
E. Martinez-Gonzalez,
J. E. Gallegos,
F. Argueso,
L. Cayon,
J. L. Sanz
Abstract:
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…
▽ More
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.
△ Less
Submitted 17 July, 2002; v1 submitted 14 November, 2001;
originally announced November 2001.
-
Spherical Mexican Hat wavelet: an application to detect non-Gaussianity in the COBE-DMR maps
Authors:
L. Cayon,
J. L. Sanz,
E. Martinez-Gonzalez,
A. J. Banday,
F. Argueso,
J. E. Gallegos,
K. M. Gorski,
G. Hinshaw
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 7 May, 2001;
originally announced May 2001.
-
Isotropic Wavelets: a Powerful Tool to Extract Point Sources from CMB Maps
Authors:
L. Cayon,
J. L. Sanz,
R. B. Barreiro,
E. Matinez-Gonzalez,
P. Vielva,
L. Toffolatti,
J. Silk,
J. M. Diego,
F. Argueso
Abstract:
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 met…
▽ More
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.
△ Less
Submitted 22 December, 1999;
originally announced December 1999.
-
Analysis of CMB maps with 2D wavelets
Authors:
J. L. Sanz,
R. B. Barreiro,
L. Cayon,
E. Martinez-Gonzalez,
G. A. Ruiz,
F. J. Diaz,
F. Argueso,
J. Silk,
L. Toffolatti
Abstract:
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.…
▽ More
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.
△ Less
Submitted 29 September, 1999;
originally announced September 1999.
-
Wavelets Applied to CMB Maps: a Multiresolution Analysis for Denoising
Authors:
J. L. Sanz,
F. Argueso,
L. Cayon,
E. Martinez-Gonzalez,
R. B. Barreiro,
L. Toffolatti
Abstract:
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 b…
▽ More
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$.
△ Less
Submitted 23 June, 1999;
originally announced June 1999.
-
On Planck simulations: towards `2-nd order' analyses
Authors:
C. Burigana,
D. Maino,
N. Mandolesi,
F. Villa,
L. Valenziano,
M. Bersanelli,
L. Danese,
L. Toffolatti,
F. Argüeso
Abstract:
We simulate Planck observations by adopting a detailed model of the microwave sky including monopole, dipole, anisotropies of the cosmic microwave background (CMB) and galactic and extragalactic foregrounds. We estimate the impact of main beam optical aberrations on CMB anisotropy measurements in presence of extragalactic source fluctuations and we discuss the main implications for the Planck te…
▽ More
We simulate Planck observations by adopting a detailed model of the microwave sky including monopole, dipole, anisotropies of the cosmic microwave background (CMB) and galactic and extragalactic foregrounds. We estimate the impact of main beam optical aberrations on CMB anisotropy measurements in presence of extragalactic source fluctuations and we discuss the main implications for the Planck telescope design. By analysing the dipole pattern, we quantify the Planck performance in the determination of CMB spectral distortion parameters in presence of foreground contaminations.
△ Less
Submitted 9 March, 1999;
originally announced March 1999.
-
Extragalactic Radio Sources and CMB anisotropies
Authors:
L. Toffolatti,
G. De Zotti,
F. Argüeso,
C. Burigana
Abstract:
Confusion noise due to extragalactic sources is a fundamental astrophysical limitation for experiments aimed at accurately determining the power spectrum of the Cosmic Microwave Background (CMB) down to arcmin angular scales and with a sensitivity $ΔT/T \simeq 10^{-6}$. At frequencies $\lsim 200-300$ GHz, the most relevant extragalactic foreground hampering the detection of intrinsic CMB anisotr…
▽ More
Confusion noise due to extragalactic sources is a fundamental astrophysical limitation for experiments aimed at accurately determining the power spectrum of the Cosmic Microwave Background (CMB) down to arcmin angular scales and with a sensitivity $ΔT/T \simeq 10^{-6}$. At frequencies $\lsim 200-300$ GHz, the most relevant extragalactic foreground hampering the detection of intrinsic CMB anisotropies is constituted by radio loud Active Galactic Nuclei (AGN), including ``flat--spectrum'' radiogalaxies, quasars, BL-LACs and blazars. We review our present understanding of astrophysical properties, spectra, and number counts of the above classes of sources. We also study the angular power spectrum of fluctuations due both to Poisson distributed and clustered radio sources and give preliminary predictions on the power spectrum of their polarized components. Furthermore, we discuss the capabilities of future space missions (NASA's MAP, Bennett et al. 1995; ESA's Planck Surveyor, Bersanelli et al. 1996) in studying bright radio sources over an almost unexplored frequency interval where spectral signatures, essential for the understanding of the physical processes, show up.
△ Less
Submitted 3 March, 1999; v1 submitted 24 February, 1999;
originally announced February 1999.
-
The Planck Surveyor mission: astrophysical prospects
Authors:
G. De Zotti,
L. Toffolatti,
F. Argüeso,
R. D. Davies,
P. Mazzotta,
R. B. Partridge,
G. F. Smoot,
N. Vittorio
Abstract:
Although the Planck Surveyor mission is optimized to map the cosmic microwave background anisotropies, it will also provide extremely valuable information on astrophysical phenomena. We review our present understanding of Galactic and extragalactic foregrounds relevant to the mission and discuss on one side, Planck's impact on the study of their properties and, on the other side, to what extent…
▽ More
Although the Planck Surveyor mission is optimized to map the cosmic microwave background anisotropies, it will also provide extremely valuable information on astrophysical phenomena. We review our present understanding of Galactic and extragalactic foregrounds relevant to the mission and discuss on one side, Planck's impact on the study of their properties and, on the other side, to what extent foreground contamination may affect Planck's ability to accurately determine cosmological parameters. Planck's multifrequency surveys will be unique in their coverage of large areas of the sky (actually, of the full sky); this will extend by two or more orders of magnitude the flux density interval over which mm/sub-mm counts of extragalactic sources can be determined by instruments already available (like SCUBA) or planned for the next decade (like the LSA-MMA or the space mission FIRST), which go much deeper but over very limited areas. Planck will thus provide essential complementary information on the epoch-dependent luminosity functions. Bright radio sources will be studied over a poorly explored frequency range where spectral signatures, essential to understand the physical processes that are going on, show up. The Sunyaev-Zeldovich effect, with its extremely rich information content, will be observed in the direction of a large number of rich clusters of Galaxies. Thanks again to its all sky coverage, Planck will provide unique information on the structure and on the emission properties of the interstellar medium in the Galaxy. At the same time, the foregrounds are unlikely to substantially limit Planck's ability to measure the cosmological signals. Even measurements of polarization of the primordial Cosmic Microwave background fluctuations appear to be feasible.
△ Less
Submitted 27 March, 1999; v1 submitted 6 February, 1999;
originally announced February 1999.