We study the impact that baryon-CDM relative density perturbations $\delta_{bc}$ have on galaxy formation using cosmological simulations with the IllustrisTNG model. These isocurvature (non-adiabatic) perturbations can be induced primordially, if multiple fields are present during inflation, and are generated before baryon-photon decoupling when baryons did not comove with CDM. The presence of long-wavelength $\delta_{bc}$ perturbations in our simulations is mimicked by modifying the ratios of the cosmic densities of baryons $\Omega_b$ and CDM $\Omega_c$, at fixed total matter density $\Omega_m$. We measure the corresponding galaxy bias parameter $b_{\delta}^{bc}$ as the response of galaxy abundances to $\delta_{bc}$. When selecting by total host halo mass, $b_{\delta}^{bc}$ is negative and it decreases with mass and redshift. Stellar-mass selected simulated galaxies show a weaker or even the opposite trend because of the competing effects of $\delta_{bc}$ on the halo mass function and stellar-to-halo-mass relations. We show that simple modeling of the latter two effects describes $b_{\delta}^{bc}$ for stellar-mass-selected objects well. We find $b_{\delta}^{bc} =0.6$ for $M_* = 10^{11}\ M_{\odot}/h$ and $z=0.5$, which is representative of BOSS DR12 galaxies. For $\delta_{bc}$ modes generated by baryon-photon interactions, we estimate the impact on the DR12 power spectrum to be below $1\%$, and shifts on inferred distance and growth rate parameters should not exceed $0.1\%$.

We study the impact that large-scale perturbations of (i) the matter density and (ii) the primordial gravitational potential with local primordial non-Gaussianity (PNG) have on galaxy formation using the IllustrisTNG model. We focus on the linear galaxy bias $b_1$ and the coefficient $b_\phi$ of the scale-dependent bias induced by PNG, which describe the response of galaxy number counts to these two types of perturbations, respectively. We perform our study using separate universe simulations, in which the effect of the perturbations is mimicked by changes to the cosmological parameters: modified cosmic matter density for $b_1$ and modified amplitude $\mathcal{A}_s$ of the primordial scalar power spectrum for $b_\phi$. We find that the widely used universality relation $b_\phi = 2\delta_c(b_1 - 1)$ is a poor description of the bias of haloes and galaxies selected by stellar mass $M_*$, which is instead described better by $b_\phi(M_*) = 2\delta_c(b_1(M_*) - p)$ with $p \in [0.4, 0.7]$. This is explained by the different impact that matter overdensities and local PNG have on the median stellar-to-halo-mass relation. A simple model of this impact allows us to describe the stellar mass dependence of $b_1$ and $b_\phi$ fairly well. Our results also show a nontrivial relation between $b_1$ and $b_\phi$ for galaxies selected by color and black hole mass accretion rate. Our results provide refined priors on $b_\phi$ for local PNG constraints and forecasts using galaxy clustering. Given that the widely used universality relation underpredicts $b_\phi(M_*)$, existing analyses may underestimate the true constraining power on local PNG.

We derive, using functional methods and the bias expansion, the conditional likelihood for observing a specific tracer field given an underlying matter field. This likelihood is necessary for Bayesian-inference methods. If we neglect all stochastic terms apart from the ones appearing in the auto two-point function of tracers, we recover the result of Schmidt et al., 2018. We then rigorously derive the corrections to this result, such as those coming from a non-Gaussian stochasticity (which include the stochastic corrections to the tracer bispectrum) and higher-derivative terms. We discuss how these corrections can affect current applications of Bayesian inference. We comment on possible extensions to our result, with a particular eye towards primordial non-Gaussianity. This work puts on solid theoretical grounds the EFT-based approach to Bayesian forward modeling.

We present updates on the cosmology inference using the effective field theory (EFT) likelihood presented previously in Schmidt et al., 2018, Elsner et al., 2019 [1,2]. Specifically, we add a cutoff to the initial conditions that serve as starting point for the matter forward model. We show that this cutoff, which was not employed in any previous related work, is important to regularize loop integrals that otherwise involve small-scale, non-perturbative modes. We then present results on the inferred value of the linear power spectrum normalization $\sigma_{8}$ from rest-frame halo catalogs using both second- and third-order bias expansions, imposing uniform priors on all bias parameters. Due to the perfect bias-$\sigma_{8}$ degeneracy at linear order, constraints on $\sigma_{8}$ rely entirely on nonlinear information. The results show the expected convergence behavior when lowering the cutoff in wavenumber, $\Lambda$. When including modes up to $k \leq \Lambda = 0.1\,h\,{\rm Mpc}^{-1}$ in the second-order case, $\sigma_{8}$ is recovered to within $\lesssim 6\,\%$ for a range of halo masses and redshifts. The systematic bias shrinks to $4\,\%$ or less for the third-order bias expansion on the same range of scales. Together with additional evidence we provide, this shows that the residual mismatch in $\sigma_{8}$ can be attributed to higher-order bias contributions. We conclude that the EFT likelihood is able to infer unbiased cosmological constraints, within expected theoretical systematic errors, from physical biased tracers on quasilinear scales

We study graviton non-Gaussianities in the EFT of Inflation. At leading (second) order in derivatives, the graviton bispectrum is fixed by Einstein gravity. There are only two contributions at third order. One of them breaks parity. They come from operators that directly involve the foliation: we then expect sizable non-Gaussianities in three-point functions involving both gravitons and scalars. However, we show that at leading order in slow roll the parity-odd operator does not modify these mixed correlators. We then identify the operators that can affect the graviton bispectrum at fourth order in derivatives. There are two operators that preserve parity. We show that one gives a scalar-tensor-tensor three-point function larger than the one computed in Maldacena, 2003 if $M^2_{\rm P}A_{\rm s}/\Lambda^2\gg 1$ (where $\Lambda$ is the scale suppressing this operator and $A_{\rm s}$ the amplitude of the scalar power spectrum). There are only two parity-odd operators at this order in derivatives.

In the EFT of biased tracers the noise field $\varepsilon_g$ is not exactly uncorrelated with the nonlinear matter field $\delta$. Its correlation with $\delta$ is effectively captured by adding stochasticities to each bias coefficient. We show that if these stochastic fields are Gaussian (the impact of their non-Gaussianity being subleading on quasi-linear scales anyway) it is possible to resum exactly their effect on the conditional likelihood ${\cal P}[\delta_g|\delta]$ to observe a galaxy field $\delta_g$ given an underlying $\delta$. This resummation allows to take them into account in EFT-based approaches to Bayesian forward modeling. We stress that the resulting corrections to a purely Gaussian conditional likelihood with white-noise covariance are the most relevant on scales where the EFT is under control: they are more important than any non-Gaussianity of the noise $\varepsilon_g$.

We investigate the potential of the galaxy power spectrum to constrain compensated isocurvature perturbations (CIPs), primordial fluctuations in the baryon density that are compensated by fluctuations in CDM density to ensure an unperturbed total matter density. We show that CIPs contribute to the galaxy overdensity at linear order, and if they are close to scale-invariant, their effects are nearly perfectly degenerate with the local PNG parameter $f_{\rm nl}$ if they correlate with the adiabatic perturbations. This degeneracy can however be broken by analyzing multiple galaxy samples with different bias parameters, or by taking CMB priors on $f_{\rm nl}$ into account. Parametrizing the amplitude of the CIP power spectrum as $P_{\sigma\sigma} = A_{\rm CIP}^2P_{\mathcal{R}\mathcal{R}}$ (where $P_{\mathcal{R}\mathcal{R}}$ is the adiabatic power spectrum) we find, for a number of fiducial galaxy samples in a simplified forecast setup, that constraints on $A_{\rm CIP}$, relative to those on $f_{\rm nl}$, of order $\sigma_{A_{\rm CIP}}/\sigma_{f_{\rm nl}} \approx 1-2$ are achievable for CIPs correlated with adiabatic perturbations, and $\sigma_{A_{\rm CIP}}/\sigma_{f_{\rm nl}} \approx 5$ for the uncorrelated case. These values are independent of survey volume, and suggest that current galaxy data are already able to improve significantly on the tightest existing constraints on CIPs from the CMB. Future galaxy surveys that aim to achieve $\sigma_{f_{\rm nl}} \sim 1$ have the potential to place even stronger bounds on CIPs.

We investigate the degree to which the Cosmic Microwave Background (CMB) can be used to constrain primordial non-Gaussianity coming from the presence of spinning particles coupled to the inflaton. We compute the $\langle TTT \rangle$ and $\langle TTTT \rangle$ correlation functions arising from the exchange of a particle with spin $s$ and generic mass, and the corresponding signal-to-noise ratios for a cosmic-variance-limited CMB experiment. We show that already with \emph{Planck} data one could improve the theoretical bounds on the amplitude of these primordial templates by an order of magnitude. We particularly emphasize the fact that the trispectrum could be sizable even if the bispectrum is not, making it a prime observable to explore the particle content during inflation.

The fact that the spatial nonlocality of galaxy formation is controlled by some short length scale like the Lagrangian radius is the cornerstone of the bias expansion for large-scale-structure tracers. However, the first sources of ionizing radiation between $z\approx 15$ and $z\approx 6$ are expected to have significant effects on the formation of galaxies we observe at lower redshift, at least on low-mass galaxies. These radiative-transfer effects introduce a new scale in the clustering of galaxies, i.e. the finite distance which ionizing radiation travels until it reaches a given galaxy. This mean free path can be very large, of order $100\,h^{-1}\,{\rm Mpc}$. Consequently, higher-derivative terms in the bias expansion could turn out to be non-negligible even on these scales: treating them perturbatively would lead to a massive loss in predictivity and, for example, could spoil the determination of the BAO feature or constraints on the neutrino mass. Here, we investigate under what assumptions an explicit non-perturbative model of radiative-transfer effects can maintain the robustness of large-scale galaxy clustering as a cosmological probe.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

Distortions of the Cosmic Microwave Background energy spectrum of the $\mu$ type are sensitive to the primordial power spectrum through the dissipation of curvature perturbations on scales $k\simeq 50$ - $10^4\,\mathrm{Mpc}^{-1}$. Their angular correlation with large-scale temperature anisotropies is then sensitive to the squeezed limit of the primordial bispectrum. For inflationary models obeying the single-field consistency relation, we show that the observed $\mu T$ angular correlation that would correspond to the local shape vanishes exactly. All leading non-primordial contributions, including all non-linear production and projection effects, are of the "equilateral shape", namely suppressed by $k^2/{\cal H}_f^2$, where ${\cal H}_f\simeq 10^{-1}\,{\rm Mpc}^{-1}$ is the Hubble radius at the end of the $\mu$-era. Therefore, these non-primordial contributions are orthogonal to a potential local primordial signal (e.g. from multi-field inflation). Moreover, they are very small in amplitude. Our results strengthen the position of $\mu$ distortions as the ultimate probe of local primordial non-Gaussianity.

Long-short mode coupling during inflation, encoded in the squeezed bispectrum of curvature perturbations, induces a dependence of the local, small-scale power spectrum on long-wavelength perturbations, leading to a scale-dependent halo bias. While this scale dependence is absent in the large-scale limit for single-field inflation models that satisfy the consistency relation, certain models such as resonant non-Gaussianity show a peculiar behavior on intermediate scales. We reconsider the predictions for the halo bias in this model by working in Conformal Fermi Coordinates, which isolate the physical effects of long-wavelength perturbations on short-scale physics. We find that the bias oscillates with scale with an envelope similar to that of equilateral non-Gaussianity. Moreover, the bias shows a peculiar modulation with the halo mass. Unfortunately, we find that upcoming surveys will be unable to detect the signal because of its very small amplitude. We also discuss non-Gaussianity due to interactions between the inflaton and massive fields: our results for the bias agree with those in the literature.

We discuss whether an unaccounted contribution to the Cosmic Microwave Background polarization $B$-mode by primordial magnetic fields (PMFs) can bias future constraints on inflationary gravitational waves. As a case-study, we consider a scale-invariant PMF spectrum with amplitude of $\approx{1}\,{\mathrm{nG}}$ on $1\,\mathrm{Mpc}$ scales, compatible with current cosmological bounds. We find a degeneracy in the $B$-mode spectra between PMFs and inflationary gravitational waves. If PMFs of this amplitude are not accounted for, future CMB experiments could claim a false detection of a tensor-to-scalar ratio $r\approx 0.007$, close to the predictions of Starobinsky and $\alpha$-attractor models. The degeneracy can be broken if $B$-modes are measured also at multipoles $\ell\gtrsim 900$: more precisely experiments like CMB-S4 or CORE-M5 would be able to discriminate PMFs from primordial GWs at high statistical significance. Experiments like LiteBIRD or PIXIE will not be able to break the degeneracy and will need complementary bounds coming, for example, from measurements of anisotropies in the Faraday rotation angle of CMB polarization. This reinforces the importance of future experimental constraints on PMFs.

We study generalized disformal transformations, including derivatives of the metric, in the context of the Effective Field Theory of Inflation. All these transformations do not change the late-time cosmological observables but change the coefficients of the operators in the action: some couplings are effectively redundant. At leading order in derivatives and up to cubic order in perturbations, one has 6 free functions that can be used to set to zero 6 of the 17 operators at this order. This is used to show that the tensor three-point function cannot be modified at leading order in derivatives, while the scalar-tensor-tensor correlator can only be modified by changing the scalar dynamics. At higher order in derivatives there are transformations that do not affect the Einstein-Hilbert action: one can find 6 additional transformations that can be used to simplify the inflaton action, at least when the dynamics is dominated by the lowest derivative terms. We also identify the leading higher-derivative corrections to the tensor power spectrum and bispectrum.

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 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)

Gravity is a non-linear theory, and hence, barring cancellations, the initial super-horizon perturbations produced by inflation must contain some minimum amount of mode coupling, or primordial non-Gaussianity. In single-field slow-roll models, where this lower bound is saturated, non-Gaussianity is controlled by two observables: the tensor-to-scalar ratio, which is uncertain by more than fifty orders of magnitude; and the scalar spectral index, or tilt, which is relatively well measured. It is well known that to leading and next-to-leading order in derivatives, the contributions proportional to the tilt disappear from any local observable, and suspicion has been raised that this might happen to all orders, allowing for an arbitrarily low amount of primordial non-Gaussianity. Employing Conformal Fermi Coordinates, we show explicitly that this is not the case. Instead, a contribution of order the tilt appears in local observables. In summary, the floor of physical primordial non-Gaussianity in our Universe has a squeezed-limit scaling of $k_\ell^2/k_s^2$, similar to equilateral and orthogonal shapes, and a dimensionless amplitude of order $0.1\times(n_\mathrm{s}-1)$.

We use the recent observations of Cosmic Microwave Background temperature and polarization anisotropies provided by the Planck satellite experiment to place constraints on the running $\alpha_\mathrm{s} = \mathrm{d}n_{\mathrm{s}} / \mathrm{d}\log k$ and the running of the running $\beta_{\mathrm{s}} = \mathrm{d}\alpha_{\mathrm{s}} / \mathrm{d}\log k$ of the spectral index $n_{\mathrm{s}}$ of primordial scalar fluctuations. We find $\alpha_\mathrm{s}=0.011\pm0.010$ and $\beta_\mathrm{s}=0.027\pm0.013$ at $68\%\,\mathrm{CL}$, suggesting the presence of a running of the running at the level of two standard deviations. We find no significant correlation between $\beta_{\mathrm{s}}$ and foregrounds parameters, with the exception of the point sources amplitude at $143\,\mathrm{GHz}$, $A^{PS}_{143}$, which shifts by half sigma when the running of the running is considered. We further study the cosmological implications of such preference for $\alpha_\mathrm{s},\beta_\mathrm{s}\sim0.01$ by including in the analysis the lensing amplitude $A_L$, the curvature parameter $\Omega_k$, and the sum of neutrino masses $\sum m_{\nu}$. We find that when the running of the running is considered, Planck data are more compatible with the standard expectations of $A_L = 1$ and $\Omega_k = 0$ but still hint at possible deviations. The indication for $\beta_\mathrm{s} > 0$ survives at two standard deviations when external datasets such as BAO and CFHTLenS are included in the analysis, and persists at $\sim 1.7$ standard deviations when CMB lensing is considered. We discuss the possibility of constraining $\beta_\mathrm{s}$ with current and future measurements of CMB spectral distortions, showing that an experiment like PIXIE could provide strong constraints on $\alpha_\mathrm{s}$ and $\beta_\mathrm{s}$.

We discuss the implications of a PIXIE-like experiment, which would measure $\mu$-type spectral distortions of the CMB at a level of $\sigma_{\mu}=(1/n)\times 10^{-8}$, with $n\geq1$ representing an improved sensitivity (e.g. $n=10$ corresponds to PRISM). Using Planck data and considering the six-parameter $\Lambda$CDM model, we compute the posterior for $\mu_8\equiv\mu\times 10^{8}$ and find $\mu_8=1.57^{+0.11}_{-0.13}$ ($68\%\,\mathrm{CL}$). This becomes $\mu_{8} = 1.28^{+0.30}_{-0.52}$ ($68\%\,\mathrm{CL}$) when the running $\alpha_\mathrm{s}$ of the spectral index is included. We point out that a sensitivity of about $3\times$ PIXIE implies a guaranteed discovery: $\mu$-distortion is detected or $\alpha_\mathrm{s}\geq 0$ is excluded (both at $95\%\,\mathrm{CL}$ or higher). This threshold sensitivity sets a clear benchmark for CMB spectrometry. For a combined analysis of PIXIE and current Planck data, we discuss the improvement on measurements of the tilt $n_\mathrm{s}$ and the running $\alpha_\mathrm{s}$ and the dependence on the choice of the pivot. A fiducial running of $\alpha_\mathrm{s}=-0.01$ (close to the Planck best-fit) leads to a detection of negative running at $2\sigma$ for $5\times$ PIXIE. A fiducial running of $\alpha_\mathrm{s}=-0.02$, still compatible with Planck, requires $3\times$ PIXIE to rule out $\alpha_\mathrm{s} = 0$ (at $95\%\,\mathrm{CL}$). We propose a convenient and compact visualization of the improving constraints on the tilt, running and tensor-to-scalar ratio.

Recent cosmic microwave background data in temperature and polarization have reached high precision in estimating all the parameters that describe the current so-called standard cosmological model. Recent results about the integrated Sachs-Wolfe effect from cosmic microwave background anisotropies, galaxy surveys, and their cross-correlations are presented. Looking at fine signatures in the cosmic microwave background, such as the lack of power at low multipoles, the primordial power spectrum and the bounds on non-Gaussianities, complemented by galaxy surveys, we discuss inflationary physics and the generation of primordial perturbations in the early Universe. Three important topics in particle physics, the bounds on neutrinos masses and parameters, on thermal axion mass and on the neutron lifetime derived from cosmological data are reviewed, with attention to the comparison with laboratory experiment results. Recent results from cosmic polarization rotation analyses aimed at testing the Einstein equivalence principle are presented. Finally, we discuss the perspectives of next radio facilities for the improvement of the analysis of future cosmic microwave background spectral distortion experiments.

We present new, tight, constraints on the cosmological background of gravitational waves (GWs) using the latest measurements of CMB temperature and polarization anisotropies provided by the Planck, BICEP2 and Keck Array experiments. These constraints are further improved when the GW contribution $N^{\rm GW}_{\rm eff}$ to the effective number of relativistic degrees of freedom $N_{\rm eff}$ is also considered. Parametrizing the tensor spectrum as a power law with tensor-to-scalar ratio $r$, tilt $n_\mathrm{t}$ and pivot $0.01\,\mathrm{Mpc}^{-1}$, and assuming a minimum value of $r=0.001$, we find $r < 0.089$, $n_\mathrm{t} = 1.7^{+2.1}_{-2.0}$ ($95\%\,\mathrm{CL}$, no $N^{\rm GW}_{\rm eff}$) and $r < 0.082$, $n_\mathrm{t} = -0.05^{+0.58}_{-0.87}$ ($95\%\,\mathrm{CL}$, with $N^{\rm GW}_{\rm eff}$). When the recently released $95\,\mathrm{GHz}$ data from Keck Array are added to the analysis, the constraints on $r$ are improved to $r < 0.067$ ($95\%\,\mathrm{CL}$, no $N^{\rm GW}_{\rm eff}$), $r < 0.061$ ($95\%\,\mathrm{CL}$, with $N^{\rm GW}_{\rm eff}$). We discuss the limits coming from direct detection experiments such as LIGO-Virgo, pulsar timing (European Pulsar Timing Array) and CMB spectral distortions (FIRAS). Finally, we show future constraints achievable from a COrE-like mission: if the tensor-to-scalar ratio is of order $10^{-2}$ and the inflationary consistency relation $n_\mathrm{t} = -r/8$ holds, COrE will be able to constrain $n_\mathrm{t}$ with an error of $0.16$ at $95\%\,\mathrm{CL}$. In the case that lensing $B$-modes can be subtracted to $10\%$ of their power, a feasible goal for COrE, these limits will be improved to $0.11$ at $95\%\,\mathrm{CL}$.

The Integrated Sachs-Wolfe (ISW) effect predicts additional anisotropies in the Cosmic Microwave Background due to time variation of the gravitational potential when the expansion of the universe is not matter dominated. The ISW effect is therefore expected in the early universe, due to the presence of relativistic particles at recombination, and in the late universe, when dark energy starts to dominate the expansion. Deviations from the standard picture can be parameterized by $A_{e\text{ISW}}$ and $A_{l\text{ISW}}$, which rescale the overall amplitude of the early and late ISW effects. Analyzing the most recent CMB temperature spectra from the Planck 2015 release, we detect the presence of the early ISW at high significance with $A_{e\mathrm{ISW}} = 1.06\pm0.04$ at 68% CL and an upper limit for the late ISW of $A_{l\mathrm{ISW}} < 1.1$ at 95% CL. The inclusion of the recent polarization data from the Planck experiment results in $A_{e\mathrm{ISW}} = 0.999\pm0.028$ at 68% CL, in better agreement with the value $A_{e\mathrm{ISW}} = 1$ of a standard cosmology. When considering the recent detections of the late ISW coming from correlations between CMB temperature anisotropies and weak lensing, a value of $A_{l\mathrm{ISW}}=0.85\pm0.21$ is predicted at 68% CL, showing a $4\sigma$ evidence. We discuss the stability of our result in the case of an extra relativistic energy component parametrized by the effective neutrino number $N_\mathrm{eff}$ and of a CMB lensing amplitude $A_\mathrm{L}$.