In order to investigate possibilities to measure non-Gaussian signatures of the non-linear iSW effect, we study in this work the family of mixed bispectra <tau^q gamma^(3-q)> and trispectra <tau^q gamma^(4-q)> between the integrated Sachs-Wolfe (iSW) temperature perturbation tau and the galaxy over-density gamma. We use standard Eulerian perturbation theory restricted to tree level expansion for predicting the cosmic matter field. As expected, the spectra are found to decrease in amplitude with increasing q. The transition scale between linear domination and the scales, on which non-linearities take over, moves to larger scales with increasing number of included iSW source fields q. We derive the cumulative signal-to-noise ratios for a combination of Planck CMB data and the galaxy sample of a Euclid-like survey. Including scales down to l_max = 1000 we find sobering values of sigma = 0.83 for the mixed bispectrum and sigma = 0.19 in case of the trispectrum for q=1. For higher values of q the polyspectra <tau^2 gamma> and <tau^3 gamma> are found to be far below the detection limit.

An accurate theoretical description of structure formation at least in the mildly non-linear regime is essential for comparison with data from next generation galaxy surveys. In a recent approach one follows the time evolution of correlators directly and finds a hierarchy of evolution equations with increasing order (Pietroni 2008). So far, in this so called time renormalisation group method the trispectrum was neglected in order to obtain a closed set of equations. In this work we study the influence of the trispectrum on the evolution of the power spectrum. In order to keep the numerical cost at a manageable level we use the tree-level trispectrum from Eulerian perturbation theory. In comparison to numerical simulations we find improvement in the mildly non-linear regime up to k = 0.25 h/Mpc. Beyond k = 0.25 h/Mpc the perturbative description of the trispectrum fails and the method performs worse than without the trispectrum included. Our results reinforce the conceptual advantage of the time renormalisation group method with respect to perturbation theory.

Topic of this article are tomographic measurements of the integrated Sachs-Wolfe effect with specifically designed, orthogonal polynomials which project out statistically independent modes of the galaxy distribution. The polynomials are contructed using the Gram-Schmidt orthogonalisation method. To quantify the power of the iSW-effect in contraining cosmological parameters we perfom a combined Fisher matrix analysis for the iSW-, galaxy- and cross-spectra for wCDM cosmologies using the survey characteristics of PLANCK and EUCLID. The signal to noise ratio has also been studied for other contemporary galaxy surveys, such as SDSS, NVSS and 2MASS. For the cross-spectra our tomographic method provides a 16% increase in the signal to noise ratio and an improvement of up to 30% in conditional errors on parameters. Including all spectra, the marginalised errors approach an inverse square-root dependence with increasing cumulative polynomial order which underlines the statistical independence of the weighted signal spectra.