A number of codes for general-relativistic simulations of cosmological structure formation have been developed in recent years. Here we demonstrate that a sample of these codes produce consistent results beyond the Newtonian regime. We simulate solutions to Einstein's equations dominated by gravitomagnetism -- a vector-type gravitational field that doesn't exist in Newtonian gravity and produces frame-dragging, the leading-order post-Newtonian effect. We calculate the coordinate-invariant effect on intersecting null geodesics by performing ray tracing in each independent code. With this observable quantity, we assess and compare each code's ability to compute relativistic effects.

We perform a reconstruction of the coupling function between vacuum energy and geodesic cold dark matter using the latest observational data. We bin the interaction in seventeen redshift bins but use a correlation prior to prevent rapid, unphysical oscillations in the coupling function. This prior also serves to eliminate any dependence of the reconstruction on the binning method. We use two different forms of the correlation prior, finding that both give similar results for the reconstruction of the dark matter -- dark energy interaction. Calculating the Bayes factor for each case, we find no meaningful evidence for deviation from the null interacting case, i.e. $\Lambda$CDM, in our reconstruction.

We investigate an interacting dark sector scenario in which the vacuum energy is free to interact with cold dark matter (CDM), which itself is assumed to cluster under the sole action of gravity, i.e. it is in free fall (geodesic), as in $\Lambda$CDM. The interaction is characterised by a dimensionless coupling $q_{\rm V}$ that we constrain using cosmic microwave background data from the Planck 2015 data release, along with baryon acoustic oscillation, redshift space distortion and Type Ia supernova measurements. We present the full linear perturbation theory of this interacting scenario and use MCMC sampling to study five different cases: two cases in which we have $\Lambda$CDM evolution in the distant past, until a set redshift $z_{\rm trans}$, below which the interaction switches on and $q_{\rm V}$ is the single sampled parameter, with $z_{\rm trans}$ fixed at $z_{\rm trans}=3000$ and $z_{\rm trans}=0.9$ respectively; a case where we allow this transition redshift to vary along with $q_{\rm V}$; a case in which the vacuum energy is zero for $z>z_{\rm trans}$ and then begins to grow once the interaction switches on; and the final case in which we bin $q_{\rm V}(z)$ in four redshift bins to investigate the possibility of a dynamical interaction, reconstructing the redshift evolution of the function using Gaussian processes. We find that, in all cases where the high redshift evolution is not modified, the results are compatible with a vanishing coupling, thus finding no significant deviation from $\Lambda$CDM.

A Bianchi IX Mixmaster spacetime is the most general spatially homogeneous solution of Einstein's equations and it can represent the space-averaged Universe. We introduce two novel mechanisms resulting in a Mixmaster Universe with non-singular bounces which are quasi-isotropic. A fluid with a non-linear equation of state allows non-singular bounces. Using negative anisotropic stresses successfully isotropises this Universe and mitigates the well known Mixmaster chaotic behaviour. Thus the Universe can be an eternal Mixmaster, going through an infinite series of different cycles separated by bounces, with a sizable fraction of cycles isotropic enough to be well approximated by a standard Friedmann-Lema\^itre-Robertson-Walker model from the radiation era onward.

We present a full-sky derivation of weak lensing observables in the Post-Friedmann (PF) formalism. Weak lensing has the characteristic of mixing small scales and large scales since it is affected by inhomogeneities integrated along the photon trajectory. With the PF formalism, we develop a modelling of lensing observables which encompasses both leading order relativistic effects and effects that are due to the fully non-linear matter distribution at small scales. We derive the reduced shear, convergence and rotation up to order $1/c^4$ in the PF approximation, accounting for scalar, vector and tensor perturbations, as well as galaxies' peculiar velocities. We discuss the various contributions that break the Kaiser-Squires relation between the shear and the convergence at different orders. We pay particular attention to the impact of the frame-dragging vector potential on lensing observables and we discuss potential ways to measure this effect in future lensing surveys.

In this paper we examine how primordial non-Gaussianity contributes to nonlinear perturbative orders in the expansion of the density field at large scales in the matter dominated era. General Relativity is an intrinsically nonlinear theory, establishing a nonlinear relation between the metric and the density field. Representing the metric perturbations with the curvature perturbation zeta, it is known that nonlinearity produces effective non-Gaussian terms in the nonlinear perturbations of the matter density field, even if the primordial zeta is Gaussian. Here we generalise these results to the case of a non-Gaussian primordial zeta. Using a standard parametrization of primordial non-Gaussianity in zeta in terms of fNL, gNL, hNL..., we show how at higher order (from third and higher) nonlinearity also produces a mixing of these contributions to the density field at large scales, e.g. both fNL and gNL contribute to the third order in the density contrast. This is the main result of this paper. Our analysis is based on the synergy between a gradient expansion (aka long-wavelength approximation) and standard perturbation theory at higher order. In essence, mathematically the equations for the gradient expansion are equivalent to those of first order perturbation theory, thus first-order results convert into gradient expansion results and, vice versa, the gradient expansion can be used to derive results in perturbation theory at higher order and large scales.

We study the nonlinear gravitational dynamics of a universe filled with a pressureless fluid and a cosmological constant $\Lambda$ in the context of Newtonian gravity, and in the relativistic post-Friedmann approach proposed in paper I [I. Milillo et al., Phys. Rev. D 92, 023519 (2015).]. The post-Friedmann approximation scheme is based on the $1/c$ expansion of the space-time metric and the energy-momentum tensor, and includes nonlinear Newtonian cosmology. Here we establish the nonlinear post-Friedmann framework in the Lagrangian-coordinates approach for structure formation. For this we first identify a Lagrangian gauge which is suitable for incorporating nonzero vorticity. We analyze our results in two limits: at the leading order we recover the fully nonlinear Newtonian cosmological equations in the Lagrangian formulation, and we provide a space-time metric consistent from the perspective of general relativity. We then linearize our expressions and recover the relativistic results at first order in cosmological perturbation theory. Therefore, the introduced approximation scheme provides a unified treatment for the two leading-order regimes, from the small scales described by Newtonian gravity to the large linear scale, where first-order relativistic cosmological perturbation theory gives a very good description of structure formation.

We construct a three-dimensional, fully relativistic numerical model of a universe filled with an inhomogeneous pressureless fluid, starting from initial data that represent a perturbation of the Einstein-de Sitter model. We then measure the departure of the average expansion rate with respect to this homogeneous and isotropic reference model, comparing local quantities to the predictions of linear perturbation theory. We find that collapsing perturbations reach the turnaround point much earlier than expected from the reference spherical top-hat collapse model and that the local deviation of the expansion rate from the homogeneous one can be as high as $28\%$ at an underdensity, for an initial density contrast of $10^{-2}$. We then study, for the first time, the exact behavior of the backreaction term ${\cal Q}_{\cal D}$. We find that, for small values of the initial perturbations, this term exhibits a $1/a$ scaling, and that it is negative with a linearly growing absolute value for larger perturbation amplitudes, thereby contributing to an overall deceleration of the expansion. Its magnitude, on the other hand, remains very small even for relatively large perturbations.

Newtonian simulations are routinely used to examine the matter dynamics on non-linear scales. However, even on these scales, Newtonian gravity is not a complete description of gravitational effects. A post-Friedmann approach shows that the leading order correction to Newtonian theory is a vector potential in the metric. This vector potential can be calculated from N-body simulations, requiring a method for extracting the velocity field. Here, we present the full details of our calculation of the post-Friedmann vector potential, using the Delauney Tesselation Field Estimator (DTFE) code. We include a detailed examination of the robustness of our numerical result, including the effects of box size and mass resolution on the extracted fields. We present the power spectrum of the vector potential and find that the power spectrum of the vector potential is $\sim 10^5$ times smaller than the power spectrum of the fully non-linear scalar gravitational potential at redshift zero. Comparing our numerical results to perturbative estimates, we find that the fully non-linear result can be more than an order of magnitude larger than the perturbative estimate on small scales. We extend the analysis of the vector potential to multiple redshifts, showing that this ratio persists over a range of scales and redshifts. We also comment on the implications of our results for the validity and interpretation of Newtonian simulations.

In General Relativity, the constraint equation relating metric and density perturbations is inherently nonlinear, leading to an effective non-Gaussianity in the dark matter density field on large scales - even if the primordial metric perturbation is Gaussian. Intrinsic non-Gaussianity in the large-scale dark matter overdensity in GR is real and physical. However, the variance smoothed on a local physical scale is not correlated with the large-scale curvature perturbation, so that there is no relativistic signature in the galaxy bias when using the simplest model of bias. It is an open question whether the observable mass proxies such as luminosity or weak lensing correspond directly to the physical mass in the simple halo bias model. If not, there may be observables that encode this relativistic signature.

In this paper we examine cosmological weak lensing on non-linear scales and show that there are Newtonian and relativistic contributions and that the latter can also be extracted from standard Newtonian simulations. We use the post-Friedmann formalism, a post-Newtonian type framework for cosmology, to derive the full weak-lensing deflection angle valid on non-linear scales for any metric theory of gravity. We show that the only contributing term that is quadratic in the first order deflection is the expected Born correction and lens-lens coupling term. We use this deflection angle to analyse the vector and tensor contributions to the E- and B- mode cosmic shear power spectra. In our approach, once the gravitational theory has been specified, the metric components are related to the matter content in a well-defined manner. Specifying General Relativity, we write down a complete set of equations for a GR$+\Lambda$CDM universe for computing all of the possible lensing terms from Newtonian N-body simulations. We illustrate this with the vector potential and show that, in a GR$+\Lambda$CDM universe, its contribution to the E-mode is negligible with respect to that of the conventional Newtonian scalar potential, even on non-linear scales. Thus, under the standard assumption that Newtonian N-body simulations give a good approximation of the matter dynamics, we show that the standard ray tracing approach gives a good description for a $\Lambda$CDM cosmology.

Many modified gravity theories are under consideration in cosmology as the source of the accelerated expansion of the universe and linear perturbation theory, valid on the largest scales, has been examined in many of these models. However, smaller non-linear scales offer a richer phenomenology with which to constrain modified gravity theories. Here, we consider the Hu-Sawicki form of $f(R)$ gravity and apply the post-Friedmann approach to derive the leading order equations for non-linear scales, i.e. the equations valid in the Newtonian-like regime. We reproduce the standard equations for the scalar field, gravitational slip and the modified Poisson equation in a coherent framework. In addition, we derive the equation for the leading order correction to the Newtonian regime, the vector potential. We measure this vector potential from $f(R)$ N-body simulations at redshift zero and one, for two values of the $f_{R_0}$ parameter. We find that the vector potential at redshift zero in $f(R)$ gravity can be close to 50\% larger than in GR on small scales for $|f_{R_0}|=1.289\times10^{-5}$, although this is less for larger scales, earlier times and smaller values of the $f_{R_0}$ parameter. Similarly to in GR, the small amplitude of this vector potential suggests that the Newtonian approximation is highly accurate for $f(R)$ gravity, and also that the non-linear cosmological behaviour of $f(R)$ gravity can be completely described by just the scalar potentials and the $f(R)$ field.

We discuss the question of gauge choice when analysing relativistic density perturbations at second order. We compare Newtonian and General Relativistic approaches. Some misconceptions in the recent literature are addressed. We show that the comoving-synchronous gauge is the unique gauge in General Relativity that corresponds to the Lagrangian frame and is entirely appropriate to describe the matter overdensity at second order. The comoving-synchronous gauge is the simplest gauge in which to describe Lagrangian bias at second order.

We present a nonlinear post-Friedmann framework for structure formation, generalizing to cosmology the weak-field (post-Minkowskian) approximation, unifying the treatment of small and large scales. We consider a universe filled with a pressureless fluid and a cosmological constant $\Lambda$, the theory of gravity is Einstein's general relativity and the background is the standard flat $\Lambda$CDM cosmological model. We expand the metric and the energy-momentum tensor in powers of $1/c$, keeping the matter density and peculiar velocity as exact fundamental variables. We assume the Poisson gauge, including scalar and tensor modes up to $1/c^4$ order and vector modes up to $1/c^5$ terms. Through a redefinition of the scalar potentials as a resummation of the metric contributions at different orders, we obtain a complete set of nonlinear equations, providing a unified framework to study structure formation from small to superhorizon scales, from the nonlinear Newtonian to the linear relativistic regime. We explicitly show the validity of our scheme in the two limits: at leading order we recover the fully nonlinear equations of Newtonian cosmology; when linearized, our equations become those for scalar and vector modes of first-order relativistic perturbation theory in the Poisson gauge. Tensor modes are nondynamical at the $1/c^4$ order we consider: they are purely nonlinear and describe a distortion of the spatial slices determined at this order by a constraint, quadratic in the scalar and vector variables. The main results of our analysis are as follows: (a) at leading order a purely Newtonian nonlinear energy current sources a frame-dragging gravitomagnetic vector potential, and (b) in the leading-order Newtonian regime and in the linear relativistic regime the two scalar metric potentials are the same, while the nonlinearity of general relativity makes them different.

We show that a general late-time interaction between cold dark matter and vacuum energy is favoured by current cosmological datasets. We characterize the strength of the coupling by a dimensionless parameter $q_V$ that is free to take different values in four redshift bins from the primordial epoch up to today. This interacting scenario is in agreement with measurements of cosmic microwave background temperature anisotropies from the Planck satellite, supernovae Ia from Union 2.1 and redshift space distortions from a number of surveys, as well as with combinations of these different datasets. We show that a non-zero interaction is very likely at late times. We then focus on the case $q_V\not=0$ in a single low-redshift bin, obtaining a nested one parameter extension of the standard $\Lambda$CDM model. We study the Bayesian evidence, with respect to $\Lambda$CDM, of this late-time interaction model, finding moderate evidence for an interaction starting at $z=0.9$, dependent upon the prior range chosen for the interaction strength parameter $q_V$. For this case the null interaction ($q_V=0$, i.e.$\Lambda$CDM) is excluded at 99% c.l..

We show how the non-linearity of general relativity generates a characteristic non-Gaussian signal in cosmological large-scale structure that we calculate at all perturbative orders in a large scale limit. Newtonian gravity and general relativity provide complementary theoretical frameworks for modelling large-scale structure in $\Lambda$CDM cosmology; a relativistic approach is essential to determine initial conditions which can then be used in Newtonian simulations studying the non-linear evolution of the matter density. Most inflationary models in the very early universe predict an almost Gaussian distribution for the primordial metric perturbation, $\zeta$. However, we argue that it is the Ricci curvature of comoving-orthogonal spatial hypersurfaces, $R$, that drives structure formation at large scales. We show how the non-linear relation between the spatial curvature, $R$, and the metric perturbation, $\zeta$, translates into a specific non-Gaussian contribution to the initial comoving matter density that we calculate for the simple case of an initially Gaussian $\zeta$. Our analysis shows the non-linear signature of Einstein's gravity in large-scale structure.

A fast transition between a standard matter-like era and a late $\Lambda$CDM-like epoch generated by a single Unified Dark Matter component can explain the observed acceleration of the Universe. UDM models with a fast transition should be clearly distinguishable from $\Lambda$CDM (and alternatives) through observations. Here we focus on a particularly simple model and analyse its viability by studying features of the background model and properties of the adiabatic UDM perturbations.

The goal of the present paper is to set initial conditions for structure formation at non-linear order, consistent with general relativity, while also allowing for primordial non-Gaussianity. We use the non-linear continuity and Raychauduri equations, which together with the non-linear energy constraint determine the evolution of the matter density fluctuation in general relativity. We solve this equations at first and second order in a perturbative expansion, recovering and extending previous results derived in the matter-dominated limit and in the Newtonian regime. We present a second-order solution for the comoving density contrast in a $\Lambda$CDM universe, identifying non-linear contributions coming from the Newtonian growing mode, primordial non-Gaussianity and intrinsic non-Gaussianity, due to the essential non-linearity of the relativistic constraint equations. We discuss the application of these results to initial conditions in N-body simulations, showing that relativistic corrections mimic a non-zero non-linear parameter $f_{\rm NL}$.

We present the first calculation of an intrinsically relativistic quantity in fully non-linear cosmolog- ical large-scale structure studies. Traditionally, non-linear structure formation in standard {\Lambda}CDM cosmology is studied using N-body simulations, based on Newtonian gravitational dynamics on an expanding background. When one derives the Newtonian regime in a way that is a consistent ap- proximation to the Einstein equations, a gravito-magnetic vector potential - giving rise to frame dragging - is present in the metric in addition to the usual Newtonian scalar potential. At leading order, this vector potential does not affect the matter dynamics, thus it can be computed from Newtonian N-body simulations. We explain how we compute the vector potential from simulations in {\Lambda}CDM and examine its magnitude relative to the scalar potential. We also discuss some possible observable effects.

A fast transition between a standard matter-like era and a late $\Lambda$ cold dark matter ($\Lambda$CDM)-like epoch (or more in general, a CDM+DE era), generated by a single Unified Dark Matter component, can provide a new interesting paradigm in the context of general relativity, alternative to $\Lambda$CDM itself or other forms of DE or modified gravity theories invoked to explain the observed acceleration of the Universe. UDM models with a fast transition have interesting features, leading to measurable predictions, thus they should be clearly distinguishable from $\Lambda$CDM (and alternatives) through observations. Here we look at different ways of prescribing phenomenological UDM models with fast transition, then focusing on a particularly simple model. We analyse the viability of this model by studying features of the background model and properties of the adiabatic UDM perturbations, which depend on the effective speed of sound and the functional form of the Jeans scale. As a result, theoretical constraints on the parameters of the models are found that allow for a behaviour compatible with observations.

Local non-Gaussianity, parametrized by $f_{\rm NL}$, introduces a scale-dependent bias that is strongest at large scales, precisely where General Relativistic (GR) effects also become significant. With future data, it should be possible to constrain $f_{\rm NL} = {\cal O}(1)$ with high redshift surveys. GR corrections to the power spectrum and ambiguities in the gauge used to define bias introduce effects similar to $f_{\rm NL}= {\cal O}(1)$, so it is essential to disentangle these effects. For the first time in studies of primordial non-Gaussianity, we include the consistent GR calculation of galaxy power spectra, highlighting the importance of a proper definition of bias. We present observable power spectra with and without GR corrections, showing that an incorrect definition of bias can mimic non-Gaussianity. However, these effects can be distinguished by their different redshift and scale dependence, so as to extract the true primordial non-Gaussianity.

Motivated by the dawn of precision cosmology and the wealth of forthcoming high precision and volume galaxy surveys, in this paper we study the effects of inhomogeneities on light propagation in a flat \Lambda CDM background. To this end we use exact solutions of Einstein's equations (Meures & Bruni 2011) where, starting from small fluctuations, inhomogeneities arise from a standard growing mode and become non-linear. While the matter distribution in these models is necessarily idealised, there is still enough freedom to assume an arbitrary initial density profile along the line of sight. We can therefore model over-densities and voids of various sizes and distributions, e.g. single harmonic sinusoidal modes, coupled modes, and more general distributions in a \Lambda CDM background. Our models allow for an exact treatment of the light propagation problem, so that the results are unaffected by approximations and unambiguous. Along lines of sight with density inhomogeneities which average out on scales less than the Hubble radius, we find the distance redshift relation to diverge negligibly from the Friedmann-Lemaitre-Robertson-Walker (FLRW) result. On the contrary, if we observe along lines of sight which do not have the same average density as the background, we find large deviations from the FLRW distance redshift relation. Hence, a possibly large systematic might be introduced into the analysis of cosmological observations, e.g. supernovae, if we observe along lines of sight which are typically more or less dense than the average density of the Universe. In turn, this could lead to wrong parameter estimation: even if the Cosmological Principle is valid, the identification of the true FLRW background in an inhomogeneous universe maybe more difficult than usually assumed.

At a time when galaxy surveys and other observations are reaching unprecedented sky coverage and precision it seems timely to investigate the effects of general relativistic nonlinear dynamics on the growth of structures and on observations. Analytic inhomogeneous cosmological models are an indispensable way of investigating and understanding these effects in a simplified context. In this paper, we develop exact inhomogeneous solutions of general relativity with pressureless matter (dust, describing cold dark matter) and cosmological constant $\Lambda$, which can be used to model an arbitrary initial matter distribution along one line of sight. In particular, we consider the second class Szekeres models with $\Lambda$ and split their dynamics into a flat $\Lambda$CDM background and exact nonlinear inhomogeneities, obtaining several new results. One single metric function $Z$ describes the deviation from the background. We show that $F$, the time dependent part of $Z$, satisfies the familiar linear differential equation for $\delta$, the first-order density perturbation of dust, with the usual growing and decaying modes. In the limit of small perturbations, $\delta \approx F$ as expected, and the growth of inhomogeneities links up exactly with standard perturbation theory. We provide analytic expressions for the exact nonlinear $\delta$ and the growth factor in our models. For the case of over-densities, we find that, depending on the initial conditions, the growing mode may or may not lead to a pancake singularity, analogous to a Zel'dovich pancake. This is in contrast with the $\Lambda=0$ pure Einstein-de-Sitter background where, at any given point in comoving (Lagrangian) coordinates pancakes will always occur.

We investigate the general properties of Unified Dark Matter (UDM) scalar field models with Lagrangians with a non-canonical kinetic term, looking specifically for models that can produce a fast transition between an early Einstein-de Sitter CDM-like era and a later Dark Energy like phase, similarly to the barotropic fluid UDM models in JCAP1001(2010)014. However, while the background evolution can be very similar in the two cases, the perturbations are naturally adiabatic in fluid models, while in the scalar field case they are necessarily non-adiabatic. The new approach to building UDM Lagrangians proposed here allows to escape the common problem of the fine-tuning of the parameters which plague many UDM models. We analyse the properties of perturbations in our model, focusing on the the evolution of the effective speed of sound and that of the Jeans length. With this insight, we can set theoretical constraints on the parameters of the model, predicting sufficient conditions for the model to be viable. An interesting feature of our models is that what can be interpreted as w_{DE} can be <-1 without violating the null energy conditions.

We investigate the general properties of Unified Dark Matter (UDM) fluid models where the pressure and the energy density are linked by a barotropic equation of state (EoS) $p = p(\rho)$ and the perturbations are adiabatic. The EoS is assumed to admit a future attractor that acts as an effective cosmological constant, while asymptotically in the past the pressure is negligible. UDM models of the dark sector are appealing because they evade the so-called "coincidence problem" and "predict" what can be interpreted as $w_{\rm DE} \approx -1$, but in general suffer the effects of a non-negligible Jeans scale that wreak havoc in the evolution of perturbations, causing a large Integrated Sachs-Wolfe effect and/or changing structure formation at small scales. Typically, observational constraints are violated, unless the parameters of the UDM model are tuned to make it indistinguishable from $\Lambda$CDM. Here we show how this problem can be avoided, studying in detail the functional form of the Jeans scale in adiabatic UDM perturbations and introducing a class of models with a fast transition between an early Einstein-de Sitter CDM-like era and a later $\Lambda$CDM-like phase. If the transition is fast enough, these models may exhibit satisfactory structure formation and CMB fluctuations. To consider a concrete case, we introduce a toy UDM model and show that it can predict CMB and matter power spectra that are in agreement with observations for a wide range of parameter values.

Bouncing cosmologies are often proposed as alternatives to standard inflation for the explanation of the homogeneity and flatness of the universe. In such scenarios, the present cosmological expansion is preceded by a contraction phase. However, during the contraction, in general the anisotropy of the universe grows and eventually leads to a chaotic mixmaster behavior. This would either be hard to reconcile with observations or even lead to a singularity instead of the bounce. In order to preserve a smooth and isotropic bounce, the source for the contraction must have a super-stiff equation of state with $P/\rho=w>1$. In this letter we propose a new mechanism to solve the anisotropy problem for any low-energy value of $w$ by arguing that high energy physics leads to a modification of the equation of state, with the introduction of non-linear terms. In such a scenario, the anisotropy is strongly suppressed during the high energy phase, allowing for a graceful isotropic bounce, even when the low-energy value of $w$ is smaller than unity.

We study a set of universe models where the dark sector is described by a perfect fluid with an affine equation of state $P=P_0+\alpha \rho$, focusing specifically on cosmological perturbations in a flat universe. We perform a Monte Carlo Markov Chain analysis spanning the full parameter space of the model using the WMAP 5 years data and the SDSS LRG4 survey. The affine fluid can either play the role of a unified dark matter (UDM), accounting for both dark matter and a cosmological constant, or work alongside cold dark matter (CDM), as a form of dark energy. A key ingredient is the sound speed, that depends on the nature of the fluid and that, for any given background model, adds a degree of freedom to the perturbations: in the barotropic case the square of the sound speed is simply equal to the affine parameter $\alpha$; if entropic perturbations are present the effective sound speed has to be specified as an additional parameter. In addition to the barotropic case, we consider the two limiting cases of effective sound speed equal to 0 or 1. For $\alpha=c_s^2=0$ our UDM model is equivalent to the standard $\Lambda$CDM with adiabatic perturbations. Apart of a trivial subcase, all models considered satisfy the data constraints, with quite standard values for the usual cosmological parameters. In general our analysis confirms that cosmological datasets require both a collisionless massive and cold component to form the potential wells that lead to structure formation, and an effective cosmological constant that drives the late accelerated expansion.

We explore the dynamics of cosmological models with two coupled dark components with energy densities $\rho_A$ and $\rho_B$. We assume that the coupling is of the form $Q=Hq(\rho_A,\rho_B)$, so that the dynamics of the two components turns out to be scale independent, i.e. does not depend explicitly on the Hubble scalar $H$. With this assumption, we focus on the general linear coupling $q=q_o+q_A\rho_A+q_B\rho_B$, which may be seen as arising from any $q(\rho_A,\rho_B)$ at late time and leads in general to an effective cosmological constant. In the second part of the paper we consider observational constraints on the form of the coupling from SN Ia data, assuming that one of the components is cold dark matter. We find that the constant part of the coupling function is unconstrained by SN Ia data and, among typical linear coupling functions, the one proportional to the dark energy density $\rho_{A}$ is preferred in the strong coupling regime, $|q_{A}|>1$. While phantom models favor a positive coupling function, in non-phantom models, not only a negative coupling function is allowed, but the uncoupled sub-case falls at the border of the likelihood.

Loop Quantum Cosmology strongly modifies the high-energy dynamics of Friedman-Robertson-Walker models and removes the big-bang singularity. We investigate how LQC corrections affect the stability properties of the Einstein static universe. In General Relativity, the Einstein static model with positive cosmological constant Lambda is unstable to homogeneous perturbations. We show that LQC modifications can lead to a centre of stability for a large enough positive value of Lambda.

We consider the hypothesis that dark energy and dark matter are the two faces of a single dark component, a unified dark matter (UDM) that we assume can be modeled by the affine equation of state (EoS) $P= p_0 +\alpha \rho$, resulting in an {\it effective cosmological constant} $\rho_\Lambda=-p_0/(1+\alpha)$. The affine EoS arises from the simple assumption that the speed of sound is constant; it may be seen as an approximation to an unknown barotropic EoS $P=P(\rho)$, and may as well represent the tracking solution for the dynamics of a scalar field with appropriate potential. Furthermore, in principle the affine EoS allows the UDM to be phantom. We constrain the parameters of the model, $\alpha$ and $\Omega_\Lambda$, using data from a suite of different cosmological observations, and perform a comparison with the standard $\Lambda$CDM model, containing both cold dark matter and a cosmological constant. First considering a flat cosmology, we find that the UDM model with affine EoS fits the joint observations very well, better than $\Lambda$CDM, with best fit values $\alpha=0.01 \pm 0.02$ and $\Omega_\Lambda=0.70 \pm 0.04$ (95% confidence intervals). The standard model (best fit $\Omega_\Lambda=0.71\pm 0.04$), having one less parameter, is preferred by a Bayesian model comparison. However, the affine EoS is at least as good as the standard model if a flat curvature is not assumed as a prior for $\Lambda$CDM. For the latter, the best fit values are $\Omega_K=-0.02^{+0.01}_{-0.02} $ and $\Omega_\Lambda=0.71 \pm 0.04$, i.e. a closed model is preferred. A phantom UDM with affine EoS is ruled out well beyond $3\sigma$.

We explore the possibility that a scalar field with appropriate Lagrangian can mimic a perfect fluid with an affine barotropic equation of state. The latter can be thought of as a generic cosmological dark component evolving as an effective cosmological constant plus a generalized dark matter. As such, it can be used as a simple, phenomenological model for either dark energy or unified dark matter. Furthermore, it can approximate (up to first order in the energy density) any barotropic dark fluid with arbitrary equation of state. We find that two kinds of Lagrangian for the scalar field can reproduce the desired behaviour: a quintessence-like with a hyperbolic potential, or a purely kinetic k-essence one. We discuss the behaviour of these two classes of models from the point of view of the cosmological background, and we give some hints on their possible clustering properties.

In general relativity, for fluids with a linear equation of state (EoS) or scalar fields, the high isotropy of the universe requires special initial conditions, and singularities are anisotropic in general. In the brane world scenario anisotropy at the singularity is suppressed by an effective quadratic equation of state. There is no reason why the effective EoS of matter should be linear at the highest energies, and a non-linear EoS may describe dark energy or unified dark matter (Paper I, astro-ph/0512224). In view of this, here we study the effects of a quadratic EoS in homogenous and inhomogeneous cosmological models in general relativity, in order to understand if in this context the quadratic EoS can isotropize the universe at early times. With respect to Paper I, here we use the simplified EoS P=alpha rho + rho^2/rho_c, which still allows for an effective cosmological constant and phantom behavior, and is general enough to analyze the dynamics at high energies. We first study anisotropic Bianchi I and V models, focusing on singularities. Using dynamical systems methods, we find the fixed points of the system and study their stability. We find that models with standard non-phantom behavior are in general asymptotic in the past to an isotropic fixed point IS, i.e. in these models even an arbitrarily large anisotropy is suppressed in the past: the singularity is matter dominated. Using covariant and gauge invariant variables, we then study linear perturbations about the homogenous and isotropic spatially flat models with a quadratic EoS. We find that, in the large scale limit, all perturbations decay asymptotically in the past, indicating that the isotropic fixed point IS is the general asymptotic past attractor for non phantom inhomogeneous models with a quadratic EoS. (Abridged)

We investigate the non-linear coupling between radial and non-radial oscillations of static spherically symmetric neutron stars as a possible mechanism for the generation of gravitational waves that may lead to observable signatures. In this paper we concentrate on the axial sector of the non-radial perturbations. By using a multi-parameter perturbative framework we introduce a complete description of the non-linear coupling between radial and axial non-radial oscillations; we study the gauge invariant character of the associated perturbative variables and develop a computational scheme to evolve the non-linear coupling perturbations in the time domain. We present results of simulations corresponding to different physical situations and discuss the dynamical behaviour of this non-linear coupling. Of particular interest is the occurrence of signal amplifications in the form of resonance phenomena when a frequency associated with the radial pulsations is close to a frequency associated with one of the axial w-modes of the star. Finally, we mention possible extensions of this work and improvements towards more astrophysically motivated scenarios.

We investigate the general relativistic dynamics of Robertson-Walker models with a non-linear equation of state (EoS), focusing on the quadratic case P = P_0 + \alpha \rho + \beta \rho^2. This may be taken to represent the Taylor expansion of any arbitrary barotropic EoS, P(\rho). With the right combination of P_0, \alpha and \beta, it serves as a simple phenomenological model for dark energy, or even unified dark matter. Indeed we show that this simple model for the EoS can produce a large variety of qualitatively different dynamical behaviors that we classify using dynamical systems theory. An almost universal feature is that accelerated expansion phases are mostly natural for these non-linear EoS's. These are often asymptotically de Sitter thanks to the appearance of an effective cosmological constant. Other interesting possibilities that arise from the quadratic EoS are closed models that can oscillate with no singularity, models that bounce between infinite contraction/expansion and models which evolve from a phantom phase, asymptotically approaching a de Sitter phase instead of evolving to a "Big Rip". In a second paper we investigate the effects of the quadratic EoS in inhomogeneous and anisotropic models, focusing in particular on singularities.

We discuss an alternative approach to quintessence modifying the usual equation of state of the cosmological fluid in order to see if going further than the approximation of perfect fluid allows to better reproduce the available data. We consider a cosmological model comprising only two fluids, namely baryons (modelled as dust) and dark matter with a Van der Waals equation of state. First, the general features of the model are presented and then the evolution of the energy density, the Hubble parameter and the scale factor are determined showing that it is possible to obtain accelerated expansion choosing suitably the model parameters. We use the estimated age of the universe and the data on the dimensionless coordinate distances to Type Ia supernovae and distant radio galaxies to see whether Van der Waals quintessence is viable to explain dark energy and to constrain its parameters.

It has recently been suggested that an isotropic singularity may be a generic feature of brane cosmologies, even in the inhomogeneous case. Using the covariant and gauge-invariant approach we present a detailed analysis of linear perturbations of the isotropic model ${\cal F}_b$ which is a past attractor in the phase space of homogeneous Bianchi models on the brane. We find that for matter with an equation of state parameter $\gamma > 1$, the dimensionless variables representing generic anisotropic and inhomogeneous perturbations decay as $t\to 0$, showing that the model ${\cal F}_b$ is asymptotically stable in the past. We conclude that brane universes are born with isotropy naturally built-in, contrary to standard cosmology. The observed large-scale homogeneity and isotropy of the universe can therefore be explained as a consequence of the initial conditions if the brane-world paradigm represents a description of the very early universe.

Linear perturbation theory is appropriate to describe small oscillations of stars, while a mild non-linearity is still tractable perturbatively but requires to consider mode coupling. It is natural to start to look at this problem by considering the coupling between linear radial and non-radial modes. Radial pulsations of a spherical compact objects do not per se emit gravitational waves but, if the coupling is efficient in driving and possibly amplifying the non-radial modes, gravitational radiation could then be produced to a significant level. In this paper we develop the relativistic formalism to study the coupling of radial and non-radial first order perturbations of a compact spherical star. From a mathematical point of view, it is convenient to treat the two sets of perturbations as separately parametrized, using a 2-parameter perturbative expansion of the metric, the energy-momentum tensor and Einstein equations in which $\lambda$ is associated with the radial modes, $\epsilon$ with the non-radial perturbations, and the $\lambda\epsilon$ terms describe the coupling. This approach provides a well-defined framework to consider the gauge dependence of perturbations, allowing us to use $\epsilon$ order gauge-invariant non-radial variables on the static background and to define new second order $\lambda\epsilon$ gauge-invariant variables describing the non-linear coupling. We present the evolution and constraint equations for our variables outlining the setup for numerical computations, and briefly discuss the surface boundary conditions in terms of the second order $\lambda\epsilon$ Lagrangian pressure perturbation.

We compare three different models of rotating neutron star spacetimes: the Hartle-Thorne (HT) slow-rotation approximation at second order in rotation, the exact analytic vacuum solution of Manko et al. and a numerical solution of the full Einstein equations. We integrate the HT structure equations for five representative equations of state. Then we match the HT models to numerical solutions of the Einstein equations, imposing that the mass and angular momentum of the models be the same. We estimate the limits of validity of the HT expansion computing relative errors in the spacetime's quadrupole moment Q and in the ISCO radii. We find that ISCO radii computed in the HT approximation are accurate to better than 1%, even for the fastest observed ms pulsar. At the same rotational rates the accuracy on Q is of order 20%. In the second part of the paper we focus on the exterior vacuum spacetimes. We introduce a physically motivated `quasi-Kinnersley' Newman-Penrose frame. In this frame we evaluate the speciality index S, a coordinate-independent quantity measuring the deviation of each model from Petrov Type D. On the equatorial plane this deviation is smaller than 5%, even for the fastest rotating models. Our main conclusion is that the HT approximation is very reliable for most astrophysical applications.

In this talk I will briefly outline work in progress in two different contexts in astrophysical relativity, i.e. the study of rotating star spacetimes and the problem of reliably extracting gravitational wave templates in numerical relativity. In both cases the use of Weyl scalars and curvature invariants helps to clarify important issues.

This is a progress report on our study of the coupling of first-order radial and non-radial relativistic perturbations of a static spherical star. Our goal is to investigate the effects of this coupling on the gravitational wave signal of neutron stars. In particular, we are looking for the existence of resonances and parametric amplifications, changes in the damping time of non-radial oscillations, etc. To that end, we have developed a formalism that introduces gauge invariant quantities to describe the coupling. Their equations have the same structure as the equations for first-order non-radial perturbations plus some source terms, which makes them very appealing for time domain studies.

We introduce N-parameter perturbation theory as a new tool for the study of non-linear relativistic phenomena. The main ingredient in this formulation is the use of the Baker-Campbell-Hausdorff formula. The associated machinery allows us to prove the main results concerning the consistency of the scheme to any perturbative order. Gauge transformations and conditions for gauge invariance at any required order can then be derived from a generating exponential formula via a simple Taylor expansion. We outline the relation between our novel formulation and previous developments.

We present the complete set of propagation and constraint equations for the kinematic and non-local first order quantities which describe general linear inhomogeneous and anisotropic perturbations of a flat FRW braneworld with vanishing cosmological constant and decompose them in the standard way into their scalar, vector and tensor contributions. A detailed analysis of the perturbation dynamics is performed using dimensionless variables that are specially tailored for the different regimes of interest; namely, the low energy GR regime, the high energy regime and the dark energy regime. Tables are presented for the evolution of all the physical quantities, making it easy to do a detailed comparison of the past asymptotic behaviour of the perturbations of these models. We find results that exactly match those obtained in the analysis of the spatially inhomogeneous $G_{2}$ braneworld cosmologies presented recently; i.e., that isotropization towards the ${\cal F}_b$ model occurs for $\gamma > 4/3$.

In this paper we consider the Long-Wavelength Approximation Scheme (LWAS) in the framework of the covariant fluid approach to general relativistic dynamics, specializing to the particular case of irrotational dust matter. We discuss the dynamics of these models during the approach to any spacelike singularity where a BKL-type evolution is expected, studying the validity of this approximation scheme and the role of the magnetic part of the Weyl tensor, $H_{ab}$. Our analytic results confirm a previous numerical analysis: it is $H_{ab}$ that destroys the pure Kasner-like approach to the singularity and eventually produces the bounce to another Kasner phase. Expanding regions evolve as separate universes where inhomogeneities and anisotropies die away.

An implicit fundamental assumption in relativistic perturbation theory is that there exists a parametric family of spacetimes that can be Taylor expanded around a background. The choice of the latter is crucial to obtain a manageable theory, so that it is sometime convenient to construct a perturbative formalism based on two (or more) parameters. The study of perturbations of rotating stars is a good example: in this case one can treat the stationary axisymmetric star using a slow rotation approximation (expansion in the angular velocity Omega), so that the background is spherical. Generic perturbations of the rotating star (say parametrized by lambda) are then built on top of the axisymmetric perturbations in Omega. Clearly, any interesting physics requires non-linear perturbations, as at least terms lambda Omega need to be considered. In this paper we analyse the gauge dependence of non-linear perturbations depending on two parameters, derive explicit higher order gauge transformation rules, and define gauge invariance. The formalism is completely general and can be used in different applications of general relativity or any other spacetime theory.

We use the Zel'dovich approximation to analyse the amplification of magnetic fields in gravitational collapse of cold dark matter during the mildly nonlinear regime, and identify two key features. First, the anisotropy of collapse effectively eliminates one of the magnetic components, confining the field in the pancake plane. Second, in agreement with recent numerical simulations, we find that the shear anisotropy can amplify the magnetic field well beyond the isotropic case. Our results suggest that the magnetic strengths observed in spiral and disk galaxies today might have originated from seeds considerably weaker than previous estimates.

Recent studies of homogeneous anisotropic universe models in the brane world scenario show that the cosmological singularity in this context is isotropic. It has therefore been suggested that this may be a generic feature of singularities on the brane, even in the inhomogeneous case. Using a perturbative approach, we show that this is not the case. As in the GR case, the presence of decaying modes in the perturbations signal the instability (in the past) of the isotropic singularity. The brane universe is therefore not born with isotropy built in: as in standard cosmology, the observed large-scale isotropy and homogeneity remains to be explained.

We summarize recent results concerning the evolution of second order perturbations in flat dust irrotational FLRW models with $\Lambda\ne 0$. We show that asymptotically these perturbations tend to constants in time, in agreement with the cosmic no-hair conjecture. We solve numerically the second order scalar perturbation equation, and very briefly discuss its all time behaviour and some possible implications for the structure formation.

We investigate how braneworld gravity affects gravitational collapse and black hole formation by studying Oppenheimer-Snyder-like collapse on a Randall-Sundrum type brane. Without making any assumptions about the bulk, we prove a no-go theorem: the exterior spacetime on the brane cannot be static, which is in stark contrast with general relativity. We also consider the role of Kaluza-Klein energy density in collapse, using a toy model.

We study the asymptotic stability of de Sitter spacetime with respect to non-linear perturbations, by considering second order perturbations of a flat Robertson-Walker universe with dust and a positive cosmological constant. Using the synchronous comoving gauge we find that, as in the case of linear perturbations, the non-linear perturbations also tend to constants, asymptotically in time. Analysing curvature and other spacetime invariants we show, however, that these quantities asymptotically tend to their de Sitter values, thus demonstrating that the geometry is indeed locally asymptotically de Sitter, despite the fact that matter inhomogeneities tend to constants in time. Our results support the inflationary picture of frozen amplitude matter perturbations that are stretched outside the horizon, and demonstrate the validity of the cosmic no-hair conjecture in the nonlinear inhomogeneous settings considered here.

We discuss the issue of observables in general-relativistic perturbation theory, adopting the view that any observable in general relativity is represented by a scalar field on spacetime. In the context of perturbation theory, an observable is therefore a scalar field on the perturbed spacetime, and as such is gauge invariant in an exact sense (to all orders), as one would expect. However, perturbations are usually represented by fields on the background spacetime, and expanded at different orders into contributions that may or may not be gauge independent. We show that perturbations of scalar quantities are observable if they are first order gauge-invariant, even if they are gauge dependent at higher order. Gauge invariance to first order plays therefore an important conceptual role in the theory, for it selects the perturbations with direct physical meaning from those having only a mathematical status. The so-called ``gauge problem'', and the relationship between measured fluctuations and gauge dependent perturbations that are computed in the theory are also clarified.

We consider the evolution of relativistic perturbations in the Einstein-de Sitter cosmological model, including second-order effects. The perturbations are considered in two different settings: the widely used synchronous gauge and the Poisson (generalized longitudinal) one. Since, in general, perturbations are gauge dependent, we start by considering gauge transformations at second order. Next, we give the evolution of perturbations in the synchronous gauge, taking into account both scalar and tensor modes in the initial conditions. Using the second-order gauge transformation previously defined, we are then able to transform these perturbations to the Poisson gauge. The most important feature of second-order perturbation theory is mode-mixing, which here also means, for instance, that primordial density perturbations act as a source for gravitational waves, while primordial gravitational waves give rise to second-order density fluctuations. Possible applications of our formalism range from the study of the evolution of perturbations in the mildly non-linear regime to the analysis of secondary anisotropies of the Cosmic Microwave Background.

Diffeomorphism freedom induces a gauge dependence in the theory of spacetime perturbations. We derive a compact formula for gauge transformations of perturbations of arbitrary order. To this end, we develop the theory of Taylor expansions for one-parameter families (not necessarily groups) of diffeomorphisms. First, we introduce the notion of knight diffeomorphism, that generalises the usual concept of flow, and prove a Taylor's formula for the action of a knight on a general tensor field. Then, we show that any one-parameter family of diffeomorphisms can be approximated by a family of suitable knights. Since in perturbation theory the gauge freedom is given by a one-parameter family of diffeomorphisms, the expansion of knights is used to derive our transformation formula. The problem of gauge dependence is a purely kinematical one, therefore our treatment is valid not only in general relativity, but in any spacetime theory.

We consider in detail the problem of gauge dependence that exists in relativistic perturbation theory, going beyond the linear approximation and treating second and higher order perturbations. We first derive some mathematical results concerning the Taylor expansion of tensor fields under the action of one-parameter families (not necessarily groups) of diffeomorphisms. Second, we define gauge invariance to an arbitrary order $n$. Finally, we give a generating formula for the gauge transformation to an arbitrary order and explicit rules to second and third order. This formalism can be used in any field of applied general relativity, such as cosmological and black hole perturbations, as well as in other spacetime theories. As a specific example, we consider here second order perturbations in cosmology, assuming a flat Robertson-Walker background, giving explicit second order transformations between the synchronous and the Poisson (generalized longitudinal) gauges.

We present the results of a study of the gauge dependence of spacetime perturbations. In particular, we consider gauge invariance in general, we give a generating formula for gauge transformations to an arbitrary order n, and explicit transformation rules at second order.

After a brief overview of the so-called silent models and their present status, we consider the subclass of Bianchi Type--I models with a magnetic field source. Due to the presence of the magnetic field, the initial singularity shows ``oscillatory'' features reminiscent of the Bianchi Type--IX case. The Bianchi Type--I models with a magnetic field are therefore a counterexample to the folklore that matter fields can be neglected in the vicinity of the singularity.

We present results on the non-linear dynamics of inhomogeneous cosmological models with irrotational dust and a positive cosmological constant, considering, in particular, a wide class with vanishing magnetic Weyl tensor. For those patches of the universe that do not recollapse to a singularity we find a unique attractor, representing a de Sitter vacuum phase. For the (re-)collapsing regions we find a family of (Kasner) attractors, so that generically these regions fall in spindle-like singularities. These results give substantial support to the idea that the universe can be very inhomogeneous on ultra-large, super-horizon scales, whit observers living in those (almost) isotropic regions that emerge from an inflationary phase.

We investigate the local non--linear dynamics of irrotational dust with vanishing magnetic part of the Weyl tensor, $H_{ab}$. Once coded in the initial conditions, this dynamical restriction is respected by the relativistic evolution equations. Thus, the outcome of the latter are {\it exact solutions} for special initial conditions with $H_{ab}=0$, but with no symmetries: they describe inhomogeneous triaxial dynamics generalizing that of a fluid element in a Tolman--Bondi, Kantowski--Sachs or Szekeres geometry. A subset of these solutions may be seen as (special) perturbations of Friedmann models, in the sense that there are trajectories in phase--space that pass arbitrarily close to the isotropic ones. We find that the final fate of ever--expanding configurations is a spherical void, locally corresponding to a Milne universe. For collapsing configurations we find a whole family of triaxial attractors, with vanishing local density parameter $\Omega$. These attractors locally correspond to Kasner vacuum solutions: there is a single physical configuration collapsing to a degenerate {\it pancake}, while the generic configuration collapses to a triaxial {\it spindle} singularity. These {\it silent universe} models may provide a fair representation of the universe on super horizon scales. Moreover, one might conjecture that the non--local information carried by $H_{ab}$ becomes negligible during the late highly non--linear stages of collapse, so that the attractors we find may give all of the relevant expansion or collapse configurations of irrotational dust.

This is a somewhat extended version of the original July 93 report. It is proved that the cosmological density perturbation is associated with a peculiar velocity field. This allows a simple formulation of cosmological perturbation theory, which works entirely with quasi-Newtonian fluid flow equations. As an illustration, the large scale cosmic microwave background anisotropy (Sachs-Wolfe effect) is calculated without any reference to the metric perturbation. In addition, assuming the usual adiabatic initial condition on the density perturbation, we show that the dipole of the anisotropy measures our peculiar velocity relative to the average peculiar velocity within the last scattering surface of the microwave background, thus defining its frame of reference.

Flat and open universe models are considered, containing a mixture of cold matter (dust) and radiation interacting only through gravity, with the aim of studying their stability with respect to linear scalar perturbations. To this end the perturbed universe is considered as a dynamical system, described by coupled differential equations for a gauge\hs invariant perturbation variable and a relevant background variable. The phase\hs space analysis of this dynamical system shows that flat dust\hs radiation models are unstable, and open models structurally unstable, with respect to adiabatic perturbations. It is shown that there are perturbations which decay even if their wavelength at equidensity is larger than the corresponding Jeans scale. Metric and curvature perturbations are also briefly discussed. We believe that this analysis gives a clearer idea of the stability properties of realistic universe models than the standard one based on the Jeans scale, despite our simplifying assumptions.