The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

The efficient classification of different types of supernova is one of the most important problems for observational cosmology. However, spectroscopic confirmation of most objects in upcoming photometric surveys, such as the The Rubin Observatory Legacy Survey of Space and Time (LSST), will be unfeasible. The development of automated classification processes based on photometry has thus become crucial. In this paper we investigate the performance of machine learning (ML) classification on the final cosmological constraints using simulated lightcurves from The Supernova Photometric Classification Challenge, released in 2010. We study the use of different feature sets for the lightcurves and many different ML pipelines based on either decision tree ensembles or automated search processes. To construct the final catalogs we propose a threshold selection method, by employing a Bias-Variance tradeoff. This is a very robust and efficient way to minimize the Mean Squared Error. With this method we were able to get very strong cosmological constraints, which allowed us to keep ~75% of the total information in the type Ia SNe when using the SALT2 feature set and ~33% for the other cases (based on either the Newling model or on standard wavelet decomposition).

We present an efficient numerical code and conduct, for the first time, a null and model-independent CMB test of statistical isotropy using Multipole Vectors (MVs) at all scales. Because MVs are insensitive to the angular power spectrum $C_\ell$, our results are independent from the assumed cosmological model. We avoid a posteriori choices and use pre-defined ranges of scales $\ell\in[2,30]$, $\ell\in[2,600]$ and $\ell\in[2,1500]$ in our analyses. We find that all four masked Planck maps, from both 2015 and 2018 releases, are in agreement with statistical isotropy for $\ell\in[2,30]$, $\ell\in[2,600]$. For $\ell\in[2,1500]$ we detect anisotropies but this is indicative of simply the anisotropy in the noise: there is no anisotropy for $\ell < 1300$ and an increasing level of anisotropy at higher multipoles. Our findings of no large-scale anisotropies seem to be a consequence of avoiding \emph{a posteriori} statistics. We also find that the degree of anisotropy in the full sky (i.e. unmasked) maps vary enormously (between less than 5 and over 1000 standard deviations) among the different mapmaking procedures and data releases.

Peculiar velocities introduce correlations between supernova magnitudes, which implies that the supernova Hubble diagram residual carries information on both the matter power spectrum at the present time and its growth rate. By a combination of brute-force exact computations of likelihoods and Fisher matrix analysis, we investigate how this information, which comes from supernova data only, depends on different survey parameters such as covered area, depth, and duration. We show that, for a survey like The Rubin Observatory Legacy Survey of Space and Time (LSST) and a fixed redshift depth, the same observing time provides the same cosmological information whether one observes a larger area, or a smaller area during more years. We also show that although the peculiar velocity information is peaked in the range $z \in [0, 0.2]$, there is yet plenty of information in $z \in [0.2, 0.5]$, and for very high supernova number densities there is even more information in the latter range. We conclude that, after 5 years, LSST could measure $\sigma_8$ with an uncertainty of 0.17 with the current strategy, and that this could be improved to 0.09 if the supernova completeness is improved to 20%. Moreover, we forecast results considering the extra parameter $\gamma$, and show that this creates a non-linear degeneracy with $\sigma_8$ that makes the Fisher matrix analysis inadequate. Finally, we discuss the possibility of achieving competitive results with the current Zwicky Transient Facility.

We present the new method DALI (Derivative Approximation for LIkelihoods) for reconstructing and forecasting posteriors. DALI extends the Fisher Matrix formalism but allows for a much wider range of posterior shapes. While the Fisher Matrix formalism is limited to yield ellipsoidal confidence contours, our method can reproduce the often observed flexed, deformed or curved shapes of known posteriors. This gain in shape fidelity is obtained by expanding the posterior to higher order in derivatives with respect to parameters, such that non-Gaussianity in the parameter space is taken into account. The resulting expansion is positive definite and normalizable at every order. Here, we present the new technique, highlight its advantages and limitations, and show a representative application to a posterior of dark energy parameters from supernovae measurements.

Supernova Ia magnitude surveys measure the dimensionless luminosity distance $H_{0}D_{L}$. However, from the distances alone one cannot obtain quantities like $H(z)$ or the dark energy equation of state, unless further cosmological assumptions are imposed. Here we show that by measuring the power spectrum of density contrast and of peculiar velocities of supernovae one can estimate also $H(z)/H_{0}$ regardless of background or linearly perturbed cosmology and of galaxy-matter bias. This method, dubbed Clustering of Standard Candles (CSC) also yields the redshift distortion parameter $\beta(k,z)$ and the galaxy power spectrum. We forecast that LSST may be able to constrain $H(z)/H_{0}$ to 3--12\% in redshift bins of $\Delta z=0.1$ up $z=0.8$.

Observational cosmology is passing through a unique moment of grandeur with the amount of quality data growing fast. However, in order to better take advantage of this moment, data analysis tools have to keep up the pace. Understanding the effect of baryonic matter on the large-scale structure is one of the challenges to be faced in cosmology. In this work, we have thoroughly studied the effect of baryonic physics on different lensing statistics. Making use of the Magneticum Pathfinder suite of simulations we show that the influence of luminous matter on the 1-point lensing statistics of point sources is significant, enhancing the probability of magnified objects with $\mu>3$ by a factor of $2$ and the occurrence of multiple-images by a factor $5-500$ depending on the source redshift and size. We also discuss the dependence of the lensing statistics on the angular resolution of sources. Our results and methodology were carefully tested in order to guarantee that our uncertainties are much smaller than the effects here presented.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]

We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $\tau=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4\sigma$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged]

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 {\mu}K$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 {\mu}K$\cdot$arcmin.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

The abundances of dark matter halos in the universe are described by the halo mass function (HMF). It enters most cosmological analyses and parametrizes how the linear growth of primordial perturbations is connected to these abundances. Interestingly, this connection can be made approximately cosmology independent. This made it possible to map in detail its near-universal behavior through large-scale simulations. However, such simulations may suffer from systematic effects, especially if baryonic physics is included. In this paper we ask how well observations can constrain directly the HMF. The observables we consider are galaxy cluster number counts, galaxy cluster power spectrum and lensing of type Ia supernovae. Our results show that DES is capable of putting the first meaningful constraints on the HMF, while both Euclid and J-PAS can give stronger constraints, comparable to the ones from state-of-the-art simulations. We also find that an independent measurement of cluster masses is even more important for measuring the HMF than for constraining the cosmological parameters, and can vastly improve the determination of the halo mass function. Measuring the HMF could thus be used to cross-check simulations and their implementation of baryon physics. It could even, if deviations cannot be accounted for, hint at new physics.

Cosmic Microwave Background (CMB) experiments, such as WMAP and Planck, measure intensity anisotropies and build maps using a linearized formula for relating them to the temperature blackbody fluctuations. However, this procedure also generates a signal in the maps in the form of y-type distortions which is degenerate with the thermal Sunyaev Zel'dovich (tSZ) effect. These are small effects that arise at second-order in the temperature fluctuations not from primordial physics but from such a limitation of the map-making procedure. They constitute a contaminant for measurements of: our peculiar velocity, the tSZ and primordial y-distortions. They can nevertheless be well-modeled and accounted for. We show that the distortions arise from a leakage of the CMB dipole into the y-channel which couples to all multipoles, mostly affecting the range $\ell$ < ~400. This should be visible in Planck's y-maps with an estimated signal-to-noise ratio of about 12. We note however that such frequency-dependent terms carry no new information on the nature of the CMB dipole. This implies that the real significance of Planck's Doppler coupling measurements is actually lower than reported by the collaboration. Finally, we quantify the level of contamination in tSZ and primordial y-type distortions and show that it is above the sensitivity of proposed next generation CMB experiments.

The aberration and Doppler coupling effects of the Cosmic Microwave Background (CMB) were recently measured by the Planck satellite. The most straightforward interpretation leads to a direct detection of our peculiar velocity $\beta$, consistent with the measurement of the well-known dipole. In this paper we discuss the assumptions behind such interpretation. We show that Doppler-like couplings appear from two effects: our peculiar velocity and a second order large-scale effect due to the dipolar part of the gravitational potential. We find that the two effects are exactly degenerate but only if we assume second-order initial conditions from single-field Inflation. Thus, detecting a discrepancy in the value of $\beta$ from the dipole and the Doppler couplings implies the presence of a primordial non-Gaussianity. We also show that aberration-like signals likewise arise from two independent effects: our peculiar velocity and lensing due to a first order large-scale dipolar gravitational potential, independently on Gaussianity of the initial conditions. In general such effects are not degenerate and so a discrepancy between the measured $\beta$ from the dipole and aberration could be accounted for by a dipolar gravitational potential. Only through a fine-tuning of the radial profile of the potential it is possible to have a complete degeneracy with a boost effect. Finally we discuss that we also expect other signatures due to integrated second order terms, which may be further used to disentangle this scenario from a simple boost.

The supernova (SN) Hubble diagram residual contains valuable information on both the present matter power spectrum and its growth history. In this paper we show that this information can be retrieved with precision by combining both peculiar velocity and weak-lensing analysis on the data. To wit, peculiar velocity induces correlations on the nearby SN while lensing induces a non-Gaussian dispersion in faraway objects. We show that both effects have almost orthogonal degeneracies and discuss how they can be extracted simultaneously from the data. We analyze the JLA supernova catalog in a 14-dimensional parameter space, assuming a flexible growth-rate index $\gamma$. We arrive at the following marginalized constraints: $\sigma_8 = 0.65^{+0.23}_{-0.37}$ and $\gamma = 1.38^{+1.7}_{-0.65}$. Assuming instead GR as the correct gravitation theory (and thus $\gamma \equiv 0.55$), the constraints in $\sigma_8$ tighten further: $\sigma_8 = 0.40^{+0.21}_{-0.23}$. We show that these constraints complement well the ones obtained from other datasets and that they could improve substantially with more SNe.

Cosmic Microwave Background (CMB) experiments, such as WMAP and Planck, measure intensity anisotropies and build maps using a \emph{linearized} formula for relating them to the temperature blackbody fluctuations. However such a procedure also generates a signal in the maps in the form of y-type distortions, and thus degenerate with the thermal SZ (tSZ) effect. These are small effects that arise at second-order in the temperature fluctuations not from primordial physics but from such a limitation of the map-making procedure. They constitute a contaminant for measurements of: our peculiar velocity, the tSZ and of primordial y-distortions, but they can nevertheless be well-modelled and accounted for. We show that the largest distortions arises at high ell from a leakage of the CMB dipole into the y-channel which couples to all multipoles, but mostly affects the range ell <~ 400. This should be visible in Planck's y-maps with an estimated signal-to-noise ratio of about 9. We note however that such frequency-dependent terms carry no new information on the nature of the CMB dipole. This implies that the real significance of Planck's Doppler coupling measurements is actually lower than quoted. Finally, we quantify the relevance of the removal of such effects in order improve future measurements of tSZ and of primordial y-type distortions.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

The Planck satellite detectors are calibrated in the 2015 release using the "orbital dipole", which is the time-dependent dipole generated by the Doppler effect due to the motion of the satellite around the Sun. Such an effect has also relativistic time-dependent corrections of relative magnitude 10^(-3), due to coupling with the "solar dipole" (the motion of the Sun compared to the CMB rest frame), which are included in the data calibration by the Planck collaboration. We point out that such corrections are subject to a frequency-dependent multiplicative factor. This factor differs from unity especially at the highest frequencies, relevant for the HFI instrument. Since currently Planck calibration errors are dominated by systematics, to the point that polarization data is currently unreliable at large scales, such a correction can in principle be highly relevant for future data releases.

It has been pointed out recently that the quadrupole-octopole alignment in the CMB data is significantly affected by the so-called kinetic Doppler quadrupole (DQ), which is the temperature quadrupole induced by our proper motion. Assuming our velocity is the dominant contribution to the CMB dipole we have v/c = beta = (1.231 +/- 0.003) * 10^{-3}, which leads to a non-negligible DQ of order beta^2. Here we stress that one should properly take into account that CMB data are usually not presented in true thermodynamic temperature, which induces a frequency dependent boost correction. The DQ must therefore be multiplied by a frequency-averaged factor, which we explicitly compute for several CMB maps finding that it varies between 1.67 and 2.47. This is often neglected in the literature and turns out to cause a small but non-negligible difference in the significance levels of some quadrupole-related statistics. For instance the alignment significance in the SMICA 2013 map goes from 2.3sigma to 3.3sigma, with the frequency dependent DQ, instead of 2.9sigma ignoring the frequency dependence in the DQ. Moreover as a result of a proper DQ removal, the agreement across different map-making techniques is improved.

A recently proposed technique allows one to constrain both the background and perturbation cosmological parameters through the distribution function of supernova Ia apparent magnitudes. Here we extend this technique to alternative cosmological scenarios, in which the growth of structure does not follow the $\Lambda$CDM prescription. We apply the method first to the supernova data provided by the JLA catalog combined with all the current independent redshift distortion data and with low-redshift cluster data from Chandra and show that although the supernovae alone are not very constraining, they help in reducing the confidence regions. Then we apply our method to future data from LSST and from a survey that approximates the Euclid satellite mission. In this case we show that the combined data are nicely complementary and can constrain the normalization $\sigma_8$ and the growth rate index $\gamma$ to within $0.6\%$ and $7\%$, respectively. In particular, the LSST supernova catalog is forecast to give the constraint $\gamma (\sigma_8/0.83)^{6.7} = 0.55 \pm 0.1$. We also report on constraints relative to a step-wise parametrization of the growth rate of structures. These results show that supernova lensing serves as a good cross-check on the measurement of perturbation parameters from more standard techniques.

We perform an analysis of the CMB temperature data taken by the Planck satellite investigating if there is any significant deviation from cosmological isotropy. We look for differences in the spectrum between two opposite hemispheres and also for dipolar modulations. We propose a new way to avoid biases due to partial-sky coverage by producing a mask symmetrized in antipodal directions, in addition to the standard smoothing procedure. We also properly take into account both Doppler and aberration effects due to our peculiar velocity and the anisotropy of the noise, since these effects induce a significant hemispherical asymmetry. We are thus able to probe scales all the way to ell = 2000. After such treatment we find no evidence for significant hemispherical anomalies (i.e. deviations are less than 1.5 sigma when summing over all scales). Although among the larger scales there are sometimes higher discrepancies, these are always less than 3 sigma. We also find results on a dipolar modulation of the power spectrum. Along the hemispheres aligned with the most asymmetric direction for 2 <= ell <= 2000 we find a 3.3 sigma discrepancy when comparing to simulations. However, if we do not restrict ourselves to Planck's maximal asymmetry axis, which can only be known a posteriori, and compare Planck data with the modulation of simulations along their respective maximal asymmetry directions, the discrepancy goes down to less than 1 sigma (with, again, almost 3 sigma discrepancies in some low-ell modes). We thus conclude that no significant power asymmetries seem to be present in the full data set. Interestingly, without proper removal of Doppler and aberration effects one would find spurious anomalies at high ell, between 3 sigma and 5 sigma. Even when considering only ell < 600 we find that the boost is non-negligible and alleviates the discrepancy by roughly half-sigma.

A method was recently proposed which allows the conversion of the weak-lensing effects in the supernova Hubble diagram from noise into signal. Such signal is sensitive to the growth of structure in the universe, and in particular can be used as a measurement of $\sigma_8$ which is independent from more traditional methods such as those based on the CMB, cosmic shear or cluster abundance. We extend here that analysis to allow for intrinsic non-Gaussianities in the supernova PDF, and discuss how this can be best modelled using the Bayes Factor. Although it was shown that a precise measurement of $\sigma_8$ requires ~$10^5$ supernovae, current data already allows an important proof of principle. In particular we make use of the 732 supernovae with z < 1 of the recent JLA catalog and show that a simple treatment of intrinsic non-Gaussianities with a couple of nuisance parameters is enough for our method to yield the values $\sigma_8 = 0.84^{+0.28}_{-0.65}$ or $\sigma_8 < 1.45$ at a $2\sigma$ confidence level. This result is consistent with mock simulations and it is also in agreement with independent measurements and presents the first ever measurement of $\sigma_8$ using supernova magnitudes alone.

Our peculiar motion with respect to the CMB rest frame represents a preferred direction in the observed CMB sky since it induces an apparent deflection of the observed CMB photons (aberration) and a shift in their frequency (Doppler). Both effects distort the multipoles $a_{\ell m}$'s at all $\ell$'s. Such effects are real as it has been recently measured for the first time by Planck according to what was forecast in some recent papers. However, the common lore when estimating a power spectrum from CMB is to consider that Doppler affects only the $\ell=1$ multipole, neglecting any other corrections. In this work we use simulations of the CMB sky in a boosted frame with a peculiar velocity $\beta = v/c = 1.23 \times 10^{-3}$ in order to assess the impact of such effect on power spectrum estimations in different regions of the sky. We show that the boost induces a north-south asymmetry in the power spectrum which is highly significant and non-negligible, of about (0.58 $\pm$ 0.10)% for half-sky cuts when going up to $\ell$ = 2500. We suggest that these effects are relevant and may account for some of the north-south asymmetries seen in the Planck data, being especially important at small scales. Finally we analyze the particular case of the ACT experiment, which observed only a small fraction of the sky and show that it suffers a bias of about 1% on the power spectrum and of similar size on some cosmological parameters: for example the position of the peaks shifts by 0.5% and the overall amplitude of the spectrum is about 0.4% lower than a full-sky case.

Soon the number of type Ia supernova (SN) measurements should exceed 100,000. Understanding the effect of weak lensing by matter structures on the supernova brightness will then be more important than ever. Although SN lensing is usually seen as a source of systematic noise, we will show that it can be in fact turned into signal. More precisely, the non-Gaussianity introduced by lensing in the SN Hubble diagram dispersion depends rather sensitively on the amplitude sigma8 of the matter power spectrum. By exploiting this relation, we are able to predict constraints on sigma8 of 7% (3%) for a catalog of 100,000 (500,000) SNe of average magnitude error 0.12 without having to assume that such intrinsic dispersion is known a priori. The intrinsic dispersion has been assumed to be Gaussian; possible intrinsic non-Gaussianities in the dataset (due to the SN themselves and/or to other transients) could be potentially dealt with by means of additional nuisance parameters describing higher moments of the intrinsic dispersion distribution function. This method is independent of and complementary to the standard methods based on CMB, cosmic shear or cluster abundance observables.

With the availability of thousands of type Ia supernovae in the near future the magnitude scatter induced by lensing will become a major issue as it affects parameter estimation. Current N-body simulations are too time consuming to be integrated in the likelihood analyses used for estimating the cosmological parameters. In this paper we show that in the weak lensing regime a statistical numerical approximation produces accurate results orders of magnitude faster. We write down simple fits to the second, third and fourth central moments of the lensing magnification probability distribution as a function of redshift, of the power spectrum normalization and of the present-day matter density. We also improve upon existing models of lensing variance and show that a shifted lognormal distribution fits well the numerical one. These fits can be easily employed in cosmological likelihood analyses. Moreover, our theoretical predictions make it possible to invert the problem and begin using supernovae lensing to constrain the cosmological parameters.

A great deal of effort is currently being devoted to understanding, estimating and removing systematic errors in cosmological data. In the particular case of type Ia supernovae, systematics are starting to dominate the error budget. Here we propose a Bayesian tool for carrying out a systematic search for systematic contamination. This serves as an extension to the standard goodness-of-fit tests and allows not only to cross-check raw or processed data for the presence of systematics but also to pin-point the data that are most likely contaminated. We successfully test our tool with mock catalogues and conclude that the Union2.1 data do not possess a significant amount of systematics. Finally, we show that if one includes in Union2.1 the supernovae that originally failed the quality cuts, our tool signals the presence of systematics at over 3.8-sigma confidence level.

In recent years the possibility of measuring the temporal change of radial and transverse position of sources in the sky in real time have become conceivable thanks to the thoroughly improved technique applied to new astrometric and spectroscopic experiments, leading to the research domain we call Real-time cosmology. We review for the first time great part of the work done in this field, analysing both the theoretical framework and some endeavor to foresee the observational strategies and their capability to constrain models. We firstly focus on real time measurements of the overall redshift drift and angular separation shift in distant source, able to trace background cosmic expansion and large scale anisotropy, respectively. We then examine the possibility of employing the same kind of observations to probe peculiar and proper acceleration in clustered systems and therefore the gravitational potential. The last two sections are devoted to the short time future change of the cosmic microwave background, as well as to the temporal shift of the temperature anisotropy power spectrum and maps. We conclude revisiting in this context the effort made to forecast the power of upcoming experiments like CODEX, GAIA and PLANCK in providing these new observational tools.

It was recently shown that our peculiar velocity \beta with respect to the CMB induces mixing among multipoles and off-diagonal correlations at all scales which can be used as a measurement of \beta, which is independent of the standard measurement using the CMB temperature dipole. The proposed techniques rely however on a perturbative expansion which breaks down for \ell \gtrsim 1/(\beta) \approx 800. Here we propose a technique which consists of deboosting the CMB temperature in the time-ordered data and show that it extends the validity of the perturbation analysis multipoles up to \ell \sim 10000. We also obtain accurate fitting functions for the mixing between multipoles valid in a full non-linear treatment. Finally we forecast the achievable precision with which these correlations can be measured in a number of current and future CMB missions. We show that Planck could measure the velocity with a precision of around 60 km/s, ACTPol in 4 years around 40 km/s, while proposed future experiments could further shrink this error bar by over a factor of around 2.

Our peculiar velocity with respect to the CMB rest frame is known to induce a large dipole in the CMB. However, the motion of an observer has also the effect of distorting the anisotropies at all scales, as shown by Challinor and Van Leeuwen (2002), due to aberration and Doppler effects. We propose to measure independently our local motion by using off-diagonal two-point correlation functions for high multipoles. We study the observability of the signal for temperature and polarization anisotropies. We point out that Planck can measure the velocity $\beta$ with an error of about 30% and the direction with an error of about 20 degrees. This method constitutes a cross-check, which can be useful to verify that our CMB dipole is due mainly to our velocity or to disentangle the velocity from other possible intrinsic sources. Although in this paper we focus on our peculiar velocity, a similar effect would result also from other intrinsic vectorial distortion of the CMB which would induce a dipolar lensing. Measuring the off-diagonal correlation terms is therefore a test for a preferred direction on the CMB sky.

In addition to shear and vorticity a homogeneous background may also exhibit anisotropic curvature. Here a class of spacetimes is shown to exist where the anisotropy is solely of the latter type, and the shear-free condition is supported by a canonical, massless 2-form field. Such spacetimes possess a preferred direction in the sky and at the same time a CMB which is isotropic at the background level. A distortion of the luminosity distances is derived and used to test the model against the CMB and supernovae (using the Union catalog), and it is concluded that the latter exhibit a higher-than-expected dependence on angular position. It is shown that future surveys could detect a possible preferred direction by observing ~ 20 / (\Omega_{k0}^2) supernovae over the whole sky.

We reanalyze the supernovae data from the Union Compilation including the weak lensing effects caused by inhomogeneities. We compute the lensing probability distribution function for each background solution described by the parameters Omega_M, Omega_L and w in the presence of inhomogeneities, approximately modeled with a single-mass population of halos. We then perform a likelihood analysis in the space of FLRW-parameters and compare our results with the standard approach. We find that the inclusion of lensing can move the best-fit model significantly towards the cosmic concordance of the flat LCDM model, improving the agreement with the constraints coming from the cosmic microwave background and baryon acoustic oscillations.

Two recently proposed techniques, involving the measurement of the cosmic parallax and redshift drift, provide novel ways of directing probing (over a time-span of several years) the background metric of the universe and therefore shed light on the dark energy conundrum. The former makes use of upcoming high-precision astrometry measurements to either observe or put tight constraints on cosmological anisotropy for off-center observers, while the latter employs high-precision spectroscopy to give an independent test of the present acceleration of the universe. In this paper, we show that both methods can break the degeneracy between LTB void models and more traditional dark energy theories. Using the near-future observational missions Gaia and CODEX we show that this distinction might be made with high confidence levels in the course of a decade.

Cosmic parallax is the change of angular separation between pair of sources at cosmological distances induced by an anisotropic expansion. An accurate astrometric experiment like Gaia could observe or put constraints on cosmic parallax. Examples of anisotropic cosmological models are Lemaitre-Tolman-Bondi void models for off-center observers (introduced to explain the observed acceleration without the need for dark energy) and Bianchi metrics. If dark energy has an anisotropic equation of state, as suggested recently, then a substantial anisotropy could arise at $z \lesssim 1$ and escape the stringent constraints from the cosmic microwave background. In this paper we show that such models could be constrained by the Gaia satellite or by an upgraded future mission.

Several authors have argued that self-consistent $f(R)$ gravity models distinct from $\Lambda $CDM are almost ruled out. Confronting such claims, we present a particular two-parameter $f(R)$ model that: (a) is cosmologically viable and distinguishable from $\Lambda $CDM; (b) is compatible with the existence of relativistic stars; (c) is free of singularities of the Ricci scalar during the cosmological evolution and (d) allows the addition of high curvature corrections that could be relevant for inflation.

Refined astrometry measurements allow us to detect large-scale deviations from isotropy through real-time observations of changes in the angular separation between sources at cosmic distances. This "cosmic parallax" effect is a powerful consistency test of FRW metric and may set independent constraints on cosmic anisotropy. We apply this novel general test to LTB cosmologies with off-center observers and show that future satellite missions such as Gaia might achieve accuracies that would put limits on the off-center distance which are competitive with CMB dipole constraints.

Cosmological models involving an interaction between dark matter and dark energy have been proposed in order to solve the so-called coincidence problem. Different forms of coupling have been studied, but there have been claims that observational data seem to narrow (some of) them down to something annoyingly close to the $\Lambda$CDM model, thus greatly reducing their ability to deal with the problem in the first place. The smallness problem of the initial energy density of dark energy has also been a target of cosmological models in recent years. Making use of a moderately general coupling scheme, this paper aims to unite these different approaches and shed some light as to whether this class of models has any true perspective in suppressing the aforementioned issues that plague our current understanding of the universe, in a quantitative and unambiguous way.

A cosmological model that aims at solving the coincidence problem should show that dark energy and dark matter follow the same scaling solution from some time onward. At the same time, the model should contain a sufficiently long matter-dominated epoch that takes place before acceleration in order to guarantee a decelerated epoch and structure formation. So a successful cosmological model requires the occurrence of a sequence of epochs, namely a radiation era, a matter-dominated era and a final accelerated scaling attractor with $\Omega_{\phi} \simeq 0.7$. In this paper we derive the generic form of a scalar-field Lagrangian that possesses scaling solutions in the case where the coupling $Q$ between dark energy and dark matter is a free function of the field $\phi$. We then show, rather surprisingly, that the aforementioned sequence of epochs cannot occur for a vast class of generalized coupled scalar field Lagrangians that includes, to our knowledge, all scaling models in the current literature.