We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

In this paper we derive a full expression for the propagation of multiplicative and additive shape measurement biases into the cosmic shear power spectrum. In doing so we identify several new terms that are associated with selection effects, as well as cross-correlation terms between the multiplicative and additive biases and the shear field. The computation of the resulting bias in the shear power spectrum scales as the fifth power of the maximum multipole considered. Consequently the calculation is unfeasible for large l-modes, and the only tractable way to assess the full impact of shape measurement biases on cosmic shear power spectrum is through forward modelling of the effects. To linear order in bias parameters the shear power spectrum is only affected by the mean of the multiplicative bias field over a survey and the cross correlation between the additive bias field and the shear field. If the mean multiplicative bias is zero then second order convolutive terms are expected to be orders of magnitude smaller.

The primordial power spectrum describes the initial perturbations that seeded the large-scale structure we observe today. It provides an indirect probe of inflation or other structure-formation mechanisms. In this letter, we recover the primordial power spectrum from the Planck PR1 dataset, using our recently published algorithm PRISM. PRISM is a sparsity-based inversion method, that aims at recovering features in the primordial power spectrum from the empirical power spectrum of the cosmic microwave background (CMB). This ill-posed inverse problem is regularised using a sparsity prior on features in the primordial power spectrum in a wavelet dictionary. Although this non-parametric method does not assume a strong prior on the shape of the primordial power spectrum, it is able to recover both its general shape and localised features. As a results, this approach presents a reliable way of detecting deviations from the currently favoured scale-invariant spectrum. We applied PRISM to 100 simulated Planck data to investigate its performance on Planck-like data. We also tested the algorithm's ability to recover a small localised feature at $k \sim 0.125$ Mpc$^{-1}$, which caused a large dip at $\ell \sim 1800$ in the angular power spectrum. We then applied PRISM to the Planck PR1 power spectrum to recover the primordial power spectrum. We find no significant departures from the fiducial Planck PR1 near scale-invariant primordial power spectrum with $A_s=2.215\times10^{-9}$ and $n_s = 0.9624$.

Large-scale anomalies have been reported in CMB data with both WMAP and Planck data. These could be due to foreground residuals and or systematic effects, though their confirmation with Planck data suggests they are not due to a problem in the WMAP or Planck pipelines. If these anomalies are in fact primordial, then understanding their origin is fundamental to either validate the standard model of cosmology or to explore new physics. We investigate three other possible issues: 1) the trade-off between minimising systematics due to foreground contamination (with a conservative mask) and minimising systematics due to masking, 2) astrophysical secondary effects (the kinetic Doppler quadrupole and kinetic Sunyaev-Zel'dovich effect), and 3) secondary cosmological signals (the integrated Sachs-Wolfe effect). We address the masking issue by considering new procedures that use both WMAP and Planck to produce higher quality full-sky maps using the sparsity methodology (LGMCA maps). We show the impact of masking is dominant over that of residual foregrounds, and the LGMCA full-sky maps can be used without further processing to study anomalies. We consider four official Planck PR1 and two LGMCA CMB maps. Analysis of the observed CMB maps shows that only the low quadrupole and quadrupole-octopole alignment seem significant, but that the planar octopole, Axis of Evil, mirror parity and cold spot are not significant in nearly all maps considered. After subtraction of astrophysical and cosmological secondary effects, only the low quadrupole may still be considered anomalous, meaning the significance of only one anomaly is affected by secondary effect subtraction out of six anomalies considered. In the spirit of reproducible research all reconstructed maps and codes will be made available for download here http://www.cosmostat.org/anomaliesCMB.html.

The primordial power spectrum is an indirect probe of inflation or other structure-formation mechanisms. We introduce a new method, named \textbf{PRISM}, to estimate this spectrum from the empirical cosmic microwave background (CMB) power spectrum. This is a sparsity-based inversion method, which leverages a sparsity prior on features in the primordial spectrum in a wavelet dictionary to regularise the inverse problem. This non-parametric approach is able to reconstruct the global shape as well as localised features of spectrum accurately and proves to be robust for detecting deviations from the currently favoured scale-invariant spectrum. We investigate the strength of this method on a set of WMAP nine-year simulated data for three types of primordial spectra and then process the WMAP nine-year data as well as the Planck PR1 data. We find no significant departures from a near scale-invariant spectrum.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

Missions such as WMAP or Planck measure full-sky fluctuations of the cosmic microwave background and foregrounds, among which bright compact source emissions cover a significant fraction of the sky. To accurately estimate the diffuse components, the point-source emissions need to be separated from the data, which requires a dedicated processing. We propose a new technique to estimate the flux of the brightest point sources using a morphological separation approach: point sources with known support and shape are separated from diffuse emissions that are assumed to be sparse in the spherical harmonic domain. This approach is compared on both WMAP simulations and data with the standard local chi2 minimization, modelling the background as a low-order polynomial. The proposed approach generally leads to 1) lower biases in flux recovery, 2) an improved root mean-square error of up to 35% and 3) more robustness to background fluctuations at the scale of the source. The WMAP 9-year point-source-subtracted maps are available online.

The primordial power spectrum describes the initial perturbations in the Universe which eventually grew into the large-scale structure we observe today, and thereby provides an indirect probe of inflation or other structure-formation mechanisms. Here, we introduce a new method to estimate this spectrum from the empirical power spectrum of cosmic microwave background (CMB) maps. A sparsity-based linear inversion method, coined \textbf{PRISM}, is presented. This technique leverages a sparsity prior on features in the primordial power spectrum in a wavelet basis to regularise the inverse problem. This non-parametric approach does not assume a strong prior on the shape of the primordial power spectrum, yet is able to correctly reconstruct its global shape as well as localised features. These advantages make this method robust for detecting deviations from the currently favoured scale-invariant spectrum. We investigate the strength of this method on a set of WMAP 9-year simulated data for three types of primordial power spectra: a nearly scale-invariant spectrum, a spectrum with a small running of the spectral index, and a spectrum with a localised feature. This technique proves to easily detect deviations from a pure scale-invariant power spectrum and is suitable for distinguishing between simple models of the inflation. We process the WMAP 9-year data and find no significant departure from a nearly scale-invariant power spectrum with the spectral index $n_s = 0.972$. A high resolution primordial power spectrum can be reconstructed with this technique, where any strong local deviations or small global deviations from a pure scale-invariant spectrum can easily be detected.

We present a novel estimate of the cosmological microwave background (CMB) map by combining the two latest full-sky microwave surveys: WMAP nine-year and Planck PR1. The joint processing benefits from a recently introduced component separation method coined "local-generalized morphological component analysis'' (LGMCA) based on the sparse distribution of the foregrounds in the wavelet domain. The proposed estimation procedure takes advantage of the IRIS 100 micron as an extra observation on the galactic center for enhanced dust removal. We show that this new CMB map presents several interesting aspects: i) it is a full sky map without using any inpainting or interpolating method, ii) foreground contamination is very low, iii) the Galactic center is very clean, with especially low dust contamination as measured by the cross-correlation between the estimated CMB map and the IRIS 100 micron map, and iv) it is free of thermal SZ contamination.

The next generation of galaxy surveys, aiming to observe millions of galaxies, are expensive both in time and cost. This raises questions regarding the optimal investment of this time and money for future surveys. In a previous work, it was shown that a sparse sampling strategy could be a powerful substitute for the contiguous observations. However, in this previous paper a regular sparse sampling was investigated, where the sparse observed patches were regularly distributed on the sky. The regularity of the mask introduces a periodic pattern in the window function, which induces periodic correlations at specific scales. In this paper, we use the Bayesian experimental design to investigate a random sparse sampling, where the observed patches are randomly distributed over the total sparsely sampled area. We find that, as there is no preferred scale in the window function, the induced correlation is evenly distributed amongst all scales. This could be desirable if we are interested in specific scales in the galaxy power spectrum, such as the Baryonic Acoustic Oscillation (BAO) scales. However, for constraining the overall galaxy power spectrum and the cosmological parameters, there is no preference over regular or random sampling. Hence any approach that is practically more suitable can be chosen and we can relax the regular-grid condition for the distribution of the observed patches.

We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher l's. For l < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz through a physically motivated Bayesian component separation technique. At l >= 50, we employ a correlated Gaussian likelihood approximation based on angular cross-spectra derived from the 100, 143 and 217 GHz channels. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on cosmological parameters. We find good internal agreement among the high-l cross-spectra with residuals of a few uK^2 at l <= 1000. We compare our results with foreground-cleaned CMB maps, and with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. The best-fit LCDM cosmology is in excellent agreement with preliminary Planck polarisation spectra. The standard LCDM cosmology is well constrained by Planck by l <= 1500. For example, we report a 5.4 sigma deviation from n_s /= 1. Considering various extensions beyond the standard model, we find no indication of significant departures from the LCDM framework. Finally, we report a tension between the best-fit LCDM model and the low-l spectrum in the form of a power deficit of 5-10% at l <~ 40, significant at 2.5-3 sigma. We do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent dataset. (Abridged)

Recovering the Cosmic Microwave Background (CMB) from WMAP data requires galactic foreground emissions to be accurately separated out. Most component separation techniques rely on second order statistics such as Internal Linear Combination (ILC) techniques. In this paper, we present a new WMAP 9-year CMB map, with 15 arcmin resolution, which is reconstructed using a recently introduced sparse component separation technique, coined Local Generalized Morphological Component Analysis (LGMCA). LGMCA emphasizes on the sparsity of the components to be retrieved in the wavelet domain. We show that although derived from a radically different separation criterion ({i.e. sparsity), the LGMCA-WMAP 9 map and its power spectrum are fully consistent with their more recent estimates from WMAP 9.

The next generation of galaxy surveys will observe millions of galaxies over large volumes of the universe. These surveys are expensive both in time and cost, raising questions regarding the optimal investment of this time and money. In this work we investigate criteria for selecting amongst observing strategies for constraining the galaxy power spectrum and a set of cosmological parameters. Depending on the parameters of interest, it may be more efficient to observe a larger, but sparsely sampled, area of sky instead of a smaller contiguous area. In this work, by making use of the principles of Bayesian Experimental Design, we will investigate the advantages and disadvantages of the sparse sampling of the sky and discuss the circumstances in which a sparse survey is indeed the most efficient strategy. For the Dark Energy Survey (DES), we find that by sparsely observing the same area in a smaller amount of time, we only increase the errors on the parameters by a maximum of 0.45%. Conversely, investing the same amount of time as the original DES to observe a sparser but larger area of sky we can in fact constrain the parameters with errors reduced by 28%.

The cosmic microwave background (CMB) power spectrum is a powerful cosmological probe as it entails almost all the statistical information of the CMB perturbations. Having access to only one sky, the CMB power spectrum measured by our experiments is only a realization of the true underlying angular power spectrum. In this paper we aim to recover the true underlying CMB power spectrum from the one realization that we have without a need to know the cosmological parameters. The sparsity of the CMB power spectrum is first investigated in two dictionaries; Discrete Cosine Transform (DCT) and Wavelet Transform (WT). The CMB power spectrum can be recovered with only a few percentage of the coefficients in both of these dictionaries and hence is very compressible in these dictionaries. We study the performance of these dictionaries in smoothing a set of simulated power spectra. Based on this, we develop a technique that estimates the true underlying CMB power spectrum from data, i.e. without a need to know the cosmological parameters. This smooth estimated spectrum can be used to simulate CMB maps with similar properties to the true CMB simulations with the correct cosmological parameters. This allows us to make Monte Carlo simulations in a given project, without having to know the cosmological parameters. The developed IDL code, TOUSI, for Theoretical pOwer spectrUm using Sparse estImation, will be released with the next version of ISAP.

We introduce a new method for reconstructing the primordial power spectrum, $P(k)$, directly from observations of the Cosmic Microwave Background (CMB). We employ Singular Value Decomposition (SVD) to invert the radiation perturbation transfer function. The degeneracy of the multipole $\ell$ to wavenumber $k$ linear mapping is thus reduced. This enables the inversion to be carried out at each point along a Monte Carlo Markov Chain (MCMC) exploration of the combined $P(k)$ and cosmological parameter space. We present best--fit $P(k)$ obtained with this method along with other cosmological parameters.

The primordial power spectrum describes the initial perturbations in the Universe which eventually grew into the large-scale structure we observe today, and thereby provides an indirect probe of inflation or other structure-formation mechanisms. In this paper we will investigate the best scales the primordial power spectrum can be probed, in accordance with the knowledge about other cosmological parameters such as $\Omega_{b}$, $\Omega_{c}$, $\Omega_{\Lambda}$, $h$ and $\tau$. The aim is to find the most informative way of measuring the primordial power spectrum at different length scales, using different types of surveys and the information they provide for the desired cosmological parameters. We will find the optimal binning of the primordial power spectrum for this purpose, by making use of the Fisher matrix formalism. We will then find a statistically orthogonal basis for a set of cosmological parameters, mentioned above, and a set of bins of the primordial power spectrum to investigate the correlation between the two sets. For this purpose we make use of principal component analysis and Hermitian square root of the Fisher matrix. The surveys used in this project are Planck and SDSS(BRG), but the formalism can easily be extended to any windowed measurements of the perturbation spectrum.