To make use of clustering statistics from large cosmological surveys, accurate and precise covariance matrices are needed. We present a new code to estimate large scale galaxy two-point correlation function (2PCF) covariances in arbitrary survey geometries that, due to new sampling techniques, runs $\sim 10^4$ times faster than previous codes, computing finely-binned covariance matrices with negligible noise in less than 100 CPU-hours. As in previous works, non-Gaussianity is approximated via a small rescaling of shot-noise in the theoretical model, calibrated by comparing jackknife survey covariances to an associated jackknife model. The flexible code, RascalC, has been publicly released, and automatically takes care of all necessary pre- and post-processing, requiring only a single input dataset (without a prior 2PCF model). Deviations between large scale model covariances from a mock survey and those from a large suite of mocks are found to be be indistinguishable from noise. In addition, the choice of input mock are shown to be irrelevant for desired noise levels below $\sim 10^5$ mocks. Coupled with its generalization to multi-tracer data-sets, this shows the algorithm to be an excellent tool for analysis, reducing the need for large numbers of mock simulations to be computed.

The relativistic generalization of the Newtonian Lagrangian perturbation theory is investigated. In previous works, the perturbation and solution schemes that are generated by the spatially projected gravitoelectric part of the Weyl tensor were given to any order of the perturbations, together with extensions and applications for accessing the nonperturbative regime. We here discuss more in detail the general first-order scheme within the Cartan formalism including and concentrating on the gravitational wave propagation in matter. We provide master equations for all parts of Lagrangian-linearized perturbations propagating in the perturbed spacetime, and we outline the solution procedure that allows one to find general solutions. Particular emphasis is given to global properties of the Lagrangian perturbation fields by employing results of Hodge-de Rham theory. We here discuss how the Hodge decomposition relates to the standard scalar-vector-tensor decomposition. Finally, we demonstrate that we obtain the known linear perturbation solutions of the standard relativistic perturbation scheme by performing two steps: first, by restricting our solutions to perturbations that propagate on a flat unperturbed background spacetime and, second, by transforming to Eulerian background coordinates with truncation of nonlinear terms.

We probe the higher-order galaxy clustering in the final data release (DR12) of the Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey (BOSS) using germ-grain Minkowski Functionals (MFs). Our data selection contains $979,430$ BOSS galaxies from both the northern and southern galactic caps over the redshift range $z=0.2 - 0.6$. We extract the higher-order part of the MFs, detecting the deviation from the purely Gaussian case with $\chi^2 \sim \mathcal{O}(10^3)$ on 24 degrees of freedom across the entire data selection. We measure significant redshift evolution in the higher-order functionals for the first time. We find $15-35\%$ growth, depending on functional and scale, between our redshift bins centered at $z=0.325$ and $z=0.525$. We show that the structure in higher order correlations grow faster than that in the two-point correlations, especially on small scales where the excess approaches a factor of $2$. We demonstrate how this trend is generalizable by finding good agreement of the data with a hierarchical model in which the higher orders grow faster than the lower order correlations. We find that the non-Gaussianity of the underlying dark matter field grows even faster than the one of the galaxies due to decreasing clustering bias. Our method can be adapted to study the redshift evolution of the three-point and higher functions individually.

We probe the higher-order clustering of the galaxies in the final data release (DR12) of the Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey (BOSS) using the method of germ-grain Minkowski Functionals (MFs). Our sample consists of 410,615 BOSS galaxies from the northern Galactic cap in the redshift range 0.450--0.595. We show the MFs to be sensitive to contributions up to the six-point correlation function for this data set. We ensure with a custom angular mask that the results are more independent of boundary effects than in previous analyses of this type. We extract the higher-order part of the MFs and quantify the difference to the case without higher-order correlations. The resulting $\chi^{2}$ value of over 10,000 for a modest number of degrees of freedom, O(200), indicates a 100-sigma deviation and demonstrates that we have a highly significant signal of the non-Gaussian contributions to the galaxy distribution. This statistical power can be useful in testing models with differing higher-order correlations. Comparing the galaxy data to the QPM and MultiDark-Patchy mocks, we find that the latter better describes the observed structure. From an order-by-order decomposition we expect that, for example, already a reduction of the amplitude of the MD-Patchy mock power spectrum by 5% would remove the remaining tension.

The relativistic generalization of the Newtonian Lagrangian perturbation theory is investigated. In previous works, the first-order trace solutions that are generated by the spatially projected gravitoelectric part of the Weyl tensor were given together with extensions and applications for accessing the nonperturbative regime. We furnish here construction rules to obtain from Newtonian solutions the gravitoelectric class of relativistic solutions, for which we give the complete perturbation and solution schemes at any order of the perturbations. By construction, these schemes generalize the complete hierarchy of solutions of the Newtonian Lagrangian perturbation theory.

We investigate the purely spatial Lagrangian coordinate transformation from the Lagrangian to the basic Eulerian frame. We demonstrate three techniques for extracting the relativistic displacement field from a given solution in the Lagrangian frame. These techniques are (a) from defining a local set of Eulerian coordinates embedded into the Lagrangian frame; (b) from performing a specific gauge transformation; and (c) from a fully non-perturbative approach based on the ADM split. The latter approach shows that this decomposition is not tied to a specific perturbative formulation for the solution of the Einstein equations. Rather, it can be defined at the level of the non-perturbative coordinate change from the Lagrangian to the Eulerian description. Studying such different techniques is useful because it allows us to compare and develop further the various approximation techniques available in the Lagrangian formulation. We find that one has to solve the gravitational wave equation in the relativistic analysis, otherwise the corresponding Newtonian limit will necessarily contain spurious non-propagating tensor artefacts at second order in the Eulerian frame. We also derive the magnetic part of the Weyl tensor in the Lagrangian frame, and find that it is not only excited by gravitational waves but also by tensor perturbations which are induced through the non-linear frame-dragging. We apply our findings to calculate for the first time the relativistic displacement field, up to second order, for a $\Lambda$CDM Universe in the presence of a local primordial non-Gaussian component. Finally, we also comment on recent claims about whether mass conservation in the Lagrangian frame is violated.

As deeper galaxy catalogues are soon to come, it becomes even more importantto measure large-scale fluctuations in the catalogues with robust statistics that cover all moments of the galaxy distribution.In this paper we reinforce a direct analysis of galaxy data by employing the Germ-Grain method to calculate thefamily of Minkowski Functionals. We introduce a new code, suitable for the analysis of large data sets without smoothingand without the construction of excursion sets. We provide new tools to measure correlation properties, putting emphasis onexplicitly isolating non-Gaussian correlations with the help of integral-geometric relations. As a first application we present the analysis of large-scale fluctuations in the luminous red galaxy sample of Sloan Digital Sky Survey data release 7 data. We findsignificant deviations from the $\Lambda$ cold dark matter mock catalogues on samples as large as $500h^{-1}$Mpc (more than $3\sigma$)and slight deviations of around $2\sigma$ on $700h^{-1}$Mpc, and we investigate possible sources of these deviations.

Kinematical and dynamical properties of a generic inhomogeneous cosmological model, spatially averaged with respect to free-falling (generalized fundamental) observers, are investigated for the matter model irrotational dust. Paraphrasing a previous Newtonian investigation, we present a relativistic generalization of a backreaction model based on volume-averaging the Relativistic Zeldovich Approximation. In this model we investigate the effect of kinematical backreaction on the evolution of cosmological parameters as they are defined in an averaged inhomogeneous cosmology, and we show that the backreaction model interpolates between orthogonal symmetry properties by covering subcases of the plane-symmetric solution, the Lemaitre-Tolman-Bondi solution and the Szekeres solution. We so obtain a powerful model that lays the foundations for quantitatively addressing curvature inhomogeneities as they would be interpreted as Dark Energy or Dark Matter in a quasi-Newtonian cosmology. The present model, having a limited architecture due to an assumed FLRW background, is nevertheless capable of replacing 1/4 of the needed amount for Dark Energy on domains of 200 Mpc in diameter for typical (one-sigma) fluctuations in a CDM initial power spectrum. However, the model is far from explaining Dark Energy on larger scales (spatially), where a 6% effect on 400 Mpc domains is identified that can be traced back to an on average negative intrinsic curvature today. One drawback of the quantitative results presented is the fact that the epoch when backreaction is effective on large scales and leads to volume acceleration lies in the future. We discuss this issue in relation to the initial spectrum, the Dark Matter problem, the coincidence problem, and the fact that large-scale Dark Energy is an effect on the past light cone (not spatial), and we pinpoint key elements of future research.

We test the isotropy of the expansion of the Universe by estimating the hemispherical anisotropy of supernova type Ia (SN Ia) Hubble diagrams at low redshifts (z<0.2). We compare the best fit Hubble diagrams in pairs of hemispheres and search for the maximal asymmetric orientation. For an isotropic Universe, we expect only a small asymmetry due to noise and the presence of nearby structures. This test does not depend on the assumed content of the Universe, the assumed model of gravity, or the spatial curvature of the Universe. The expectation for possible fluctuations due to large scale structure is evaluated for the \Lambda cold dark matter (\Lambda CDM) model and is compared to the supernova data from the Constitution set for four different light curve fitters, thus allowing a study of the systematic effects. The expected order of magnitude of the hemispherical asymmetry of the Hubble expansion agrees with the observed one. The direction of the Hubble asymmetry is established at 95% confidence level (C.L.) using both, the MLCS2k2 and the SALT II light curve fitter. The highest expansion rate is found towards (l, b) ~ (-35{\deg},-19{\deg}), which agrees with directions reported by other studies. Its amplitude is not in contradiction to expectations from the \Lambda CDM model. The measured Hubble anisotropy is \Delta H/H ~ 0.026. With 95% C.L. the expansion asymmetry is \Delta H/H<0.038.

Modern cosmology relies on the assumption of large-scale isotropy and homogeneity of the Universe. However, locally the Universe is inhomogeneous and anisotropic. So, how can local measurements (at the 100 Mpc scale) be used to determine global cosmological parameters (defined at the 10 Gpc scale)? We use Buchert's averaging formalism and determine a set of locally averaged cosmological parameters in the context of the flat Lambda cold dark matter model. We calculate their ensemble means (i.e. their global values) and variances (i.e. their cosmic variances). We apply our results to typical survey geometries and focus on the study of the effects of local fluctuations of the curvature parameter. By this means we show, that in the linear regime cosmological backreaction and averaging can be reformulated as the issue of cosmic variance. The cosmic variance is found largest for the curvature parameter and discuss some of its consequences. We further propose to use the observed variance of cosmological parameters to measure the growth factor. [abbreviated]

The backreaction of inhomogeneities on the global expansion history of the Universe suggests a possible link of the formation of structures to the recent accelerated expansion. In this paper, the origin of this conjecture is illustrated and a model without Dark Energy that allows for a more explicit investigation of this link is discussed. Additionally to this conceptually interesting feature, the model leads to a LCDM-like distance-redshift relation that is consistent with SN data.

Cosmological backreaction suggests a link between structure formation and the expansion history of the Universe. In order to quantitatively examine this connection, we dynamically investigate a volume partition of the Universe into over-- and underdense regions. This allows us to trace structure formation using the volume fraction of the overdense regions $\lambda_{\CM}$ as its characterizing parameter. Employing results from cosmological perturbation theory and extrapolating the leading mode into the nonlinear regime, we construct a three--parameter model for the effective cosmic expansion history, involving $\lambda_{\CM_{0}}$, the matter density $\Omega_{m}^{\CD_{0}}$, and the Hubble rate $H_{\CD_{0}}$ of today's Universe. Taking standard values for $\Omega_{m}^{\CD_{0}}$ and $H_{\CD_{0}}$ as well as a reasonable value for $\lambda_{\CM_{0}}$, that we derive from $N$--body simulations, we determine the corresponding amounts of backreaction and spatial curvature. We find that the obtained values that are sufficient to generate today's structure also lead to a $\Lambda$CDM--like behavior of the scale factor, parametrized by the same parameters $\Omega_{m}^{\CD_{0}}$ and $H_{\CD_{0}}$, but without a cosmological constant. However, the temporal behavior of $\lambda_{\CM}$ does not faithfully reproduce the structure formation history. Surprisingly, however, the model matches with structure formation with the assumption of a low matter content, $\Omega_{m}^{\CD_{0}}\approx3\%$, a result that hints to a different interpretation of part of the backreaction effect as kinematical Dark Matter. (truncated)