We reconstruct the 3D structure of magnetic fields, which were seeded by density perturbations during the radiation dominated epoch of the Universe and later on were evolved by structure formation. To achieve this goal, we rely on three dimensional initial density fields inferred from the 2M++ galaxy compilation via the Bayesian $\texttt{BORG}$ algorithm. Using those, we estimate the magnetogenesis by the so called Harrison mechanism. This effect produced magnetic fields exploiting the different photon drag on electrons and ions in vortical motions, which are exited due to second order perturbation effects in the Early Universe. Subsequently we study the evolution of these seed fields through the non-linear cosmic structure formation by virtue of a MHD simulation to obtain a 3D estimate for the structure of this primordial magnetic field component today. At recombination we obtain a reliable lower limit on the large scale magnetic field strength around $10^{-23} \mathrm{G}$, with a power spectrum peaking at about $ 2\, \mathrm{Mpc}^{-1}h$ in comoving scales. At present we expect this evolved primordial field to have strengthts above $\approx 10^{-27}\, \mathrm{G}$ and $\approx 10^{-29}\, \mathrm{G}$ in clusters of galaxies and voids, respectively. We also calculate the corresponding Faraday rotation measure map and show the magnetic field morphology and strength for specific objects of the Local Universe.

We present charm (cosmic history agnostic reconstruction method), a novel inference algorithm that reconstructs the cosmic expansion history as encoded in the Hubble parameter $H(z)$ from SNe Ia data. The novelty of the approach lies in the usage of information field theory, a statistical field theory that is very well suited for the construction of optimal signal recovery algorithms. The charm algorithm infers non-parametrically $s(a)=\ln(\rho(a)/\rho_{\mathrm{crit}0})$, the density evolution which determines $H(z)$, without assuming an analytical form of $\rho(a)$ but only its smoothness with the scale factor $a=(1+z)^{-1}$. The inference problem of recovering the signal $s(a)$ from the data is formulated in a fully Bayesian way. In detail, we have rewritten the signal as the sum of a background cosmology and a perturbation. This allows us to determine the maximum a posteriory estimate of the signal by an iterative Wiener filter method. Applying charm to the Union2.1 supernova compilation, we have recovered a cosmic expansion history that is fully compatible with the standard $\Lambda$CDM cosmological expansion history with parameter values consistent with the results of the Planck mission.

Matrix determinants play an important role in data analysis, in particular when Gaussian processes are involved. Due to currently exploding data volumes, linear operations - matrices - acting on the data are often not accessible directly but are only represented indirectly in form of a computer routine. Such a routine implements the transformation a data vector undergoes under matrix multiplication. While efficient probing routines to estimate a matrix's diagonal or trace, based solely on such computationally affordable matrix-vector multiplications, are well known and frequently used in signal inference, there is no stochastic estimate for its determinant. We introduce a probing method for the logarithm of a determinant of a linear operator. Our method rests upon a reformulation of the log-determinant by an integral representation and the transformation of the involved terms into stochastic expressions. This stochastic determinant determination enables large-size applications in Bayesian inference, in particular evidence calculations, model comparison, and posterior determination.

The calibration of a measurement device is crucial for every scientific experiment, where a signal has to be inferred from data. We present CURE, the calibration uncertainty renormalized estimator, to reconstruct a signal and simultaneously the instrument's calibration from the same data without knowing the exact calibration, but its covariance structure. The idea of CURE, developed in the framework of information field theory, is starting with an assumed calibration to successively include more and more portions of calibration uncertainty into the signal inference equations and to absorb the resulting corrections into renormalized signal (and calibration) solutions. Thereby, the signal inference and calibration problem turns into solving a single system of ordinary differential equations and can be identified with common resummation techniques used in field theories. We verify CURE by applying it to a simplistic toy example and compare it against existent self-calibration schemes, Wiener filter solutions, and Markov Chain Monte Carlo sampling. We conclude that the method is able to keep up in accuracy with the best self-calibration methods and serves as a non-iterative alternative to it.

An essential quantity required to understand the physics of the early Universe, in particular the inflationary epoch, is the primordial scalar potential $\Phi$ and its statistics. We present for the first time an all-sky reconstruction of $\Phi$ with corresponding $1\sigma$-uncertainty from WMAP's cosmic microwave background (CMB) temperature data -- a map of the very early Universe right after the inflationary epoch. This has been achieved by applying a Bayesian inference method that separates the whole inverse problem of the reconstruction into many independent ones, each of them solved by an optimal linear filter (Wiener filter). In this way, the three-dimensional potential $\Phi$ gets reconstructed slice by slice resulting in a thick shell of nested spheres around the comoving distance to the last scattering surface. Each slice represents the primordial scalar potential $\Phi$ projected onto a sphere with corresponding distance. Furthermore, we present an advanced method for inferring $\Phi$ and its power spectrum simultaneously from data, but argue that applying it requires polarization data with high signal-to-noise levels not available yet. Future CMB data should improve results significantly, as polarization data will fill the present $\ell-$blind gaps of the reconstruction.

The presence of multiple fields during inflation might seed a detectable amount of non-Gaussianity in the curvature perturbations, which in turn becomes observable in present data sets like the cosmic microwave background (CMB) or the large scale structure (LSS). Within this proceeding we present a fully analytic method to infer inflationary parameters from observations by exploiting higher-order statistics of the curvature perturbations. To keep this analyticity, and thereby to dispense with numerically expensive sampling techniques, a saddle-point approximation is introduced whose precision has been validated for a numerical toy example. Applied to real data, this approach might enable to discriminate among the still viable models of inflation.

We present a generic inference method for inflation models from observational data by the usage of higher-order statistics of the curvature perturbation on uniform density hypersurfaces. This method is based on the calculation of the posterior for the primordial non-Gaussianity parameters $f_\text{NL}$ and $g_\text{NL}$, which in general depend on specific parameters of inflation and reheating models, and enables to discriminate among the still viable inflation models. To keep analyticity as far as possible to dispense with numerically expensive sampling techniques a saddle-point approximation is introduced, whose precision is validated for a numerical toy example. The mathematical formulation is done in a generic way so that the approach remains applicable to cosmic microwave background data as well as to large scale structure data. Additionally, we review a few currently interesting inflation models and present numerical toy examples thereof in two and three dimensions to demonstrate the efficiency of the higher-order statistics method. A second quantity of interest is the primordial power spectrum. Here, we present two Bayesian methods to infer it from observational data, the so called critical filter and an extension thereof with smoothness prior, both allowing for a non-parametric spectrum reconstruction. These methods are able to reconstruct the spectra of the observed perturbations and the primordial ones of curvature perturbation even in case of non-Gaussianity and partial sky coverage. We argue that observables like $T-$ and $B-$modes permit to measure both spectra. This also allows to infer the level of non-Gaussianity generated since inflation.

We present an approximate calculation of the full Bayesian posterior probability distribution for the local non-Gaussianity parameter $f_{\text{nl}}$ from observations of cosmic microwave background anisotropies within the framework of information field theory. The approximation that we introduce allows us to dispense with numerically expensive sampling techniques. We use a novel posterior validation method (DIP test) in cosmology to test the precision of our method. It transfers inaccuracies of the calculated posterior into deviations from a uniform distribution for a specially constructed test quantity. For this procedure we study toy cases that use one- and two-dimensional flat skies, as well as the full spherical sky. We find that we are able to calculate the posterior precisely under a flat-sky approximation, albeit not in the spherical case. We argue that this is most likely due to an insufficient precision of the used numerical implementation of the spherical harmonic transform, which might affect other non-Gaussianity estimators as well. Furthermore, we present how a nonlinear reconstruction of the primordial gravitational potential on the full spherical sky can be obtained in principle. Using the flat-sky approximation, we find deviations for the posterior of $f_{\text{nl}}$ from a Gaussian shape that become more significant for larger values of the underlying true $f_{\text{nl}}$. We also perform a comparison to the well-known estimator of Komatsu et al. [Astrophys. J. 634, 14 (2005)] and finally derive the posterior for the local non-Gaussianity parameter $g_{\text{nl}}$ as an example of how to extend the introduced formalism to higher orders of non-Gaussianity.

We present an error-diagnostic validation method for posterior distributions in Bayesian signal inference, an advancement of a previous work. It transfers deviations from the correct posterior into characteristic deviations from a uniform distribution of a quantity constructed for this purpose. We show that this method is able to reveal and discriminate several kinds of numerical and approximation errors, as well as their impact on the posterior distribution. For this we present four typical analytical examples of posteriors with incorrect variance, skewness, position of the maximum, or normalization. We show further how this test can be applied to multidimensional signals.