We investigate effective equations governing the volume expansion of spatially averaged portions of inhomogeneous cosmologies in spacetimes filled with an arbitrary fluid. This work is a follow-up to previous studies focused on irrotational dust models (Paper I) and irrotational perfect fluids (Paper II) in flow-orthogonal foliations of spacetime. It complements them by considering arbitrary foliations, arbitrary lapse and shift, and by allowing for a tilted fluid flow with vorticity. As for the first studies, the propagation of the spatial averaging domain is chosen to follow the congruence of the fluid, which avoids unphysical dependencies in the averaged system that is obtained. We present two different averaging schemes and corresponding systems of averaged evolution equations providing generalizations of Papers I and II. The first one retains the averaging operator used in several other generalizations found in the literature. We extensively discuss relations to these formalisms and pinpoint limitations, in particular regarding rest mass conservation on the averaging domain. The alternative averaging scheme that we subsequently introduce follows the spirit of Papers I and II and focuses on the fluid flow and the associated 1+3 threading congruence, used jointly with the 3+1 foliation that builds the surfaces of averaging. This results in compact averaged equations with a minimal number of cosmological backreaction terms. We highlight that this system becomes especially transparent when applied to a natural class of foliations which have constant fluid proper time slices.

Recently it has been noted by Di Valentino, Melchiorri and Silk (2019) that the enhanced lensing signal relative to that expected in the spatially flat $\Lambda$CDM model poses a possible crisis for the Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) class of models usually used to interpret cosmological data. The 'crisis' amounts to inconsistencies between cosmological datasets arising when the FLRW curvature parameter $\Omega_k$ is determined from the data rather than constrained to be zero a priori. Moreover, the already substantial discrepancy between the Hubble parameter as determined by Planck and local observations increases to the level of $5\sigma$. While such inconsistencies might arise from systematic effects of astrophysical origin affecting the Planck Cosmic Microwave Background (CMB) power spectra at small angular scales, it is an option that the inconsistencies are due to the failure of the FLRW assumption. In this paper we recall how the FLRW curvature ansatz is expected to be violated for generic relativistic spacetimes. We explain how the FLRW conservation equation for volume-averaged spatial curvature is modified through structure formation, and we illustrate in a simple framework how the curvature tension in a FLRW spacetime can be resolved - and is even expected to occur - from the point of view of general relativity. Requiring early-time convergence towards a Friedmannian model with a spatial curvature parameter $\Omega_{k0}$ equal to that preferred from the Planck power spectra resolves the Hubble tension within our dark energy-free model.

Few statements in cosmology can be made without assuming a cosmological model within which to interpret data. Statements about cosmic acceleration are no exception to this rule, and the inferred positive volume acceleration of our Universe often quoted in the literature is valid in the context of the standard Friedmann-Lemaitre-Robertson-Walker (FLRW) class of space-times. Using the Joint Light-curve Analysis (JLA) catalogue of type Ia supernovae (SNIa), we examine the fit of a class of exact scaling solutions with dynamical spatial curvature formulated in the framework of a scalar averaging scheme for relativistic inhomogeneous space-times. In these models, global volume acceleration may emerge as a result of the non-local variance between expansion rates of clusters and voids, the latter gaining volume dominance in the late-epoch Universe. We find best-fit parameters for a scaling model of backreaction that are reasonably consistent with previously found constraints from SNIa, CMB, and baryon acoustic oscillations data. The quality of fit of the scaling solutions is indistinguishable from that of the LCDM model and the timescape cosmology from an Akaike Information Criterion (AIC) perspective. This indicates that a broad class of models can account for the $z\lesssim 1$ expansion history.

The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We introduce a generalization of the 4-dimensional averaging window function of Gasperini, Marozzi and Veneziano (2010) that may prove useful for a number of applications. The covariant nature of spatial scalar averaging schemes to address the averaging problem in relativistic cosmology is an important property that is implied by construction, but usually remains implicit. We employ here the approach of Gasperini et al. for two reasons. First, the formalism and its generalization presented here are manifestly covariant. Second, the formalism is convenient for disentangling the dependencies on foliation, volume measure, and boundaries in the averaged expressions entering in scalar averaging schemes. These properties will prove handy for simplifying expressions, but also for investigating extremal foliations and for comparing averaged properties of different foliations directly. The proposed generalization of the window function allows for choosing the most appropriate averaging scheme for the physical problem at hand, and for distinguishing between the role of the foliation itself and the role of the volume measure in averaged dynamic equations. We also show that one particular window function obtained from this generalized class results in an averaging scheme corresponding to that of a recent investigation by Buchert, Mourier and Roy (2018) and, as a byproduct, we explicitly show that the general equations for backreaction derived therein are covariant.

This study presents a numerical analysis of the topology of a set of cosmologically interesting three-dimensional Gaussian random fields in terms of their Betti numbers $\beta_0$, $\beta_1$ and $\beta_2$. We show that Betti numbers entail a considerably richer characterization of the topology of the primordial density field. Of particular interest is that Betti numbers specify which topological features - islands, cavities or tunnels - define its spatial structure. A principal characteristic of Gaussian fields is that the three Betti numbers dominate the topology at different density ranges. At extreme density levels, the topology is dominated by a single class of features. At low levels this is a \emph{Swiss-cheeselike} topology, dominated by isolated cavities, at high levels a predominantly \emph{Meatball-like} topology of isolated objects. At moderate density levels, two Betti number define a more \emph{Sponge-like} topology. At mean density, the topology even needs three Betti numbers, quantifying a field consisting of several disconnected complexes, not of one connected and percolating overdensity. A {\it second} important aspect of Betti number statistics is that they are sensitive to the power spectrum. It reveals a monotonic trend in which at a moderate density range a lower spectral index corresponds to a considerably higher (relative) population of cavities and islands. We also assess the level of complementary information that Betti numbers represent, in addition to conventional measures such as Minkowski functionals. To this end, we include an extensive description of the Gaussian Kinematic Formula (GKF), which represents a major theoretical underpinning for this discussion.

We study the topology generated by the temperature fluctuations of the Cosmic Microwave Background (CMB) radiation, as quantified by the number of components and holes, formally given by the Betti numbers, in the growing excursion sets. We compare CMB maps observed by the Planck satellite with a thousand simulated maps generated according to the LCDM paradigm with Gaussian distributed fluctuations. The survey of the CMB over $\mathbb{S}^2$ is incomplete due to obfuscation effects by bright point sources and other extended foreground objects like our own galaxy. To deal with such situations, where analysis in the presence of "masks" is of importance, we introduce the concept of relative homology. The parametric $\chi^2$-test shows differences between observations and simulations, yielding $p$-values at per-cent to less than per-mil levels roughly between 2 to 7 degrees. The highest observed deviation for $b_0$ and $b_1$ is approximately between $3\sigma$-4$\sigma$ at scales of 3 to 7 degrees. There are reports of mildly unusual behaviour of the Euler characteristic at 3.66 degrees in the literature, computed from independent measurements of the CMB temperature fluctuations by Planck's predecessor WMAP satellite. The mildly anomalous behaviour of Euler characteristic is related to the strongly anomalous behaviour of components and holes. These are also the scales at which the observed maps exhibit low variance compared to the simulations. Non-parametric tests show even stronger differences at almost all scales. Regardless, beyond the trivial possibility that this may still be a manifestation of an extreme Gaussian case, these observations, along with the super-horizon scales involved, may motivate to look at primordial non-Gaussianity. Alternative scenarios worth exploring may be models with non-trivial topology.

The subject of cosmological backreaction in General Relativity is often approached by coordinate-dependent and metric-based analyses. We present in this letter an averaging formalism for the scalar parts of Einstein's equations that is coordinate-independent and only functionally depends on a metric. This formalism is applicable to general 3+1 foliations of spacetime for an arbitrary fluid with tilted flow. We clarify the dependence on spacetime foliation and argue that this dependence is weak in cosmological settings. We also introduce a new set of averaged equations that feature a global cosmological time despite the generality of the setting, and we put the statistical nature of effective cosmologies into perspective.

The standard model of cosmology assumes that the Universe can be described to hover around a homogeneous-isotropic solution of Einstein's general theory of relativity. This description needs (sometimes hidden) hypotheses that restrict the generality, and relaxing these restrictions is the headline of a new physical approach to cosmology that refurnishes the cosmological framework. Considering a homogeneous geometry as a template geometry for the in reality highly inhomogeneous Universe must be considered a strong idealization. Unveiling the limitations of the standard model opens the door to rich consequences of general relativity, giving rise to effective (i.e. spatially averaged) cosmological models that may even explain the longstanding problems of dark energy and dark matter. We explore in this talk the influence of structure formation on average properties of the Universe by discussing: (i) general thoughts on why considering average properties, on the key-issue of non-conserved curvature, and on the global gravitational instability of the standard model of cosmology; (ii) the general set of cosmological equations arising from averaging the scalar parts of Einstein's equations, the generic property of structure formation interacting with the average properties of the Universe in a scale-dependent way, and the description of cosmological backreaction in terms of an effective scalar field.

We extend the general relativistic Lagrangian perturbation theory, recently developed for the formation of cosmic structures in a dust continuum, to the case of model universes containing a single fluid with a single-valued analytic equation of state. Using a coframe-based perturbation approach, we investigate evolution equations for structure formation in pressure-supported irrotational fluids that generate their rest-frame spacetime foliation. We provide master equations to first order for the evolution of the trace and traceless parts of barotropic perturbations that evolve in the perturbed space, where the latter describes the propagation of gravitational waves in the fluid. We illustrate the trace evolution for a linear equation of state and for a model equation of state describing isotropic velocity dispersion, and we discuss differences to the dust matter model, to the Newtonian case, and to standard perturbation approaches.

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 clarify that a result recently stated by Kaiser is contained in a theorem of Buchert and Ehlers that is widely known for its main result: that there is no global kinematical backreaction in Newtonian cosmology. Kaiser cites this paper, re-derives parts of the theorem, but incompletely restates its content. He makes further claims, which cannot be proven beyond the limited context of Newtonian cosmology. We also discuss recent papers of R\'acz et al. and Roukema who claim the existence of global backreaction within the Newtonian framework.

Despite the wealth of $Planck$ results, there are difficulties in disentangling the primordial non-Gaussianity of the Cosmic Microwave Background (CMB) from the secondary and the foreground non-Gaussianity (NG). For each of these forms of NG the lack of complete data introduces model-dependencies. Aiming at detecting the NGs of the CMB temperature anisotropy $\delta T$, while paying particular attention to a model-independent quantification of NGs, our analysis is based upon statistical and morphological univariate descriptors, respectively: the probability density function $P(\delta T)$, related to ${\mathrm v}_{0}$, the first Minkowski Functional (MF), and the two other MFs, ${\mathrm v}_{1}$ and ${\mathrm v}_{2}$. From their analytical Gaussian predictions we build the discrepancy functions $\Delta_{k}$ ($k=P,0,1,2$) which are applied to an ensemble of $10^{5}$ CMB realization maps of the $\Lambda$CDM model and to the $Planck$ CMB maps. In our analysis we use general Hermite expansions of the $\Delta_{k}$ up to the $12^{th}$ order, where the coefficients are explicitly given in terms of cumulants. Assuming hierarchical ordering of the cumulants, we obtain the perturbative expansions generalizing the $2^{nd}$ order expansions of Matsubara to arbitrary order in the standard deviation $\sigma_0$ for $P(\delta T)$ and ${\mathrm v}_0$, where the perturbative expansion coefficients are explicitly given in terms of complete Bell polynomials. The comparison of the Hermite expansions and the perturbative expansions is performed for the $\Lambda$CDM map sample and the $Planck$ data. We confirm the weak level of non-Gaussianity ($1$-$2$)$\sigma$ of the foreground corrected masked $Planck$ $2015$ maps.

Self-consistent treatment of cosmological structure formation and expansion within the context of classical general relativity may lead to "extra" expansion above that expected in a structureless universe. We argue that in comparison to an early-epoch, extrapolated Einstein-de Sitter model, about 10-15% "extra" expansion is sufficient at the present to render superfluous the "dark energy" 68% contribution to the energy density budget, and that this is observationally realistic.

The current cosmological model ($\Lambda$CDM) with the underlying FLRW metric relies on the assumption of local isotropy, hence homogeneity of the Universe. Difficulties arise when one attempts to justify this model as an average description of the Universe from first principles of general relativity, since in general, the Einstein tensor built from the averaged metric is not equal to the averaged stress--energy tensor. In this context, the discrepancy between these quantities is called "cosmological backreaction" and has been the subject of scientific debate among cosmologists and relativists for more than $20$ years. Here we present one of the methods to tackle this problem, i.e. averaging the scalar parts of the Einstein equations, together with its application, the cosmological mass function of galaxy clusters.

In relativistic inhomogeneous cosmology, structure formation couples to average cosmological expansion. A conservative approach to modelling this assumes an Einstein--de Sitter model (EdS) at early times and extrapolates this forward in cosmological time as a "background model" against which average properties of today's Universe can be measured. This requires adopting an early-epoch--normalised background Hubble constant $H_1^{bg}$. Here, we show that the $\Lambda$CDM model can be used as an observational proxy to estimate $H_1^{bg}$ rather than choose it arbitrarily. We assume (i) an EdS model at early times; (ii) a zero dark energy parameter; (iii) bi-domain scalar averaging---division of the spatial sections into over- and underdense regions; and (iv) virialisation (stable clustering) of collapsed regions. We find $H_1^{bg}= 37.7 \pm 0.4$ km/s/Mpc (random error only) based on a Planck $\Lambda$CDM observational proxy. Moreover, since the scalar-averaged expansion rate is expected to exceed the (extrapolated) background expansion rate, the expected age of the Universe should be much less than $2/(3 H_1^{bg}) = 17.3$ Gyr. The maximum stellar age of Galactic Bulge microlensed low-mass stars (most likely: 14.7 Gyr; 68\% confidence: 14.0--15.0 Gyr) suggests an age about a Gyr older than the (no-backreaction) $\Lambda$CDM estimate.

We summarise some of the main observational challenges for the standard Friedmann-Lemaitre-Robertson-Walker cosmological model and describe how results recently presented in the parallel session `Large--scale Structure and Statistics' (DE3) at the `Fourteenth Marcel Grossman Meeting on General Relativity' are related to these challenges.

The general relativistic description of cosmological structure formation is an important challenge from both the theoretical and the numerical point of views. In this paper we present a brief prescription for a general relativistic treatment of structure formation and a resulting mass function on galaxy cluster scales in a highly generic scenario. To obtain this we use an exact scalar averaging scheme together with the relativistic generalization of Zel'dovich's approximation (RZA) that serves as a closure condition for the averaged equations.

In the framework of Lagrangian perturbation theory in general relativity we discuss the possibility to split the Einstein equations, written in terms of spatial Cartan coframes within a 3+1 foliation of spacetime, into gravitoelectric and gravitomagnetic parts. While the former reproduces the full hierarchy of the Newtonian perturbation solutions, the latter contains non-Newtonian aspects like gravitational waves. This split can be understood and made unique through the Hodge decomposition of Cartan coframe fields.

No. In a number of papers Green and Wald argue that the standard FLRW model approximates our Universe extremely well on all scales, except close to strong field astrophysical objects. In particular, they argue that the effect of inhomogeneities on average properties of the Universe (backreaction) is irrelevant. We show that this latter claim is not valid. Specifically, we demonstrate, referring to their recent review paper, that (i) their two-dimensional example used to illustrate the fitting problem differs from the actual problem in important respects, and it assumes what is to be proven; (ii) the proof of the trace-free property of backreaction is unphysical and the theorem about it fails to be a mathematically general statement; (iii) the scheme that underlies the trace-free theorem does not involve averaging and therefore does not capture crucial non-local effects; (iv) their arguments are to a large extent coordinate-dependent, and (v) many of their criticisms of backreaction frameworks do not apply to the published definitions of these frameworks. It is therefore incorrect to infer that Green and Wald have proven a general result that addresses the essential physical questions of backreaction in cosmology.

In the standard model of cosmology, the Universe is static in comoving coordinates; expansion occurs homogeneously and is represented by a global scale factor. The baryon acoustic oscillation (BAO) peak location is a statistical tracer that represents, in the standard model, a fixed comoving-length standard ruler. Recent gravitational collapse should modify the metric, rendering the effective scale factor, and thus the BAO standard ruler, spatially inhomogeneous. Using the Sloan Digital Sky Survey, we show to high significance (P < 0.001) that the spatial compression of the BAO peak location increases as the spatial paths' overlap with superclusters increases. Detailed observational and theoretical calibration of this BAO peak location environment dependence will be needed when interpreting the next decade's cosmological surveys.

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.

The Friedmann-Lemaitre-Robertson-Walker (FLRW) metric assumes comoving spatial rigidity of metrical properties. The curvature term in comoving coordinates is environment-independent and cannot evolve. In the standard model, structure formation is interpreted accordingly: structures average out on the chosen metrical background, which remains rigid in comoving coordinates despite nonlinear structure growth. The latter claim needs to be tested, since it is a hypothesis that is not derived using general relativity. We introduce a test of the comoving rigidity assumption by measuring the two-point auto-correlation function on comoving scales---assuming FLRW comoving spatial rigidity---in order to detect shifts in the baryon acoustic oscillation (BAO) peak location for Large Red Galaxy (LRG) pairs of the Sloan Digital Sky Survey Data Release 7. In tangential directions, subsets of pairs overlapping with superclusters or voids show the BAO peak. The tangential BAO peak location for overlap with Nadathur & Hotchkiss superclusters is $4.3\pm1.6$ Mpc/h less than that for LRG pairs unselected for supercluster overlap, and $6.6\pm2.8$ Mpc/h less than that of the complementary pairs. Liivamagi et al. superclusters give corresponding differences of $3.7\pm2.9$ Mpc/h and $6.3\pm2.6$ Mpc/h, respectively. We have found moderately significant evidence (Kolmogorov--Smirnov tests suggest very significant evidence) that the BAO peak location for supercluster-overlapping pairs is compressed by about 6% compared to that of the complementary sample, providing a potential challenge to FLRW models and a benchmark for predictions from models based on an averaging approach that leaves the spatial metric a priori unspecified.

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.

The concordance model is successful in explaining numerous observable phenomena at the price of introducing an exotic source of unknown origin: dark energy. Dark energy dominance occurs at recent epochs, when we expect most cosmological structures to have already formed, and thus, when the error induced by forcing the homogeneous FLRW metric onto the data is expected to be the most significant. We propose a way to quantify the impact of deviations from homogeneity on the evolution of cosmological parameters. Using a multi-scale partitioning approach and the virialization fraction estimated from numerical simulations in an Einstein-de Sitter model, we obtain an observationally realistic distance modulus over redshifts 0 < z < 3 by a relativistic correction of the FLRW metric.

The 3-torus (T^3) FLRW model better fits the nearly zero large-scale auto-correlation of the WMAP CMB sky maps than the infinite flat model. The T^3 model's WMAP parameters imply approximately equal-redshift topological lensing at z \sim 6. We investigate observational strategies for rejecting the T^3 solution or providing candidate topologically lensed galaxy pairs. T^3 holonomies are applied to (i) existing z \sim 6 observations and (ii) simulated observations, creating multiply connected catalogues. Corresponding simply connected catalogues are generated. Each catalogue is analysed using a successive filter method and collecting matched quadruples. Quadruple statistics between the multiply and simply connected catalogues are compared. The expected rejection of the hypothesis, or detection of candidate topologically lensed galaxies, is possible at a significance of 5% for a pair of T^3 axis-centred northern and southern surveys if photometric redshift accuracy is \sigma(\zphot) < 0.01 for a pair of nearly complete 100 deg^2 surveys with a total of > 500 galaxies over 4.3 < z < 6.6, or for a pair of 196 deg^2 surveys with > 400 galaxies and \sigma(\zphot) < 0.02 over 4<z<7. Dropping the maximum time interval in a pair from \Delta t =1 Gyr/h to \Delta t =0.1 Gyr/h yields a requirement of \sigma(\zphot) < 0.005 or \sigma(\zphot) < 0.01, respectively. Millions of z \sim 6 galaxies will be observed over fields of these sizes during the coming decades, implying much stronger constraints. The question is not if the hypothesis will be rejected or confirmed, it is when.

The geometry of the dark energy and cold dark matter dominated cosmological model (LambdaCDM) is commonly assumed to be given by a Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, i.e. it assumes homogeneity in the comoving spatial section. The homogeneity assumption fails most strongly at (i) small distance scales and (ii) recent epochs, implying that the FLRW approximation is most likely to fail at these scales. We use the virialisation fraction to quantify (i) and (ii), which approximately coincide with each other on the observational past light cone. For increasing time, the virialisation fraction increases above 10% at about the same redshift (sim 1-3) at which Omega_Lambda grows above 10% (approx 1.8). Thus, instead of non-zero Omega_Lambda, we propose an approximate, general-relativistic correction to the matter-dominated (Omega_m =1, Omega_Lambda=0), homogeneous metric (Einstein de Sitter, EdS). A low-redshift effective matter-density parameter of Omega_m^eff(0) = 0.26 pm 0.05 is inferred. Over redshifts 0 < z < 3, the distance modulus of the virialisation-corrected EdS model approximately matches the LambdaCDM distance modulus. This rough approximation assumes "old physics" (general relativity), not "new physics". Thus, pending more detailed calculations, we strengthen the claim that "dark energy" should be considered as an artefact of emerging average curvature in the void-dominated Universe, via a novel approach that quantifies the relation between virialisation and average curvature evolution.

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 review the effect of the formation of nonlinear structures on the expansion rate, spatial curvature and light propagation in the universe, focusing on the possibility that it could explain cosmological observations without the introduction of dark energy or modified gravity. We concentrate on explaining the relevant physics and highlighting open questions.

We investigate the Lagrangian perturbation theory of a homogeneous and isotropic universe in the non-relativistic limit, and derive the solutions up to the fourth order. These solutions are needed for example for the next-to-leading order correction of the (resummed) Lagrangian matter bispectrum, which we study in an accompanying paper. We focus on flat cosmologies with a vanishing cosmological constant, and provide an in-depth description of two complementary approaches used in the current literature. Both approaches are solved with two different sets of initial conditions---both appropriate for modelling the large-scale structure. Afterwards we consider only the fastest growing mode solution, which is not affected by either of these choices of initial conditions. Under the reasonable approximation that the linear density contrast is evaluated at the initial Lagrangian position of the fluid particle, we obtain the nth-order displacement field in the so-called initial position limit: the nth order displacement field consists of 3(n-1) integrals over n linear density contrasts, and obeys self-similarity. Then, we find exact relations between the series in Lagrangian and Eulerian perturbation theory, leading to identical predictions for the density contrast and the peculiar-velocity divergence up to the fourth order.

The standard model of cosmology is based on homogeneous-isotropic solutions of Einstein's equations. These solutions are known to be gravitationally unstable to local inhomogeneous perturbations, commonly described as evolving on a background given by the same solutions. In this picture, the FLRW backgrounds are taken to describe the average over inhomogeneous perturbations for all times. We study in the present article the (in)stability of FLRW dust backgrounds within a class of averaged inhomogeneous cosmologies. We examine the phase portraits of the latter, discuss their fixed points and orbital structure and provide detailed illustrations. We show that FLRW cosmologies are unstable in some relevant cases: averaged models are driven away from them through structure formation and accelerated expansion. We find support for the proposal that the dark components of the FLRW framework may be associated to these instability sectors. Our conclusion is that FLRW cosmologies have to be considered critically as for their role to serve as reliable models for the physical background.

We outline the key-steps toward the construction of a physical, fully relativistic cosmology. The influence of inhomogeneities on the effective evolution history of the Universe is encoded in backreaction terms and expressed through spatially averaged geometrical invariants. These are absent and potential candidates for the missing dark sources in the standard model. Since they can be interpreted as energies of an emerging scalar field (the morphon), we are in the position to propose a strategy of how phenomenological scalar field models for Dark Energy, Dark Matter and Inflation, that are usually added as fundamental sources to a homogeneous-geometry (FLRW) cosmology, can be potentially traced back to inhomogeneous geometrical properties of space and its embedding into spacetime. We lay down a line of arguments that is - thus far only qualitatively - conclusive, and we address open problems of quantitative nature, related to the interpretation of observations. We discuss within a covariant framework ... [abridged]

We consider spatially averaged inhomogeneous universe models and argue that, already in the absence of sources, an effective scalar field arises through foliating and spatially averaging inhomogeneous geometrical curvature invariants of the Einstein vacuum. This scalar field (the `morphon') acts as an inflaton, if we prescribe a potential of some generic form. We show that, for any initially negative average spatial curvature, the morphon is driven through an inflationary phase and leads - on average - to a spatially flat, homogeneous and isotropic universe model, providing initial conditions for pre-heating and, by the same mechanism, a possibly natural self-exit.

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.

We outline the key-steps towards the construction of a physical, fully relativistic cosmology, in which we aim to trace Dark Energy and Dark Matter back to physical properties of space. The influence of inhomogeneities on the effective evolution history of the Universe is encoded in backreaction terms and expressed through spatially averaged geometrical invariants. These are absent and interpreted as missing dark fundamental sources in the standard model. In the inhomogeneous case they can be interpreted as energies of an emerging scalar field (the morphon). These averaged invariants vanish for a homogeneous geometry, where the morphon is in an unstable equilibrium state. If this state is perturbed, the morphon can act as a classical inflaton in the Early Universe and its de-balanced energies can mimic the dark sources in the Late Universe, depending on spatial scale as Dark Energy or as Dark Matter, respectively. We lay down a line of arguments that is qualitatively conclusive, and we outline open problems of quantitative nature, related to the interpretation of observations.

In this contribution we summarize two recent applications of a correspondence between backreaction terms in averaged inhomogeneous cosmologies and an effective scalar field (the `morphon'). Backreaction terms that add to the standard sources of Friedmannian kinematical laws and that emerge from geometrical curvature invariants built from inhomogeneities, can be interpreted in terms of a minimally coupled scalar field in the case of a dust matter source. We consider closure conditions of the averaged equations that lead to different evolution scenarii: a) the standard Chaplygin equation of state imposed as an effective relation between kinematical fluctuations and intrinsic curvature of space sections, and b) an inflationary scenario that emerges out of inhomogeneities of the Einstein vacuum, where averaged curvature inhomogeneities model the potential of an effective classical inflaton.

In the context of averaging an inhomogeneous cosmological model, we propose a natural measure identical to the Kullback-Leibler relative information entropy, which expresses the distinguishability of the local inhomogeneous density field from its spatial average on arbitrary compact domains. This measure is expected to be an increasing function in time and thus to play a significant role in studying gravitational entropy. To verify this conjecture, we explore the time evolution of the measure using the linear perturbation theory of a spatially flat FLRW model and a spherically symmetric nonlinear solution. We discuss the generality and conditions for the time-increasing nature of the measure, and also the connection to the backreaction effect caused by inhomogeneities.

In the framework of spatially averaged inhomogeneous cosmologies in classical general relativity, effective Einstein equations govern the dynamics of averaged scalar variables in a scale--dependent way. A particular cosmology may be characterized by a cosmic equation of state, closing the hierarchy of effective equations. In this context a natural candidate is provided by the Chaplygin gas, standing for a unified description of dark energy and dark matter. In this paper, we suppose that the inhomogeneous properties of matter and geometry obey the Chaplygin equation of state. The most extreme interpretation assumes that both dark energy and dark matter are not included as additional sources, but are both manifestations of spatial geometrical properties. This feature is an important conceptual difference in comparison with the standard approach of a Friedmann-Lema\^itre-Robertson-Walker universe filled with dust and another fundamental field characterized by the Chaplygin equation of state. We finally discuss the consequences of the resulting scenario for effective cosmological parameters in order to establish the framework of a future confrontation with observations, and we note that the standard Chaplygin gas may not be ruled out by them.

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)

In this note we wish to complement some recent work in the cosmological literature concerning the Weyl conformal curvature tensor and its parts. In particular, we shall give a clear-cut definition of the Newtonian limits of electric and magnetic parts of the Weyl tensor. We also discuss that only a subset of the relativistic equations is needed to obtain a closed system of equations in the Newtonian limit.

The thoughts expressed in this article are based on remarks made by J\"urgen Ehlers at the Albert-Einstein-Institut, Golm, Germany in July 2007. The main objective of this article is to demonstrate, in terms of plausible order-of-magnitude estimates for geometrical scalars, the relevance of spatial curvature in realistic models of the Universe that describe the dynamics of structure formation since the epoch of matter-radiation decoupling. We introduce these estimates with a commentary on the use of a quasi-Newtonian metric form in this context.

In order to quantitatively test the ability of averaged inhomogeneous cosmologies to correctly describe observations of the large scale properties of the Universe, we introduce a smoothed template metric corresponding to a constant spatial curvature model at any time, but with an evolving curvature parameter. This metric is used to compute quantities along an approximate effective lightcone of the averaged model of the Universe. As opposed to the standard Friedmann model, we parameterize this template metric by exact scaling properties of an averaged inhomogeneous cosmology, and we also motivate this form of the metric by results on a geometrical smoothing of inhomogeneous cosmological hypersurfaces. We test our hypothesis for the template metric against supernova data and the position of the CMB peaks, and infer the goodness--of--fit and parameter uncertainties. We find that averaged inhomogeneous models can reproduce the observations without requiring an additional Dark Energy component (though a volume acceleration is still needed), and that current data do not disfavour our main assumption on the effective lightcone structure. We also show that the experimental uncertainties on the angular diameter distance and the Hubble parameter from Baryon Acoustic Oscillations measurements -- forseen in future surveys like the proposed EUCLID satellite project -- are sufficiently small to distinguish between a FLRW template geometry and the template geometry with consistently evolving curvature. (Abridged)

We discuss the effect of curvature and matter inhomogeneities on the averaged scalar curvature of the present-day Universe. Motivated by studies of averaged inhomogeneous cosmologies, we contemplate on the question whether it is sensible to assume that curvature averages out on some scale of homogeneity, as implied by the standard concordance model of cosmology, or whether the averaged scalar curvature can be largely negative today, as required for an explanation of Dark Energy from inhomogeneities. We confront both conjectures with a detailed analysis of the kinematical backreaction term and estimate its strength for a multi-scale inhomogeneous matter and curvature distribution. Our main result is a formula for the spatially averaged scalar curvature involving quantities that are all measurable on regional (i.e. up to 100 Mpc) scales. We propose strategies to quantitatively evaluate the formula, and pinpoint the assumptions implied by the conjecture of a small or zero averaged curvature. We reach the conclusion that the standard concordance model needs fine-tuning in the sense of an assumed equipartition law for curvature in order to reconcile it with the estimated properties of the averaged physical space, whereas a negative averaged curvature is favoured, independent of the prior on the value of the cosmological constant.

The effective evolution of an inhomogeneous universe model in any theory of gravitation may be described in terms of spatially averaged variables. In Einstein's theory, restricting attention to scalar variables, this evolution can be modeled by solutions of a set of Friedmann equations for an effective volume scale factor, with matter and backreaction source terms. The latter can be represented by an effective scalar field (`morphon field') modeling Dark Energy. The present work provides an overview over the Dark Energy debate in connection with the impact of inhomogeneities, and formulates strategies for a comprehensive quantitative evaluation of backreaction effects both in theoretical and observational cosmology. We recall the basic steps of a description of backreaction effects in relativistic cosmology that lead to refurnishing the standard cosmological equations, but also lay down a number of challenges and unresolved issues in connection with their observational interpretation. The present status of this subject is intermediate: we have a good qualitative understanding of backreaction effects pointing to a global instability of the standard model of cosmology; exact solutions and perturbative results modeling this instability lie in the right sector to explain Dark Energy from inhomogeneities. It is fair to say that, even if backreaction effects turn out to be less important than anticipated by some researchers, the concordance high-precision cosmology, the architecture of current N-body simulations, as well as standard perturbative approaches may all fall short in correctly describing the Late Universe.

The effective evolution of an inhomogeneous universe model in Einstein's theory of gravitation may be described in terms of spatially averaged scalar variables. This evolution can be modeled by solutions of a set of Friedmann equations for an effective scale factor, with matter and backreaction source terms, where the latter can be represented by a minimally coupled scalar field (`morphon field'). We review the basic steps of a description of backreaction effects in relativistic cosmology that lead to refurnishing the standard cosmological equations, but also lay down a number of unresolved issues in connection with their interpretation within observational cosmology.

In the context of averaged cosmologies, the effective equations can be written in the form of "regional" Friedmannian equations with additional sources arising from the so-called backreaction of inhomogeneities. We propose a mean field description of this backreaction in terms of a regionally homogeneous scalar field: this provides a physical motivation to the phenomenological scalar fields generically called quintessence fields. We explicitly reconstruct the potential of the scalar field for a one-parameter family of scaling solutions to the backreaction problem, showing that it entails most of the standard scalar fields including e.g. standard and phantom quintessence scenarii.

Spatially averaged inhomogeneous cosmologies in classical general relativity can be written in the form of effective Friedmann equations with sources that include backreaction terms. In this paper we propose to describe these backreaction terms with the help of a homogeneous scalar field evolving in a potential; we call it the `morphon field'. This new field links classical inhomogeneous cosmologies to scalar field cosmologies, allowing to reinterpret, e.g., quintessence scenarios by routing the physical origin of the scalar field source to inhomogeneities in the Universe. We investigate a one-parameter family of scaling solutions to the backreaction problem. Subcases of these solutions (all without an assumed cosmological constant) include scale-dependent models with Friedmannian kinematics that can mimic the presence of a cosmological constant or a time-dependent cosmological term. We explicitly reconstruct the scalar field potential for the scaling solutions, and discuss those cases that provide a solution to the Dark Energy and coincidence problems. In this approach, Dark Energy emerges from morphon fields, a mechanism that can be understood through the proposed correspondence: the averaged cosmology is characterized by a weak decay (quintessence) or growth (phantom quintessence) of kinematical fluctuations, fed by `curvature energy' that is stored in the averaged 3-Ricci curvature. We find that the late-time trajectories of those models approach attractors that lie in the future of a state that is predicted by observational constraints.

Regionally averaged relativistic cosmologies have recently been considered as a possible explanation for the apparent late time acceleration of the Universe. This contribution reports on a mean field description of the backreaction in terms of a minimally coupled regionally homogeneous scalar field evolving in a potential, then giving a physical origin to the various phenomenological scalar fields generically called quintessence fields. As an example, the correspondence is then applied to scaling solutions.

This paper focusses on the barely understood gap between the weakly nonlinear regime of structure formation and the onset of the virialized regime. While the former is accessed through perturbative calculations and the latter through virialization conditions incorporating dynamical stresses that arise in collisionless self-gravitating systems due to velocity dispersion forces, the addressed regime can only be understood through non-perturbative models. We here present an exact Lagrangian integral that provides a tool to access this regime. We derive a transport equation for the peculiar-gravitational field strength and integrate it along comoving trajectories of fluid elements. The so-obtained integral provides an exact expression that solves the longitudinal gravitational field equation in general. We argue that this integral provides a powerful approximation beyond the Lagrangian perturbative regime, and discuss its relation to known approximations, among them Lagrangian perturbation solutions including the Zel'dovich approximation and approximations for adhesive gravitational clustering, including the adhesion approximation. Furthermore, we propose an iteration scheme for a systematic analytical and numerical construction of trajectory fields. The integral may also be employed to improve inverse reconstruction techniques.

An exact expression for the gravitational field strength in a self-gravitating dust continuum is derived within the Lagrangian picture of continuum mechanics. From the Euler-Newton system a transport equation for the gravitational field strength is formulated and then integrated along trajectories of continuum elements. The resulting integral solves one of the Lagrangian equations of the corresponding Lagrange-Newton system in general. Relations to known exact solutions without symmetry in Newtonian gravity are discussed. The presented integral may be employed to access the non-perturbative regime of structure formation in Newtonian cosmology, and to apply iterative Lagrangian schemes to solve the Lagrange-Newton system.

In the framework of spatially averaged inhomogeneous cosmologies in classical General Relativity, effective Einstein equations govern the regional and the global dynamics of averaged scalar variables of cosmological models. A particular solution may be characterized by a cosmic equation of state. In this paper it is pointed out that a globally static averaged dust model is conceivable without employing a compensating cosmological constant. Much in the spirit of Einstein's original model we discuss consequences for the global, but also for the regional properties of this cosmology. We then consider the wider class of globally stationary cosmologies that are conceivable in the presented framework. All these models are based on exact solutions of the averaged Einstein equations and provide examples of cosmologies in an out-of-equilibrium state, which we characterize by an information-theoretical measure. It is shown that such cosmologies preserve high-magnitude kinematical fluctuations and so tend to maintain their global properties. The same is true for a $\Lambda-$driven cosmos in such a state despite of exponential expansion. We outline relations to inflationary scenarios, and put the Dark Energy problem into perspective. Here, it is argued, on the grounds of the discussed cosmologies, that a classical explanation of Dark Energy through backreaction effects is theoretically conceivable, if the matter-dominated Universe emerged from a non-perturbative state in the vicinity of the stationary solution. We also discuss a number of caveats that furnish strong counter arguments in the framework of structure formation in a perturbed Friedmannian model.

A system of effective Einstein equations for spatially averaged scalar variables of inhomogeneous cosmological models can be solved by providing a `cosmic equation of state'. Recent efforts to explain Dark Energy focus on `backreaction effects' of inhomogeneities on the effective evolution of cosmological parameters in our Hubble volume, avoiding a cosmological constant in the equation of state. In this Letter it is argued that, if kinematical backreaction effects are indeed of the order of the averaged density (or larger as needed for an accelerating domain of the Universe), then the state of our regional Hubble volume would have to be in the vicinity of a far-from-equilibrium state that balances kinematical backreaction and average density. This property, if interpreted globally, is shared by a stationary cosmos with effective equation of state $p_{\rm eff} = -1/3 \rho_{\rm eff}$. It is concluded that a confirmed explanation of Dark Energy by kinematical backreaction may imply a paradigmatic change of cosmology.

A large-scale smoothed-out model of the universe ignores small-scale inhomogeneities, but the averaged effects of those inhomogeneities may alter both observational and dynamical relations at the larger scale. This article discusses these effects, and comments briefly on the relation to gravitational entropy.

The notion of `adhesion' has been advanced for the phenomenon of stabilization of large-scale structure emerging from gravitational instability of a cold medium. Recently, the physical origin of adhesion has been identified: a systematic derivation of the equations of motion for the density and the velocity fields leads naturally to the key equation of the `adhesion approximation' - however, under a set of strongly simplifying assumptions. In this work, we provide an evaluation of the current status of adhesive gravitational clustering and a clear explanation of the assumptions involved. Furthermore, we propose systematic generalizations with the aim to relax some of the simplifying assumptions. We start from the general Newtonian evolution equations for self-gravitating particles on an expanding Friedmann background and recover the popular `dust model' (pressureless fluid), which breaks down after the formation of density singularities; then we investigate, in a unified framework, two other models which, under the restrictions referred to above, lead to the `adhesion approximation'. We apply the Eulerian and Lagrangian perturbative expansions to these new models and, finally, we discuss some non-perturbative results that may serve as starting points for workable approximations of non-linear structure formation in the multi-stream regime. In particular, we propose a new approximation that includes, in limiting cases, the standard `adhesion model' and the Eulerian as well as Lagrangian first-order approximations.

We present a first morphometric investigation of a preliminary sample from the SDSS of 154287 galaxies with apparent magnitude 14.5<m_r<17.5 and redshift 0.001<z<0.4. We measure the Minkowski Functionals, which are a complete set of morphological descriptors. To account for the complicated wedge--like geometry of the present survey data, we construct isodensity contour surfaces from the galaxy positions in redshift space and employ two complementary methods of computing the Minkowski Functionals. We find that the observed Minkowski Functionals for SDSS galaxies are consistent with the prediction of a Lambda--dominated spatially--flat Cold Dark Matter model with random--Gaussian initial conditions, within the cosmic variance estimated from the corresponding mock catalogue. We expect that future releases of the SDSS survey will allow us to distinguish morphological differences in the galaxy distribution with regard to different morphological type and luminosity ranges.

The effective evolution of an inhomogeneous cosmological model may be described in terms of spatially averaged variables. We point out that in this context, quite naturally, a measure arises which is identical to a fluid model of the `Kullback-Leibler Relative Information Entropy', expressing the distinguishability of the local inhomogeneous mass density field from its spatial average on arbitrary compact domains. We discuss the time-evolution of `effective information' and explore some implications. We conjecture that the information content of the Universe -- measured by Relative Information Entropy of a cosmological model containing dust matter -- is increasing.

We discuss the relation between `bare' cosmological parameters as the true spatial average characteristics that determine the cosmological model, and the parameters interpreted by observers with a `Friedmannian bias', i.e., within a homogeneous space geometry. We may say that the latter are `dressed' by the smoothed-out geometrical inhomogeneities of the surveyed spatial region. We identify two effects that quantify the difference between `bare' and `dressed' parameters: `curvature backreaction' and `volume effect'. An estimate of the latter is given in terms of a simple geometrical example.

We numerically investigate the long--time evolution of density perturbations after the first appearance of caustics in an expanding cosmological model with one--dimensional `single--wave' initial conditions. Focussing on the time--intervals of caustic appearances and the spatial distribution of caustics at subsequent times, we find that the time--intervals of caustic appearances approach a constant, i.e., their time--subsequent ratio converges to 1; it is also found that the spatial distribution of caustics at a given time features some universality rules, e.g., the ratio between the position of the nearest caustic from the center and that of the second nearest caustic from the center approaches a constant. Furthermore we find some rules for the mass distribution for each caustic. Using these universality constants we are in the position to predict the spatial distribution of caustics at an arbitrary time in order to give an estimate for the power spectral index in the fully--developed non--dissipative turbulent (`virialized') regime.

Averaged inhomogeneous cosmologies lie at the forefront of interest, since cosmological parameters like the rate of expansion or the mass density are to be considered as volume-averaged quantities and only these can be compared with observations. For this reason the relevant parameters are intrinsically scale-dependent and one wishes to control this dependence without restricting the cosmological model by unphysical assumptions. In the latter respect we contrast our way to approach the averaging problem in relativistic cosmology with shortcomings of averaged Newtonian models. Explicitly, we investigate the scale-dependence of Eulerian volume averages of scalar functions on Riemannian three-manifolds. We propose a complementary view of a Lagrangian smoothing of (tensorial) variables as opposed to their Eulerian averaging on spatial domains. This program is realized with the help of a global Ricci deformation flow for the metric. We explain rigorously the origin of the Ricci flow which, on heuristic grounds, has already been suggested as a possible candidate for smoothing the initial data set for cosmological spacetimes. The smoothing of geometry implies a renormalization of averaged spatial variables. We discuss the results in terms of effective cosmological parameters that would be assigned to the smoothed cosmological spacetime.

In the context of the averaging problem in relativistic cosmology, we provide a key to the interpretation of cosmological parameters by taking into account the actual inhomogeneous geometry of the Universe. We discuss the relation between `bare' cosmological parameters determining the cosmological model, and the parameters interpreted by observers with a ``Friedmannian bias'', which are `dressed' by the smoothed-out geometrical inhomogeneities of the surveyed spatial region.

We explore the morphological and dynamical evolution of galaxy clusters in simulations using scalar and vector-valued Minkowski valuations and the concept of fundamental plane relations. In this context, three questions are of fundamental interest: 1. How does the average cluster morphology depend on the cosmological background model? 2. Is it possible to discriminate between different cosmological models using cluster substructure in a statistically significant way? 3. How is the dynamical state of a cluster, especially its distance from a virial equilibrium, correlated to its visual substructure? To answer these questions, we quantify cluster substructure using a set of morphological order parameters constructed on the basis of the Minkowski valuations (MVs). The dynamical state of a cluster is described using global cluster parameters: in certain spaces of such parameters fundamental band-like structures are forming indicating the emergence of a virial equilibrium. We find that the average distances from these fundamental structures are correlated to the average amount of cluster substructure for our cluster samples during the time evolution. Furthermore, significant differences show up between the high- and the low-Omega models. We pay special attention to the redshift evolution of morphological characteristics and find large differences between the cosmological models even for higher redshifts.

The redshift dependence of the abundance of collapsed objects places strong constraints on cosmological models of structure formation. We apply a recently proposed model describing the anisotropic collapse of inhomogeneous spatial domains. Compared with the spherical top-hat model, this generic model leads to significantly more collapsed objects at high redshifts: at redshift one and on the scale of rich clusters a factor of 65. Furthermore, for a fixed normalization of the initial fluctuation spectrum (\sigma_8=1), we predict four times as much presently collapsed objects on the mass-scale of rich clusters within the standard CDM cosmogony, compared to the spherical collapse.

For general relativistic spacetimes filled with an irrotational perfect fluid a generalized form of Friedmann's equations governing the expansion factor of spatially averaged portions of inhomogeneous cosmologies is derived. The averaging problem for scalar quantities is condensed into the problem of finding an `effective equation of state' including kinematical as well as dynamical `backreaction' terms that measure the departure from a standard FLRW cosmology. Applications of the averaged models are outlined including radiation-dominated and scalar field cosmologies (inflationary and dilaton/string cosmologies). In particular, the averaged equations show that the averaged scalar curvature must generically change in the course of structure formation, that an averaged inhomogeneous radiation cosmos does not follow the evolution of the standard homogeneous-isotropic model, and that an averaged inhomogeneous perfect fluid features kinematical `backreaction' terms that, in some cases, act like a free scalar field source. The free scalar field (dilaton) itself, modelled by a `stiff' fluid, is singled out as a special inhomogeneous case where the averaged equations assume a simple form.

In a follow-up study to a previous analysis of the IRAS 1.2Jy catalogue, we quantify the morphological fluctuations in the PSCz survey. We use a variety of measures, among them the family of scalar Minkowski functionals. We confirm the existence of significant fluctuations that are discernible in volume-limited samples out to 200Mpc/h. In contrast to earlier findings, comparisons with cosmological N-body simulations reveal that the observed fluctuations roughly agree with the cosmic variance found in corresponding mock samples. While two-point measures, e.g. the variance of count-in-cells, fluctuate only mildly, the fluctuations in the morphology on large scales indicate the presence of coherent structures that are at least as large as the sample.

Minkowski functionals constitute a family of order parameters which discriminate spatial patterns according to size, shape and connectivity. Here we point out, that these scalar descriptors can be complemented by vector-valued curvature measures also known as Quermass vectors. Using examples of galaxy clusters, we demonstrate that the Quermass vectors provide additional morphological information on directional features and symmetries displayed by spatial data.

Averaging and evolving inhomogeneities are non-commuting operations. This implies the existence of deviations of an averaged model from the standard Friedmann-Lemaitre cosmologies. We quantify these deviations, encoded in a backreaction parameter, in the framework of Newtonian cosmology. We employ the linear theory of gravitational instability in the Eulerian and Lagrangian approaches, as well as the spherically- and plane-symmetric solutions as standards of reference. We propose a model for the evolution of the average characteristics of a spatial domain for generic initial conditions that contains the spherical top-hat model and the planar collapse model as exact sub cases. A central result is that the backreaction term itself, calculated on sufficiently large domains, is small but, still, its presence can drive the cosmological parameters on the averaging domain far away from their global values of the standard model. We quantify the variations of these parameters in terms of the fluctuations in the initial data as derived from the power spectrum of initial cold dark matter density fluctuations. E.g. in a domain with a radius of 100Mpc today and initially one-sigma fluctuations, the density parameters deviate from their homogeneous values by 15%; three-sigma fluctuations lead to deviations larger than 100%.

The present talk summarizes the recently progressed state of a systematic re-evaluation of cosmological models that respect the presence of inhomogeneities. Emphasis is given to identifying the basic steps towards an effective (i.e. spatially averaged) description of structural evolution, also unfolding the various facets of a `smoothed--out' cosmology. We shall highlight some results obtained within Newtonian cosmology, discuss expansion laws in general relativity within a covariant fluid approach, and put forward some promising directions of future research.

For general relativistic spacetimes filled with irrotational `dust' a generalized form of Friedmann's equations for an `effective' expansion factor $a_D (t)$ of inhomogeneous cosmologies is derived. Contrary to the standard Friedmann equations, which hold for homogeneous-isotropic cosmologies, the new equations include the `backreaction effect' of inhomogeneities on the average expansion of the model. A universal relation between `backreaction' and average scalar curvature is also given. For cosmologies whose averaged spatial scalar curvature is proportional to $a_D^{-2}$, the expansion law governing a generic domain can be found. However, as the general equations show, `backreaction' acts as to produce average curvature in the course of structure formation, even when starting with space sections that are spatially flat on average.

The interplay between gravitational and dispersive forces in a multi-streamed medium leads to an effect which is exposed in the present note as the genuine driving force of stabilization of large-scale structure. The conception of `adhesive gravitational clustering' is advanced to interlock the fairly well-understood epoch of formation of large-scale structure and the onset of virialization into objects that are dynamically in equilibrium with their large-scale structure environment. The classical `adhesion model' is opposed to a class of more general models traced from the physical origin of adhesion in kinetic theory.

We apply Minkowski functionals and various derived measures to decipher the morphological properties of large-scale structure seen in simulations of gravitational evolution. Minkowski functionals of isodensity contours serve as tools to test global properties of the density field. Furthermore, we identify coherent objects at various threshold levels and calculate their partial Minkowski functionals. We propose a set of two derived dimensionless quantities, planarity and filamentarity, which reduce the morphological information in a simple and intuitive way. Several simulations of the gravitational evolution of initial power-law spectra provide a framework for systematic tests of our method.

Most cosmological models studied today are based on the assumption of homogeneity and isotropy. Observationally one can find evidence that supports these assumptions on very large scales, the strongest being the almost isotropy of the Cosmic Microwave Background radiation after assigning the whole dipole to our proper motion relative to this background. However, on small and on intermediate scales up to several hundreds of Mpcs, there are strong deviations from homogeneity and isotropy. Here the problem arises how to relate the observations with the homogeneous and isotropic models. The usual proposal for solving this problem is to assume that Friedmann-Lemaitre models describe the mean observables. Such mean values may be identified with spatial averages. For Newtonian fluid dynamics the averaging procedure has been discussed in detail in Buchert and Ehlers (1997), leading to an additional backreaction term in the Friedmann equation. We use the Eulerian linear approximation and the `Zel'dovich approximation' to estimate the effect of the backreaction term on the expansion. Our results indicate that even for domains matching the background density in the mean, the evolution of the scale factor strongly deviates from the Friedmann solution, critically depending on the velocity field inside.

We propose a phenomenological generalization of the models of large-scale structure formation in the Universe by gravitational instability in two ways: we include pressure forces to model multi-streaming, and noise to model fluctuations due to neglected short-scale physical processes. We show that pressure gives rise to a viscous-like force of the same character as that one introduced in the ``adhesion model'', while noise leads to a roughening of the density field yielding a scaling behavior of its correlations.

We discuss implications of the fundamental plane parameters of clusters of galaxies derived from combined optical and X-ray data of a sample of 78 nearby clusters. In particular, we investigate the dependence of these parameters on the dynamical state of the cluster. We introduce a new concept of allocation of the fundamental plane of clusters derived from their intrinsic morphological properties, and put some theoretical implications of the existence of a fundamental plane into perspective.

The Lagrangian theory of structure formation in cosmological fluids, restricted to the matter model ``dust'', provides successful models of large-scale structure in the Universe in the laminar regime, i.e., where the fluid flow is single-streamed and ``dust''-shells are smooth. Beyond the epoch of shell-crossing a qualitatively different behavior is expected, since in general anisotropic stresses powered by multi-stream forces arise in collisionless matter. In this paper we provide the basic framework for the modeling of pressure-supported fluids, thus restricting attention to isotropic stresses and to the cases where pressure can be given as a function of the density. We derive the governing set of Lagrangian evolution equations and study the resulting system using Lagrangian perturbation theory. We discuss the first-order equations and compare them to the Eulerian theory of gravitational instability, as well as to the case of plane-symmetric collapse. We obtain a construction rule that allows to derive first-order solutions of the Lagrangian theory from known first-order solutions of the Eulerian theory and so extend Zel'dovich's extrapolation idea into the multi-streamed regime. These solutions can be used to generalize current structure formation models in the spirit of the ``adhesion approximation''.

We present the function J as a morphological descriptor for point patterns formed by the distribution of galaxies in the Universe. This function was recently introduced in the field of spatial statistics, and is based on the nearest neighbor distribution and the void probability function. The J descriptor allows to distinguish clustered (i.e. correlated) from ``regular'' (i.e. anti-correlated) point distributions. We outline the theoretical foundations of the method, perform tests with a Matern cluster process as an idealised model of galaxy clustering, and apply the descriptor to galaxies and loose groups in the Perseus-Pisces Survey. A comparison with mock-samples extracted from a mixed dark matter simulation shows that the J descriptor can be profitably used to constrain (in this case reject) viable models of cosmic structure formation.

The generally held view that a model of large-scale structure, formed by collisionless matter in the Universe, can be based on the matter model ``dust'' fails in the presence of multi-stream flow, i.e., velocity dispersion. We argue that models for large-scale structure should rather be constructed for a flow which describes the average motion of a multi-stream system. We present a clearcut reasoning how to approach the problem and derive an evolution equation for the mean peculiar-velocity relative to background solutions of Friedmann-Lema\^\i tre type. We consider restrictions of the nonlinear problem and show that the effect of velocity dispersion gives rise to an effective viscosity of non-dissipative gravitational origin. We discuss subcases which arise naturally from this approach: the ``sticky particle model'' and the ``adhesion approximation''. We also construct a novel approximation that features adhesive action in the multi-stream regime while conserving momentum, which was considered a drawback of the standard approximation based on Burger's equation. We finally argue that the assumptions made to obtain these models should be relaxed and we discuss how this can be achieved.

An analysis of the IRAS 1.2 Jy redshift catalogue with emphasis on the separate examination of northern and southern parts (in galactic coordinates) is performed using a complete set of morphological descriptors (Minkowski functionals), nearest neighbour distributions and the variance of the galaxy counts. We find large fluctuations in the clustering properties as seen in a large difference between the northern and southern parts of the catalogue on scales of 100Mpc/h. These fluctuations remain discernible even on the scale of 200Mpc/h. We also identify sparse sampling as a major source of ``apparent homogenization''. Tests on observational selection effects concerning luminosity, colour and redshift-space distortion support the significance of these large-scale fluctuations.

I report on a systematic derivation of the phenomenological ``adhesion approximation'' from gravitational instability together with a brief evaluation of the related status of analytical modeling of large-scale structure.

The morphology of galaxy clusters is quantified using Minkowski functionals, especially the vector-valued ones, which contain directional information and are related to curvature centroids. The asymmetry of clusters and the amount of their substructure can be characterized in a unique way using these measures. -- We briefly introduce vector-valued Minkowski functionals (also known as Querma\ss vectors) and suggest their application to cluster data in terms of a morphological characterization of excursion sets. Furthermore, we develop robust structure functions which describe the dynamical state of a cluster and study the evolution of clusters using numerical simulations.

A first- and second-order relation between cosmic density and peculiar-velocity fields is presented. The calculation is purely Lagrangian and it is derived using the second-order solutions of the Lagrange-Newton system obtained by Buchert & Ehlers. The procedure is applied to two particular solutions given generic initial conditions. In this approach, the continuity equation yields a relation between the over-density and peculiar-velocity fields that automatically satisfies Euler's equation because the orbits are derived from the Lagrange-Newton system. This scheme generalizes some results obtained by Nusser et al. (1991) in the context of the Zel'dovich approximation. As opposed to several other reconstruction schemes, in this approach it is not necessary to truncate the expansion of the Jacobian given by the continuity equation in order to calculate a first- or second-order expression for the density field. In these previous schemes, the density contrast given by (a) the continuity equation and (b) Euler's equation are mutually incompatible. This inconsistency arises as a consequence of an improper handling of Lagrangian and Eulerian coordinates in the analysis. Here, we take into account the fact that an exact calculation of the density is feasible in the Lagrangian picture and therefore an accurate and consistent description is obtained.

The averaging problem for inhomogeneous cosmologies is discussed in the form of a disputation between two cosmologists, one of them (RED) advocating the standard model, the other (GREEN) advancing some arguments against it. Technical explanations of these arguments as well as the conclusions of this debate are given by BLUE.

In an analysis of the IRAS 1.2 Jy redshift catalogue, with emphasis on the separate examination of northern and southern parts (in galactic coordinates), we found that the clustering of galaxies differs significantly between north and south, showing fluctuations in the clustering properties at least on scales of 100 Mpc/h (Kerscher et al. 1997, astro-ph/9704028). We give a brief description of our morphological method which is based on Minkowski functionals and show the results obtained from the IRAS 1.2 Jy galaxy catalogue. We discuss several error estimates and select different subsamples from the 1.2 Jy catalogue according to flux and colour and validate the results for the whole sample. Furthermore we look closer at the spatial origin of the fluctuation and we compare with optically selected galaxies from the CfA1 survey.

The genus statistics of isodensity contours has become a well-established tool in cosmology. In this Letter we place the genus in the wider framework of a complete family of morphological descriptors. These are known as the Minkowski functionals, and we here apply them for the first time to isodensity contours of a continuous random field. By taking two equivalent approaches, one through differential geometry, the other through integral geometry, we derive two complementary formulae suitable for numerically calculating the Minkowski functionals. As an example we apply them to simulated Gaussian random fields and compare the outcome to the analytically known results, demonstrating that both are indeed well suited for numerical evaluation. The code used for calculating all Minkowski functionals is available from the authors.

The Lagrangian perturbation theory on Friedmann-Lemaitre cosmologies is compared with numerical simulations (tree-, adaptive P$^3$M- and PM codes). In previous work we have probed the large-scale performance of the Lagrangian perturbation solutions up to the third order by studying their cross- correlations with N-body simulations for various power spectra (Buchert etal 1994, Melott etal 1995, Weiss etal 1996). Thereby, spatial optimization techniques were applied by (high-frequency-)filtering of the initial power spectra. In this work the novel method of temporal optimization [Shifted-Time- Approximation (STA) and Frozen-Time-Approximation (FTA)] is investigated and used. The method is designed to compensate the native property of Lagrangian perturbation solutions to delay the collapse of structures. The method can be treated analytically. Applying the STA and FTA prescriptions a significant improvement of the performance of Lagrangian perturbation schemes (as measured by cross-correlation statistics) is observed. Using this tool we investigate a local study of special clustering models of dark matter as candidates for typical elements of the large-scale structure in the Universe, and so also focus on the performance of the perturbation solutions on smaller scales at high-spatial resolution. The models analyzed were presented in (Buchert etal 1996) and allow studying typical features of the clustering process in the weakly non-linear regime. The spatial and temporal limits of applicability of the solutions at second and third order are determined and compared with the first-order solution, which is equivalent to the ``Zel'dovich approximation'' (Zel'dovich 1970, 1973) for the type of initial data analyzed.

Idealizing matter as a pressureless fluid and representing its motion by a peculiar--velocity field superimposed on a homogeneous and isotropic Hubble expansion, we apply (Lagrangian) spatial averaging on an arbitrary domain $\cal D$ to the (nonlinear) equations of Newtonian cosmology and derive an exact, general equation for the evolution of the (domain dependent) scale factor $a_{\cal D}(t)$. We consider the effect of inhomogeneities on the average expansion and discuss under which circumstances the standard description of the average motion in terms of Friedmann's equation holds. We find that this effect vanishes for spatially compact models if one averages over the whole space. For spatially infinite inhomogeneous models obeying the cosmological principle of large--scale isotropy and homogeneity, Friedmann models may provide an approximation to the average motion on the largest scales, whereas for hierarchical (Charlier--type) models the general expansion equation shows how inhomogeneities might appreciably affect the expansion at all scales. An averaged vorticity evolution law is also given. Since we employ spatial averaging, the problem of justifying ensemble averaging does not arise. A generalization of the expansion law to general relativity is straightforward for the case of irrotational flows and will be discussed. The effect may have important consequences for a variety of problems in large--scale structure modeling as well as for the interpretation of observations.

The ``Newtonian'' theory of spatially unbounded, self--gravitating, pressureless continua in Lagrangian form is reconsidered. Following a review of the pertinent kinematics, we present alternative formulations of the Lagrangian evolution equations and establish conditions for the equivalence of the Lagrangian and Eulerian representations. We then distinguish open models based on Euclidean space $\R^3$ from closed models based (without loss of generality) on a flat torus $\T^3$. Using a simple averaging method we show that the spatially averaged variables of an inhomogeneous toroidal model form a spatially homogeneous ``background'' model and that the averages of open models, if they exist at all, in general do not obey the dynamical laws of homogeneous models. We then specialize to those inhomogeneous toroidal models whose (unique) backgrounds have a Hubble flow, and derive Lagrangian evolution equations which govern the (conformally rescaled) displacement of the inhomogeneous flow with respect to its homogeneous background. Finally, we set up an iteration scheme and prove that the resulting equations have unique solutions at any order for given initial data, while for open models there exist infinitely many different solutions for given data.

We determine the Minkowski functionals for a sample of Abell/ACO clusters, 401 with measured and 16 with estimated redshifts. The four Minkowski functionals (including the void probability function and the mean genus) deliver a global description of the spatial distribution of clusters on scales from $10$ to $60\hMpc$ with a clear geometric interpretation. Comparisons with mock catalogues of N--body simulations using different variants of the CDM model demonstrate the discriminative power of the description. The standard CDM model and the model with tilted perturbation spectrum cannot generate the Minkowski functionals of the cluster data, while a model with a cosmological constant and a model with breaking of the scale invariance of perturbations (BSI) yield compatible results.

Average properties of general inhomogeneous cosmological models are discussed in the Newtonian framework. It is shown under which circumstances the average flow reduces to a member of the standard Friedmann--Lema\^\i tre cosmologies. Possible choices of global boundary conditions of inhomogeneous cosmologies as well as consequences for the interpretation of cosmological parameters are put into perspective.

We discuss the application of Minkowski functionals to galaxy catalogues with emphasis on the boundary corrections and point out the relevance for other statistical methods. The analysis of a cubic subsample of the CfA--1 data serves as an example.

We report on the performance of Lagrangian perturbation theory up to the second order for the standard cold dark matter (SCDM) and broken scale invariance (BSI) scenarios. We normalize both models to the COBE data, the BSI model serves as an example of models which fit the small-scale power of galaxy surveys. We optimize Lagrangian perturbation solutions by removing small-scale power from the initial data and compare the results with those of numerical simulations. We find an excellent performance of the optimized Lagrangian schemes down to scales around the correlation length or smaller, depending on the statistics used for the comparison. The optimization scheme can be expressed in a way which is independent of the type of fluctuation spectrum and of the size of the simulations.

We present high--spatial resolution studies of the density field as predicted by Lagrangian perturbation approximations up to the third order. The first--order approximation is equivalent to the ``Zel'dovich approximation'' for the type of initial data analyzed. The study is performed for two simple models which allow studying of typical features of the clustering process in the early non--linear regime. We calculate the initial perturbation potentials as solutions of Poisson equations algebraically, and automate this calculation for a given initial random density field. The presented models may also be useful for other questions addressed to Lagrangian perturbation solutions and for the comparison of different approximation schemes. In an accompanying paper we investigate a detailed comparison with various N--body integrators using these models (Karakatsanis \& Buchert 1995). Results of the present paper include the following: 1. The collapse is accelerated significantly by the higher--order corrections confirming previous results by Moutarde \etal (1991); 2. the spatial structure of the density patterns predicted by the ``Zel'dovich approximation'' differs much from those predicted by the second-- and third--order Lagrangian approximations; 3. Second--order effects amount to internal substructures such as ``second generation'' --pancakes, --filaments and --clusters, as are also observed in N--body simulations; 4. The third--order effect gives rise to substructuring of the secondary mass--shells. The hierarchy of shell--crossing singularities that form features small high--density clumps at the intersections of caustics which we interprete as gravitational fragmentation.

Minkowski functionals provide a novel tool to characterize the large-scale galaxy distribution in the Universe. Here we give a brief tutorial on the basic features of these morphological measures and indicate their practical application for simulation data and galaxy redshift catalogues as examples.

A new method for the statistical analysis of 3D point processes, based on the family of Minkowski functionals, is explained and applied to modelled galaxy distributions generated by a toy-model and cosmological simulations of the large-scale structure in the Universe. These measures are sensitive to both, geometrical and topological properties of spatial patterns and appear to be very effective in discriminating different point processes. Moreover by the means of conditional subsampling, different building blocks of large-scale structures like sheets, filaments and clusters can be detected and extracted from a given distribution.

The present lecture notes address three columns on which the Lagrangian perturbation approach to cosmological dynamics is based: 1. the formulation of a Lagrangian theory of self--gravitating flows in which the dynamics is described in terms of a single field variable; 2. the procedure, how to obtain the dynamics of Eulerian fields from the Lagrangian picture, and 3. a precise definition of a Newtonian cosmology framework in which Lagrangian perturbation solutions can be studied. While the first is a discussion of the basic equations obtained by transforming the Eulerian evolution and field equations to the Lagrangian picture, the second exemplifies how the Lagrangian theory determines the evolution of Eulerian fields including kinematical variables like expansion, vorticity, as well as the shear and tidal tensors. The third column is based on a specification of initial and boundary conditions, and in particular on the identification of the average flow of an inhomogeneous cosmology with a ``Hubble--flow''. Here, we also look at the limits of the Lagrangian perturbation approach as inferred from comparisons with N--body simulations and illustrate some striking properties of the solutions.

We investigate the performance of Lagrangian perturbation theory up to the second order for two scenarios of cosmological large-scale structure formation, SCDM (standard cold dark matter) and BSI (broken scale invariance). The latter model we study as a representative of COBE-normalized CDM models which fit the small-scale power of galaxy surveys. In this context we optimize the performance of the Lagrangian perturbation schemes by smoothing the small-scale fluctuations in the initial data. The results of the so obtained Lagrangian mappings are computed for a set of COBE-normalized SCDM and BSI initial data of different sizes and at different times. We compare these results against those obtained with a numerical PM-code. We find an excellent performance of the optimized Lagrangian schemes down to scales close to the correlation length. This is explained by the counterintuitive fact that nonlinearities in the model can produce more small-scale power, if initially such power is removed. The optimization scheme can be expressed in a way which is independent of the type of fluctuation spectrum and of the size of the simulations.

Approximations to the exact solutions for gravitational instability in the expanding Universe are extremely useful for understanding the evolution of large--scale structure. We report on a series of tests of Newtonian Lagrangian perturbation schemes using N--body simulations for various power--spectra with scale--independent indices in the range $-3$ to $+1$. The models have been evolved deeply into the non--linear regime of structure formation in order to probe the dynamical and statistical performance of the Lagrangian perturbation schemes (whose first--order solution contains as a subset the celebrated ``Zel'dovich--approximation'', hereafter ZA). These tests reveal properties of the approximations at stages beyond the obvious validity of perturbation theory. Recently, another series of tests of different analytical and semi--numerical approximations for large--scale structure was conducted with the result that ZA displays the best dynamical performance in comparison with the N--body simulations, if the initial data were smoothed before evolving the model, i.e., a truncated form of ZA (TZA). We show in this contribution that the excellent performance of TZA can be further improved by going to second order in the Lagrangian perturbation approach. The truncated second--order Lagrangian scheme provides a useful improvement over TZA especially for negative power indices, which suggests it will be very useful for modelling standard scenarios such as ``Cold--'', ``Hot--'' and ``Mixed--Dark--Matter''.

The problem of formation of generic structures in the Universe is addressed, whereby first the kinematics of inertial continua for coherent initial data is considered. The generalization to self--gravitating continua is outlined focused on the classification problem of singularities and metamorphoses arising in the density field. Self--gravity gives rise to an internal hierarchy of structures, and, dropping the assumption of coherence, also to an external hierarchy of structures dependent on the initial power spectrum of fluctuations.

A complete family of statistical descriptors for the morphology of large--scale structure based on Minkowski--Functionals is presented. These robust and significant measures can be used to characterize the local and global morphology of spatial patterns formed by a coverage of point sets which represent galaxy samples. Basic properties of these measures are highlighted and their relation to the `genus statistics' is discussed. Test models like a Poissonian point process and samples generated from a Voronoi--model are put into perspective.

We report on a series of tests of Newtonian Lagrangian perturbation schemes using N--body simulations for various power--spectra with scale--independent indices in the range $-3$ to $+1$. The models have been evolved deeply into the non--linear regime of structure formation in order to probe the dynamical and statistical performance of the Lagrangian perturbation schemes (whose first--order solution contains as a subset the celebrated ``Zel'dovich--approximation'' (hereafter ZA). These tests reveal properties of the approximations at stages beyond the obvious validity of perturbation theory. Recently, another series of tests of different analytical and semi--numerical approximations for large--scale structure was conducted with the result that ZA displays the best dynamical performance in comparison with the N--body simulations, if the initial data were smoothed before evolving the model, i.e., a truncated form of ZA (TZA). We show in this Letter that the excellent performance of TZA can be further improved by going to second order in the Lagrangian perturbation approach. The truncated second--order Lagrangian scheme provides a useful improvement over TZA especially for negative power indices, which suggests it will be very useful for modelling standard scenarios such as ``Cold--'', ``Hot--'' and ``Mixed--Dark--Matter''.

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann--Lema\^\i tre cosmogonies investigated and solved up to the third order in the series of papers by Buchert (1989, 1992, 1993), Buchert \& Ehlers (1993), Buchert (1994), Ehlers \& Buchert (1994), is compared with numerical simulations. In this paper we study the dynamics of hierarchical models as a second step. In the first step (Buchert, Melott and Wei{\ss} 1994) we analyzed the performance of the Lagrangian schemes for pancake models, i.e., models which initially have a truncated power spectrum. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non--linear regime............We find that for spectra with negative power--index the second--order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low--order statistics like the power--spectrum. In cases with much small--scale power the gain from the higher--order schemes is small, but still measurable. However, in contrast to the results found for pancake models, where the higher--order schemes get worse than TZA at late non--linear stages and on small scales, we here find that the second--order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these

We report on a recent paper (Wei\ss\ \& Buchert 1993), where we carry out pencil beam constructions in a high-resolution simulation of the large-scale structure of galaxies. As an example we present the results for the case of ``Hot-Dark-Matter" (HDM) initial conditions (with scale-free $n = 1$ power index on large scales and $\Omega = 1$) as a representative of models with sufficient large-scale power. We use an analytic approximation for particle trajectories of a self-gravitating dust continuum and apply a local dynamical biasing of volume elements to identify luminous matter in the model. Using this method, we are able to resolve formally a simulation box of 1200h^-1 Mpc (e.g. for HDM initial conditions) down to the scale of galactic halos using $2160^3$ particles. Pencil beam probes are taken for a given epoch using the parameters of observed beams. In particular, our analysis concentrates on the detection of a quasi-periodicity in the beam probes. The simulation is designed for application to parameter studies which prepare future observational projects. We find that a large percentage of the beams shows quasi-periodicities with periods which cluster at a certain length scale. The periods found range between one and eight times the cutoff length in the initial fluctuation spectrum.

Fundamental assumptions which form the basis of models for large-scale structure in the Universe are sketched in light of a Lagrangian description of inhomogeneities. This description is introduced for Newtonian self-gravitating flows. On its basis a Lagrangian perturbation approach is discussed and compared with the standard Eulerian theory of gravitational instability. The performance of Lagrangian perturbation solutions up to the third order is demonstrated in comparison with numerical N-body simulations. First results of this comparison are presented for large scales (PM-code) and for small scales (tree-code).

We propose a novel method for the description of spatial patterns formed by a coverage of point sets representing galaxy samples. This method is based on a complete family of morphological measures known as Minkowski functionals, which includes the topological Euler characteristic and geometric descriptors to specify the content, shape and connectivity of spatial sets. The method is numerically robust even for small samples, independent of statistical assumptions, and yields global as well as local morphological information. We illustrate the method by applying it to a Poisson process, a `double-Poisson' process, and to the Abell catalogue of galaxy clusters.

The Lagrangian theory of gravitational instability of homogeneous-isotropic Friedman-Lemaitre cosmogonies investigated and solved in the series of papers by Buchert (1989), (1992), Buchert & Ehlers (1993), Buchert (1993a,b), Ehlers & Buchert (1993), is illustrated. The third-order solution of this theory for generic initial conditions is presented and realized in a special case by employing methods of high-spatial resolution of the density field.

The Lagrangian perturbation theory on Friedman-Lemaitre cosmologies investigated and solved up to the second order in earlier papers (Buchert 1992, Buchert \& Ehlers 1993) is evaluated up to the third order. On its basis a model for non-linear clustering applicable to the modeling of large-scale structure in the Universe for generic initial conditions is formulated. A truncated model is proposed which represents the ``main body'' of the perturbation sequence in the early non-linear regime by neglecting all gravitational sources which describe interaction of the perturbations. However, I also give the irrotational solutions generated by the interaction terms to the third order, which induce vorticity in Lagrangian space. The consequences and applicability of the solutions are put into perspective. In particular, the model presented enables the study of previrialization effects in gravitational clustering and the onset of non-dissipative gravitational turbulence within the cluster environment.

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasilinear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order in the series of papers by Buchert (1989, 1992, 1993a), Buchert \& Ehlers (1993), Buchert (1993b), Ehlers \& Buchert (1993), is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work (Coles \etal 1993, Melott \etal 1993, Melott 1993) the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the ``Zel'dovich approximation'' (Zel'dovich 1970, 1973, hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of $\sigma \approx 2$). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-index $n=-1$) using cross-correlation statistics employed in