The subject of this paper is a quantification of the information content of cosmological probes of the large-scale structures, specifically of temperature and polarisation anisotropies in the cosmic microwave background, CMB-lensing, weak cosmic shear and galaxy clustering, in terms of Information theory measures like information entropies. We aim to establish relationships for Gaussian likelihoods, between conventional measures of statistical uncertainties and information entropies. Furthermore, we extend these studies to the computation of (Bayesian) evidences and the power of measurement to distinguish between competing models. We investigate in detail how cosmological data decreases information entropy by reducing statistical errors and by breaking degeneracies. In addition, we work out how tensions between data sets increase information entropy and quantify this effect in three examples: the discrepancy in $\Omega_m$ and $\sigma_8$ between the CMB and weak lensing, the role of intrinsic alignments in weak lensing data when attempting the dark energy equation of state parameters, and the famous $H_0$-tension between Cepheids in the Hubble keystone project and the cosmic microwave background as observed by Planck.

Statistical inference more often than not involves models which are non-linear in the parameters thus leading to non-Gaussian posteriors. Many computational and analytical tools exist that can deal with non-Gaussian distributions, and empirical Gaussianisation transforms can reduce the amount of non-Gaussianity in a distribution. Alternatively, in this work, we employ methods from information geometry. The latter formulates a set of probability distributions for some given model as a manifold employing a Riemannian structure, equipped with a metric, the Fisher information. In this framework we study the differential geometrical meaning of non-Gaussianities in a higher order Fisher approximation, and their respective transformation behaviour under re-parameterisation, which corresponds to a chart transition on the statistical manifold. While weak non-Gaussianities vanish in normal coordinates in a first order approximation, one can in general not find transformations that discard non-Gaussianities globally. As an application we consider the likelihood of the supernovae distance-redshift relation in cosmology for the parameter pair ($\Omega_{\mathrm{m_0}}$, $w$). We demonstrate the connection between confidence intervals and geodesic length and demonstrate how the Lie-derivative along the degeneracy directions gives hints at possible isometries of the Fisher metric.

Atomic nuclei appearing in cosmic rays are typically classified as primary or secondary. However, a better understanding of their origin and propagation properties is still necessary. We analyse the flux of primary (He, C, O) and secondary nuclei (Li, Be, B) detected with rigidity (momentum/charge) between 2 GV and 3 TV by the Alpha Magnetic Spectrometer (AMS) on the International Space Station. We show that $q$-exponential distribution functions, as motivated by generalized versions of statistical mechanics with temperature fluctuations, provide excellent fits for the measured flux of all nuclei considered. Primary and secondary fluxes reveal a universal dependence on kinetic energy per nucleon for which the underlying energy distribution functions are solely distinguished by their effective degrees of freedom. All given spectra are characterized by a universal mean temperature parameter $\sim$ 200 MeV which agrees with the Hagedorn temperature. Our analysis suggests that QCD scattering processes together with nonequilibrium temperature fluctuations imprint heavily and universally onto the measured cosmic ray spectra, and produce a similar shape of energy spectra as high energy collider experiments on the Earth.

We present constraints on Horndeski gravity from a combined analysis of cosmic shear, galaxy-galaxy lensing and galaxy clustering from $450\,\mathrm{deg}^2$ of the Kilo-Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey. The Horndeski class of dark energy/modified gravity models includes the majority of universally coupled extensions to $\Lambda$CDM with one scalar field in addition to the metric. We study the functions of time that fully describe the evolution of linear perturbations in Horndeski gravity. Our results are compatible throughout with a $\Lambda$CDM model. By imposing gravitational wave constraints, we fix the tensor speed excess to zero and consider a subset of models including e.g. quintessence and $f(R)$ theories. Assuming proportionality of the Horndeski functions $\alpha_B$ and $\alpha_M$ (kinetic braiding and the Planck mass run rate, respectively) to the dark energy density fraction $\Omega_{\mathrm{DE}}(a) = 1 - \Omega_{\mathrm{m}}(a)$, we find for the proportionality coefficients $\hat{\alpha}_B = 0.20_{-0.33}^{+0.20} \,$ and $\, \hat{\alpha}_M = 0.25_{-0.29}^{+0.19}$. Our value of $S_8 \equiv \sigma_8 \sqrt{\Omega_{\mathrm{m}}/0.3}$ is in better agreement with the $Planck$ estimate when measured in the enlarged Horndeski parameter space than in a pure $\Lambda$CDM scenario. In our joint three-probe analysis we report a downward shift of the $S_8$ best fit value from the $Planck$ measurement of $\Delta S_8 = 0.016_{-0.046}^{+0.048}$ in Horndeski gravity, compared to $\Delta S_8 = 0.059_{-0.039}^{+0.040}$ in $\Lambda$CDM. Our constraints are robust to the modelling uncertainty of the non-linear matter power spectrum in Horndeski gravity. Our likelihood code for multi-probe analysis in both $\Lambda$CDM and Horndeski gravity is publicly available at http://github.com/alessiospuriomancini/KiDSHorndeski .

We study the parameter estimation bias induced by intrinsic alignments on a Euclid-like weak lensing survey. For the intrinsic alignment signal we assume a composite alignment model for elliptical and spiral galaxies using tidal shearing and tidal torquing as the alignment generating mechanism, respectively. The parameter estimation bias is carried out analytically with a Gaussian bias model and through running Monte-Carlo-Markov-chains on synthetic data including the alignment signal with a likelihood only including the cosmic shear signal. In particular, we study the impact of $II$ and $GI$ alignment terms individually as well as the more realistic situation where both types of alignment are present, and investigate the scaling of the estimation biases with varying strength of the alignment signal. Our results show that intrinsic alignments can cause substantial biases in cosmological parameters even if the coupling of galaxies to the ambient large is small. Especially $GI$-contributions strongly bias key cosmological parameters such as the dark energy equation of state. We also correct the analytic expression for the Gaussian bias model and find that the biases induced by intrinsic alignments are not accurately recovered by the simple analytic model.

Observations of the large-scale structure (LSS) implicitly assume an ideal FLRW observer with the ambient structure having no influence on the observer. However, due to correlations in the LSS, cosmological observables are dependent on the position of an observer. We investigate this influence in full generality for a weakly non-Gaussian random field, for which we derive expressions for angular spectra of large-scale structure observables conditional on a property of the large-scale structure that is typical for the observer's location. As an application, we then apply to the formalism to angular spectra of the weak gravitational lensing effect and provide numerical estimates for the resulting change on the spectra using linear structure formation. For angular weak lensing spectra we find the effect to be of order of a few percent, for instance we estimate for an overdensity of $\delta =0.5$ and multipoles up to $\ell=100$ the change in the weak lensing spectra to be approximately 4 percent. We show that without accounting for correlation between the density at observer's location and the weak gravitational lensing spectra, the values of the parameters $\Omega_m$ and $\sigma_8$ are underestimated by a few percent. Thus, this effect will be important when analysing data from future surveys such as Euclid, which aim at the percent-level precision. The effect is difficult to capture in simulations, as estimates of the number of numerical simulations necessary to quantify the effect are high.

Future observations of the large-scale structure have the potential to investigate cosmological models with a high degree of complexity, including the properties of gravity on large scales, the presence of a complicated dark energy component, and the addition of neutrinos changing structures on small scales. Here we study Horndeski theories of gravity, the most general minimally coupled scalar-tensor theories of second order. While the cosmological background evolution can be described by an effective equation of state, the perturbations are characterised by four free functions of time. We consider a specific parametrisation of these functions tracing the dark energy component. The likelihood of the full parameter set resulting from combining cosmic microwave background primary anisotropies including their gravitational lensing signal, tomographic angular galaxy clustering and weak cosmic shear, together with all possible non-vanishing cross-correlations is evaluated; both with the Fisher-formalism as well as without the assumption of a specific functional form of the posterior through Monte-Carlo Markov-chains (MCMCs). Our results show that even complex cosmological models can be constrained and could exclude variations of the effective Newtonian gravitational coupling larger than 10\% over the age of the Universe. In particular, we confirm strong correlations between parameter groups. Furthermore, we find that the expected contours from MCMC are significantly larger than those from the Fisher analysis even with the vast amount of signal provided by stage IV experiments, illustrating the importance of a proper treatment of non-Gaussian likelihoods and the high level of precision needed for unlocking the sensitivity on gravitational parameters.

In this work we investigate the environmental dependence of the intrinsic ellipticity correlations in cosmic shear surveys. We use the quadratic and linear alignment model to describe the contributions by spiral and elliptical galaxies, respectively. The density field is in both cases described by a Gaussian random field and ellipticity correlation functions that are conditional on the environment of the galaxies are constructed by sampling random values for the tidal tensor and inertial tensor. The covariance of the Gaussian random process from which the tensor entries are drawn is decomposed by means of a spherical Fourier-Bessel transformation of the density field. The dependence on environment is modelled by the number of positive eigenvalues of the tidal tensor, which allows a differentiation between voids, sheets, filaments and superclusters. We find that elliptical galaxies align strongest in elongated structures such as sheets and filaments with an amplitude almost an order of magnitude higher compared to the alignment in clusters or voids. In contrast to this, spiral galaxies align equally strong in all environments. Cross-alignments between different environments are smaller than the respective auto-correlations subject to the Cauchy-Schwarz inequality which is an effective bound on their amplitude. Furthermore, we find misalignment between inertial and tidal tensor to be stronger in anisotropic regions compared to clusters or voids. While the imprint of weak lensing on galaxy ellipticities is agnostic about the environment, using environment information can help to distinguish between lensing and the intrinsic alignment signal.

Cosmic shear is one of the primary probes to test gravity with current and future surveys. There are two main techniques to analyse a cosmic shear survey; a tomographic method, where correlations between the lensing signal in different redshift bins are used to recover redshift information, and a 3D approach, where the full redshift information is carried through the entire analysis. Here we compare the two methods, by forecasting cosmological constraints for future surveys like Euclid. We extend the 3D formalism for the first time to theories beyond the standard model, belonging to the Horndeski class. This includes the majority of universally coupled extensions to $\Lambda$CDM with one scalar degree of freedom in addition to the metric, still in agreement with current observations. Given a fixed background, the evolution of linear perturbations in Horndeski gravity is described by a set of four functions of time only. We model their time evolution assuming proportionality to the dark energy density fraction and place Fisher matrix constraints on the proportionality coefficients. We find that a 3D analysis can constrain Horndeski theories better than a tomographic one, in particular with a decrease in the errors of the order of 20$\%$. This paper shows for the first time a quantitative comparison on an equal footing between Fisher matrix forecasts for both a fully 3D and a tomographic analysis of cosmic shear surveys. The increased sensitivity of the 3D formalism comes from its ability to retain information on the source redshifts along the entire analysis.

Cosmic shear - the weak gravitational lensing effect generated by fluctuations of the gravitational tidal fields of the large-scale structure - is one of the most promising tools for current and future cosmological analyses. The spherical-Bessel decomposition of the cosmic shear field ("3D cosmic shear") is one way to maximise the amount of redshift information in a lensing analysis and therefore provides a powerful tool to investigate in particular the growth of cosmic structure that is crucial for dark energy studies. However, the computation of simulated 3D cosmic shear covariance matrices presents numerical difficulties, due to the required integrations over highly oscillatory functions. We present and compare two numerical methods and relative implementations to perform these integrations. We then show how to generate 3D Gaussian random fields on the sky in spherical coordinates, starting from the 3D cosmic shear covariances. To validate our field-generation procedure, we calculate the Minkowski functionals associated with our random fields, compare them with the known expectation values for the Gaussian case and demonstrate parameter inference from Minkowski functionals from a cosmic shear survey. This is a first step towards producing fully 3D Minkowski functionals for a lognormal field in 3D to extract Gaussian and non-Gaussian information from the cosmic shear field, as well as towards the use of Minkowski functionals as a probe of cosmology beyond the commonly used two-point statistics.

The Heraeus Summer School series "Astronomy from four perspectives", funded by the WE Heraeus Foundation, draws together teachers and teacher students, astronomers, physicists and astronomy students from Germany and Italy. For each summer school, participants gather at one of the four participating nodes: Heidelberg, Padua, Jena, and Florence. The main goal of the series is to bring astronomy into schools, which is achieved by educating and training the teachers and teacher students. In this e-print, we present the exercises, tutorials, and high-school classroom materials developed during the fifth summer school of the series, which took place at Haus der Astronomie in Heidelberg August 26 -- September 2, 2017. "The tutorials" were prepared beforehand for the participants of the Summer schools, and are suitable for use in teacher training. "Classroom materials" were developed mainly during the summer school itself, and are suitable for high-school level teaching. They include question sheets for pupils, and some pointers on where to use the material in the German high school curriculum. Both sets of materials address the summer school's four main topics: Supernova cosmology, the virial theorem, rotation curves of galaxies, and the temperature of the cosmic microwave background (CMB).

Traditionally the spherical collapse of objects is studied with respect to a uniform background density, yielding the critical over-density $\delta_\mathrm{c}$ as key ingredient to the mass function of virialized objects. Here we investigate the shear and rotation acting on a peak in a Gaussian random field. By assuming that collapsing objects mainly form at those peaks, we use this shear and rotation as external effects changing the dynamics of the spherical collapse, which is described by the Raychaudhuri equation. We therefore assume that the shear and rotation have no additional dynamics on top of their cosmological evolution and thus only appear as inhomogeneities in the differential equation.

Tidal gravitational forces can modify the shape of galaxies and clusters of galaxies, thus correlating their orientation with the surrounding matter density field. We study the dependence of this phenomenon, known as intrinsic alignment (IA), on the mass of the dark matter haloes that host these bright structures, analysing the Millennium and Millennium-XXL $N$-body simulations. We closely follow the observational approach, measuring the halo position-halo shape alignment and subsequently dividing out the dependence on halo bias. We derive a theoretical scaling of the IA amplitude with mass in a dark matter universe, and predict a power-law with slope $\beta_{\mathrm{M}}$ in the range $1/3$ to $1/2$, depending on mass scale. We find that the simulation data agree with each other and with the theoretical prediction remarkably well over three orders of magnitude in mass, with the joint analysis yielding an estimate of $\beta_{\mathrm{M}} = 0.36^{+0.01}_{-0.01}$. This result does not depend on redshift or on the details of the halo shape measurement. The analysis is repeated on observational data, obtaining a significantly higher value, $\beta_{\mathrm{M}} = 0.56^{+0.05}_{-0.05}$. There are also small but significant deviations from our simple model in the simulation signals at both the high- and low-mass end. We discuss possible reasons for these discrepancies, and argue that they can be attributed to physical processes not captured in the model or in the dark matter-only simulations.

We extend a previous work by Reischke et al., 2016 by studying the effects of tidal shear on clustering dark energy models within the framework of the extended spherical collapse model and using the Zel'dovich approximation. As in previous works on clustering dark energy, we assumed a vanishing effective sound speed describing the perturbations in dark energy models. To be self-consistent, our treatment is valid only on linear scales since we do not intend to introduce any heuristic models. This approach makes the linear overdensity $\delta_{\rm c}$ mass dependent and similarly to the case of smooth dark energy, its effects are predominant at small masses and redshifts. Tidal shear has effects of the order of percent or less, regardless of the model and preserves a well known feature of clustering dark energy: When dark energy perturbations are included, the models resemble better the $\Lambda$CDM evolution of perturbations. We also showed that effects on the comoving number density of halos are small and qualitatively and quantitatively in agreement with what previously found for smooth dark energy models.

The likelihood function for cosmological parameters, given by e.g. weak lensing shear measurements, depends on contributions to the covariance induced by the nonlinear evolution of the cosmic web. As nonlinear clustering to date has only been described by numerical $N$-body simulations in a reliable and sufficiently precise way, the necessary computational costs for estimating those covariances at different points in parameter space are tremendous. In this work we describe the change of the matter covariance and of the weak lensing covariance matrix as a function of cosmological parameters by constructing a suitable basis, where we model the contribution to the covariance from nonlinear structure formation using Eulerian perturbation theory at third order. We show that our formalism is capable of dealing with large matrices and reproduces expected degeneracies and scaling with cosmological parameters in a reliable way. Comparing our analytical results to numerical simulations we find that the method describes the variation of the covariance matrix found in the SUNGLASS weak lensing simulation pipeline within the errors at one-loop and tree-level for the spectrum and the trispectrum, respectively, for multipoles up to $\ell\leq 1300$. We show that it is possible to optimize the sampling of parameter space where numerical simulations should be carried out by minimising interpolation errors and propose a corresponding method to distribute points in parameter space in an economical way.

Forecasts in cosmology, both with Monte-Carlo Markov-chain methods and with the Fisher matrix formalism, depend on the choice of the fiducial model because both the signal strength of any observable as well as the model nonlinearities linking observables to cosmological parameters vary in the general case. In this paper we propose a method for extrapolating Fisher-forecasts across the space of cosmological parameters by constructing a suitable ba- sis. We demonstrate the validity of our method with constraints on a standard dark energy model extrapolated from a {\Lambda}CDM-model, as can be expected from 2-bin weak lensing to- mography with a Euclid-like survey, in the parameter pairs $(\Omega_\text{m},\sigma_8)$, $(\Omega_\text{m}, w_0)$ and $(w_0, w_\text{a})$. Our numerical results include very accurate extrapolations across a wide range of cosmo- logical parameters in terms of shape, size and orientation of the parameter likelihood, and a decomposition of the change of the likelihood contours into modes, which are straightforward to interpret in a geometrical way. We find that in particular the variation of the dark energy figure of merit is well captured by our formalism.

We study the spherical collapse model in the presence of external gravitational tidal shear fields for different dark energy scenarios and investigate the impact on the mass function and cluster number counts. While previous studies of the influence of shear and rotation on $\delta_\mathrm{c}$ have been performed with heuristically motivated models, we try to avoid this model dependence and sample the external tidal shear values directly from the statistics of the underlying linearly evolved density field based on first order Lagrangian perturbation theory. Within this self-consistent approach, in the sense that we restrict our treatment to scales where linear theory is still applicable, only fluctuations larger than the scale of the considered objects are included into the sampling process which naturally introduces a mass dependence of $\delta_\mathrm{c}$. We find that shear effects are predominant for smaller objects and at lower redshifts, i. e. the effect on $\delta_\mathrm{c}$ is at or below the percent level for the $\Lambda$CDM model. For dark energy models we also find small but noticeable differences, similar to $\Lambda$CDM. The virial overdensity $\Delta_\mathrm{V}$ is nearly unaffected by the external shear. The now mass dependent $\delta_c$ is used to evaluate the mass function for different dark energy scenarios and afterwards to predict cluster number counts, which indicate that ignoring the shear contribution can lead to biases of the order of $1\sigma$ in the estimation of cosmological parameters like $\Omega_\mathrm{m}$, $\sigma_8$ or $w$.

Galaxy shapes are not randomly oriented, rather they are statistically aligned in a way that can depend on formation environment, history and galaxy type. Studying the alignment of galaxies can therefore deliver important information about the physics of galaxy formation and evolution as well as the growth of structure in the Universe. In this review paper we summarise key measurements of galaxy alignments, divided by galaxy type, scale and environment. We also cover the statistics and formalism necessary to understand the observations in the literature. With the emergence of weak gravitational lensing as a precision probe of cosmology, galaxy alignments have taken on an added importance because they can mimic cosmic shear, the effect of gravitational lensing by large-scale structure on observed galaxy shapes. This makes galaxy alignments, commonly referred to as intrinsic alignments, an important systematic effect in weak lensing studies. We quantify the impact of intrinsic alignments on cosmic shear surveys and finish by reviewing practical mitigation techniques which attempt to remove contamination by intrinsic alignments.

The alignments between galaxies, their underlying matter structures, and the cosmic web constitute vital ingredients for a comprehensive understanding of gravity, the nature of matter, and structure formation in the Universe. We provide an overview on the state of the art in the study of these alignment processes and their observational signatures, aimed at a non-specialist audience. The development of the field over the past one hundred years is briefly reviewed. We also discuss the impact of galaxy alignments on measurements of weak gravitational lensing, and discuss avenues for making theoretical and observational progress over the coming decade.

The shapes of galaxies are not randomly oriented on the sky. During the galaxy formation and evolution process, environment has a strong influence, as tidal gravitational fields in the large-scale structure tend to align nearby galaxies. Additionally, events such as galaxy mergers affect the relative alignments of both the shapes and angular momenta of galaxies throughout their history. These "intrinsic galaxy alignments" are known to exist, but are still poorly understood. This review will offer a pedagogical introduction to the current theories that describe intrinsic galaxy alignments, including the apparent difference in intrinsic alignment between early- and late-type galaxies and the latest efforts to model them analytically. It will then describe the ongoing efforts to simulate intrinsic alignments using both N-body and hydrodynamic simulations. Due to the relative youth of this field, there is still much to be done to understand intrinsic galaxy alignments and this review summarises the current state of the field, providing a solid basis for future work.

Assuming a Euclid-like weak lensing data set, we compare different methods of dealing with its inherent parameter degeneracies. Including priors into a data analysis can mask the information content of a given data set alone. However, since the information content of a data set is usually estimated with the Fisher matrix, priors are added in order to enforce an approximately Gaussian likelihood. Here, we compare priorless forecasts to more conventional forecasts that use priors. We find strongly non-Gaussian likelihoods for 2d-weak lensing if no priors are used, which we approximate with the DALI-expansion. Without priors, the Fisher matrix of the 2d-weak lensing likelihood includes unphysical values of $\Omega_m$ and $h$, since it does not capture the shape of the likelihood well. The Cramer-Rao inequality then does not need to apply. We find that DALI and Monte Carlo Markov Chains predict the presence of a dark energy with high significance, whereas a Fisher forecast of the same data set also allows decelerated expansion. We also find that a 2d-weak lensing analysis provides a sharp lower limit on the Hubble constant of $h > 0.4$, even if the equation of state of dark energy is jointly constrained by the data. This is not predicted by the Fisher matrix and usually masked in other works by a sharp prior on $h$. Additionally, we find that DALI estimates Figures of Merit in the presence of non-Gaussianities better than the Fisher matrix. We additionally demonstrate how DALI allows switching to a Hamiltonian Monte Carlo sampling of a highly curved likelihood with acceptance rates of $\approx 0.5$, an effective covering of the parameter space, and numerically effectively costless leapfrog steps. This shows how quick forecasts can be upgraded to accurate forecasts whenever needed. Results were gained with the public code from http://lnasellentin.github.io/DALI/

Many weak lensing calculations make use of the Born approximation where the light ray is approximated by a straight path. We examine the effect of Born-corrections for lensing of the cosmic microwave background in an analytical approach by taking perturbative corrections to the geodesic into account. The resulting extra power in the lensing potential spectrum is comparable to the power generated by nonlinear structure formation and affects especially the polarisation spectra, leading to relative changes of the order of one per cent for the $E$-mode spectrum and up to 10 per cent on all scales to the $B$-mode spectrum. In contrast, there is only little change of spectra involving the CMB temperature. Additionally, the corrections excite one more degree of freedom resulting in a deflection component which can not be described as a gradient of the lensing potential as it is related to image rotation in lens-lens coupling. We estimate the magnitude of this effect on the CMB-spectra and find it to be negligible.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

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.

We introduce a novel method for weak-lensing measurements, which is based on a mathematically exact deconvolution of the moments of the apparent brightness distribution of galaxies from the telescope's PSF. No assumptions on the shape of the galaxy or the PSF are made. The (de)convolution equations are exact for unweighted moments only, while in practice a compact weight function needs to be applied to the noisy images to ensure that the moment measurement yields significant results. We employ a Gaussian weight function, whose centroid and ellipticity are iteratively adjusted to match the corresponding quantities of the source. The change of the moments caused by the application of the weight function can then be corrected by considering higher-order weighted moments of the same source. Because of the form of the deconvolution equations, even an incomplete weighting correction leads to an excellent shear estimation if galaxies and PSF are measured with a weight function of identical size. We demonstrate the accuracy and capabilities of this new method in the context of weak gravitational lensing measurements with a set of specialized tests and show its competitive performance on the GREAT08 challenge data. A complete C++ implementation of the method can be requested from the authors.

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.

As recently shown by Viola et al., the common (KSB) method for measuring weak gravitational shear creates a non-linear relation between the measured and the true shear of objects. We investigate here what effect such a non-linear calibration relation may have on cosmological parameter estimates from weak lensing if a simpler, linear calibration relation is assumed. We show that the non-linear relation introduces a bias in the shear-correlation amplitude and thus a bias in the cosmological parameters Omega_matter and sigma_8. Its direction and magnitude depends on whether the point-spread function is narrow or wide compared to the galaxy images from which the shear is measured. Substantial over- or underestimates of the cosmological parameters are equally possible, depending also on the variant of the KSB method. Our results show that for trustable cosmological-parameter estimates from measurements of weak lensing, one must verify that the method employed is free from ellipticity-dependent biases or monitor that the calibration relation inferred from simulations is applicable to the survey at hand.

In this paper I investigate the family of mixed three-point correlation functions <tau^q gamma^{3-q}>, q=0,1,2, between the integrated Sachs-Wolfe (iSW) temperature perturbation tau and the galaxy overdensity gamma as a tool for detecting the gravitational interaction of cosmic microwave background (CMB) photons with the potentials of nonlinearly evolving cosmological structures. Both the iSW-effect as well as the galaxy overdensity are derived in hyper-extended perturbation theory to second order and I emphasise the different parameter sensitivities of the linear and non-linear iSW-effect. I examine the configuration dependence of the relevant bispectra, quantify their sensitivities and discuss their degeneracies with respect to the cosmological parameters Omega_m, sigma_8, h and the dark energy equation of state parameter $w$. I give detection sigificances for combining PLANCK CMB data and the DUNE galaxy sample, by using a quadratic approximation for the likelihood with LambdaCDM as the fiducial cosmology: The combination of PLANCK with DUNE should be able to reach a cumulative signal to noise ratio of ~0.6sigma for the bispectrum <tau gamma^2> up to l=2000, where the most important noise source are the primary CMB fluctuations.