We examine the implementation of Buchert's and Green & Wald's averaging formalisms in exact spherically symmetric and plane symmetric dust-filled cosmological models. We find that, given a cosmological space-time, Buchert's averaging scheme gives a faithful way of interpreting the large-scale expansion of space, and explicit terms that precisely quantify deviations from the behaviour expected from the Friedmann equations of homogeneous and isotropic cosmological models. The Green & Wald formalism, on the other hand, does not appear to yield any information about the large-scale properties of a given inhomogeneous space-time. Instead, this formalism is designed to calculate the back-reaction effects of short-wavelength fluctuations around a given "background" geometry. We find that the inferred expansion of space in this approach is entirely dependent on the choice of this background, which is not uniquely specified for any given inhomogeneous space-time, and that the "back-reaction" from small-scale structures vanishes in every case we study. This would appear to limit the applicability of Green & Wald's formalism to the study of large-scale expansion in the real Universe, which also has no pre-defined background. Further study is required to enhance the evaluation and comparison of these averaging formalisms, and determine whether the same difficulties exist, in less idealized space-time geometries.

We propose in this letter a relativistic coordinate independent interpretation for Milgrom's acceleration $a_{0}=1.2 \times 10^{-8} \hbox{cm/s}^{2}$ through a geometric constraint obtained from the product of the Kretschmann invariant scalar times the surface area of 2--spheres defined through suitable characteristic length scales for local and cosmic regimes, described by Schwarzschild and Friedman--Lema\^\i tre--Robertson--Walker (FLRW) geometries, respectively. By demanding consistency between these regimes we obtain an appealing expression for the empirical (so far unexplained) relation between the accelerations $a_0$ and $c H_0$. Imposing this covariant geometric criterion upon a FLRW model, yields a dynamical equation for the Hubble scalar whose solution matches, to a very high accuracy, the cosmic expansion rate of the $\Lambda$CDM concordance model fit for cosmic times close to the present epoch. We believe that this geometric interpretation of $a_0$ could provide relevant information for a deeper understanding of gravity

We examine the fully relativistic evolution of cosmic voids constituted by baryons and cold dark matter (CDM), represented by two non-comoving dust sources in a $\Lambda$CDM background. For this purpose, we consider numerical solutions of Einstein's field equations in a fluid-flow representation adapted to spherical symmetry and multiple components. We present a simple example that explores the frame-dependence of the local expansion and the Hubble flow for this mixture of two dusts, revealing that the relative velocity between the sources yields a significantly different evolution in comparison with that of the two sources in a common 4-velocity (which reduces to a Lemaitre-Tolman-Bondi model). In particular, significant modifications arise for the density contrast depth and void size, as well as in the amplitude of the surrounding over-densities. We show that an adequate model of a frame-dependent evolution that incorporates initial conditions from peculiar velocities and large-scale density contrast observations may contribute to understand the discrepancy between the local value of $H_0$ and that inferred from the CMB.

Understanding the observations of dynamical tracers and the trajectories of lensed photons at galactic scales within the context of General Relativity (GR), requires the introduction of a hypothetical dark matter dominant component. The onset of these gravitational anomalies, where the Schwarzschild solution no longer describes observations, closely corresponds to regions where accelerations drop below the characteristic $a_{0}$ acceleration of MOND, which occur at a well established mass-dependent radial distance, $R_{c}\propto (GM/a_{0})^{1/2}$. At cosmological scales, inferred dynamics are also inconsistent with GR and the observed distribution of mass. The current accelerated expansion rate requires the introduction of a hypothetical dark energy dominant component. We here show that for a Schwarzschild metric at galactic scales, the scalar curvature, K, multiplied by $(r^{4}/M)$ at the critical MOND transition radius, $r=R_{c}$, has an invariant value of $\kappa_{B}=K(r^{4}/M)=28Ga_{0}/c^{4}$. Further, assuming this condition holds for $r>R_{c}$, is consistent with the full spacetime which under GR corresponds to a dominant isothermal dark matter halo, to within observational precision at galactic level. For a FLRW metric, this same constant bounding curvature condition yields for a spatially flat spacetime a cosmic expansion history which agrees with the $\Lambda$CDM empirical fit for recent epochs, and which similarly tends asymptotically to a de Sitter solution. Thus, a simple covariant purely geometric condition identifies the low acceleration regime of observed gravitational anomalies, and can be used to guide the development of { extended} gravity theories at both galactic and cosmological scales.

We examine the gravitational collapse and black hole formation of multiple non--spherical configurations constructed from Szekeres dust models with positive spatial curvature that smoothly match to a Schwarzschild exterior. These configurations are made of an almost spherical central core region surrounded by a network of "pancake-like" overdensities and voids with spatial positions prescribed through standard initial conditions. We show that a full collapse into a focusing singularity, without shell crossings appearing before the formation of an apparent horizon, is not possible unless the full configuration becomes exactly or almost spherical. Seeking for black hole formation, we demand that shell crossings are covered by the apparent horizon. This requires very special fine-tuned initial conditions that impose very strong and unrealistic constraints on the total black hole mass and full collapse time. As a consequence, non-spherical non-rotating dust sources cannot furnish even minimally realistic toy models of black hole formation at astrophysical scales: demanding realistic collapse time scales yields huge unrealistic black hole masses, while simulations of typical astrophysical black hole masses collapse in unrealistically small times. We note, however, that the resulting time--mass constraint is compatible with early Universe models of primordial black hole formation, suitable in early dust-like environments. Finally, we argue that the shell crossings appearing when non-spherical dust structures collapse are an indicator that such structures do not form galactic mass black holes but virialise into stable stationary objects.

Recently we have studied in great detail a model of Hybrid Natural Inflation (HNI) by constructing two simple effective field theories. These two versions of the model allow inflationary energy scales as small as the electroweak scale in one of them or as large as the Grand Unification scale in the other, therefore covering the whole range of possible energy scales. In any case the inflationary sector of the model is of the form $V(\phi)=V_0 \left(1+a \cos(\phi/f)\right)$ where $0\leq a<1$ and the end of inflation is triggered by an independent waterfall field. One interesting characteristic of this model is that the slow-roll parameter $\epsilon(\phi)$ is a non-monotonic function of $\phi$ presenting a {\it maximum} close to the inflection point of the potential. Because the scalar spectrum $\mathcal{P}_s(k)$ of density fluctuations when written in terms of the potential is inversely proportional to $\epsilon(\phi)$ we find that $\mathcal{P}_s(k)$ presents a {\it minimum} at $\phi_{min}$. The origin of the HNI potential can be traced to a symmetry breaking phenomenon occurring at some energy scale $f$ which gives rise to a (massless) Goldstone boson. Non-perturbative physics at some temperature $T<f$ might occur which provides a potential (and a small mass) to the originally massless boson to become the inflaton (a pseudo-Nambu-Goldstone boson). Thus the inflaton energy scale $\Delta$ is bounded by the symmetry breaking scale, $\Delta\equiv V_H^{1/4} <f.$ To have such a well defined origin and hierarchy of scales in inflationary models is not common. We use this property of HNI to determine bounds for the inflationary energy scale $\Delta$ and for the tensor-to-scalar ratio $r$.

We derive a class of non-static inhomogeneous dust solutions in f(R) gravity described by the Lemaitre-Tolman-Bondi (LTB) metric. The field equations are fully integrated for all parameter subcases and compared with analogous subcases of LTB dust solutions of GR. Since the solutions do not admit regular symmetry centres, we have two possibilities: (i) a spherical dust cloud with angle deficit acting as the source of a vacuum Schwarzschild-like solution associated with a global monopole, or (ii) fully regular dust wormholes without angle deficit, whose rest frames are homeomorphic to the Schwarzschild-Kruskal manifold or to a 3d torus. The compatibility between the LTB metric and generic f(R) ansatzes furnishes an "inverse procedure" to generate LTB solutions whose sources are found from the f(R) geometry. While the resulting fluids may have an elusive physical interpretation, they can be used as exact non--perturbative toy models in theoretical and cosmological applications of f(R) theories.

We study the differences and equivalences between the non-perturbative description of the evolution of cosmic structure furnished by the Szekeres dust models (a non-spherical exact solution of Einstein's equations) and the dynamics of Cosmological Perturbation Theory (CPT) for dust sources in a $\Lambda$CDM background. We show how the dynamics of Szekeres models can be described by evolution equations given in terms of "exact fluctuations" that identically reduce (at all orders) to evolution equations of CPT in the comoving isochronous gauge. We explicitly show how Szekeres linearised exact fluctuations are specific (deterministic) realisations of standard linear perturbations of CPT given as random fields but, as opposed to the latter perturbations, they can be evolved exactly into the full non-linear regime. We prove two important results: (i) the conservation of the curvature perturbation (at all scales) also holds for the appropriate approximation of the exact Szekeres fluctuations in a $\Lambda$CDM background, and (ii) the different collapse morphologies of Szekeres models yields, at nonlinear order, different functional forms for the growth factor that follows from the study of redshift space distortions. The metric based potentials used in linear CPT are computed in terms of the parameters of the linearised Szekeres models, thus allowing us to relate our results to linear CPT results in other gauges. We believe that these results provide a solid starting stage to examine the role of non-perturbative General Relativity in current cosmological research.

We take a pragmatic, model independent approach to single field slow-roll canonical inflation by imposing conditions, not on the potential, but on the slow-roll parameter $\epsilon(\phi)$ and its derivatives $\epsilon^{\prime }(\phi)$ and $\epsilon^{\prime\prime }(\phi)$, thereby extracting general conditions on the tensor-to-scalar ratio $r$ and the running $n_{sk}$ at $\phi_{H}$ where the perturbations are produced, some $50$ $-$ $60$ $e$-folds before the end of inflation. We find quite generally that for models where $\epsilon(\phi)$ develops a maximum, a relatively large $r$ is most likely accompanied by a positive running while a negligible tensor-to-scalar ratio implies negative running. The definitive answer, however, is given in terms of the slow-roll parameter $\xi_2(\phi)$. To accommodate a large tensor-to-scalar ratio that meets the limiting values allowed by the Planck data, we study a non-monotonic $\epsilon(\phi)$ decreasing during most part of inflation. Since at $\phi_{H}$ the slow-roll parameter $\epsilon(\phi)$ is increasing, we thus require that $\epsilon(\phi)$ develops a maximum for $\phi > \phi_{H}$ after which $\epsilon(\phi)$ decrease to small values where most $e$-folds are produced. The end of inflation might occur trough a hybrid mechanism and a small field excursion $\Delta\phi_e\equiv |\phi_H-\phi_e |$ is obtained with a sufficiently thin profile for $\epsilon(\phi)$ which, however, should not conflict with the second slow-roll parameter $\eta(\phi)$. As a consequence of this analysis we find bounds for $\Delta \phi_e$, $r_H$ and for the scalar spectral index $n_{sH}$. Finally we provide examples where these considerations are explicitly realised.

We show that the full dynamical freedom of the well known Szekeres models allows for the description of elaborated 3--dimensional networks of cold dark matter structures (over--densities and/or density voids) undergoing "pancake" collapse. By reducing Einstein's field equations to a set of evolution equations, which themselves reduce in the linear limit to evolution equations for linear perturbations, we determine the dynamics of such structures, with the spatial comoving location of each structure uniquely specified by standard early Universe initial conditions. By means of a representative example we examine in detail the density contrast, the Hubble flow and peculiar velocities of structures that evolved, from linear initial data at the last scattering surface, to fully non--linear 10--20 Mpc. scale configurations today. To motivate further research, we provide a qualitative discussion on the connection of Szekeres models with linear perturbations and the pancake collapse of the Zeldovich approximation. This type of structure modelling provides a coarse grained -- but fully relativistic non--linear and non--perturbative -- description of evolving large scale cosmic structures before their virialisation, and as such it has an enormous potential for applications in cosmological research.

We undertake a hydrodynamical study of a mixture of tightly coupled primordial radiation, neutrinos, baryons, electrons and positrons, together with a gas of already decoupled dark matter WIMPS and an already existing "frozen" magnetic field in the infinite conductivity regime. Considering this cosmic fluid as the source of a homogeneous but anisotropic Bianchi I model, we describe its interaction with the magnetic field by means of suitable equations of state that are appropriate for the particle species of the mixture between the end of the leptonic era and the beginning of the radiation-dominated epoch. Fulfilment of observational bounds on the magnetic field intensity yields a "near FLRW" (but strictly non-perturbative) evolution of the geometric, kinematic and thermodynamical variables. This evolution is roughly comparable to the weak field approximation in linear perturbations on a spatially flat FLRW background of sources in which the frozen magnetic fields are coherent over very large supra-horizon scales. Our approach and results may provide interesting guidelines in potential situations in which non-perturbative methods are required to study the interaction between magnetic fields and the cosmic fluid.

We examine the spatial extrema (local maxima, minima and saddle points) of the covariant scalars (density, Hubble expansion, spatial curvature and eigenvalues of the shear and electric Weyl tensors) of the quasi-spherical Szekeres dust models. Sufficient conditions are obtained for the existence of distributions of multiple extrema in spatial comoving locations that can be prescribed through initial conditions. These distributions evolve without shell crossing singularities at least for ever expanding models (with or without cosmological constant) in the full evolution range where the models are valid. By considering the local maxima and minima of the density, our results allow for setting up elaborated networks of "pancake" shaped evolving cold dark matter over-densities and density voids whose spatial distribution and amplitudes can be controlled from initial data compatible with standard early Universe initial conditions. We believe that these results have an enormous range of potential application by providing a fully relativistic non-perturbative coarse grained modelling of cosmic structure at all scales.

We undertake a non-perturbative study of the evolution of the "gravitational entropy" proposed by Clifton, Ellis and Tavakol (CET) on local expanding cosmic CDM voids of $\sim 50-100$ Mpc size described as spherical under-dense regions with negative spatial curvature, whose dynamics is determined by Lemaitre-Tolman-Bondi (LTB) dust models asymptotic to three different types of FLRW background: $\Lambda$CDM, Einstein de Sitter and "open" FLRW with $\Lambda=0$ and negative spatial curvature. By assuming generic nearly spatially flat and linear initial conditions at the last scattering time, we examine analytically and numerically the CET entropy evolution into a fully non-linear regime in our present cosmic time and beyond. Both analytic and numerical analysis reveal that the late time CET entropy growth is determined by the amplitude of initial fluctuations of spatial curvature at the last scattering time. This entropy growth decays to zero in the late asymptotic time range for all voids, but at a faster rate in voids with $\Lambda$CDM and open FLRW backgrounds. However, only for voids in a $\Lambda$CDM background this suppression is sufficiently rapid for the CET entropy itself to reach a terminal equilibrium (or "saturation") value. The CET gravitational temperature vanishes asymptotically if $\Lambda=0$ and becomes asymptotically proportional to $\Lambda$ for voids in a $\Lambda$CDM background. In the linear regime of the LTB evolution our results coincide, qualitatively and quantitatively, with previous results based on linear perturbation theory.

We examine the relation between the dynamics of Lema\^{\i}tre-Tolman-Bondi (LTB) dust models (with and without $\Lambda$) and the dynamics of dust perturbations in two of the more familiar formalisms used in cosmology: the metric based Cosmological Perturbation Theory (CPT) and the Covariant Gauge Invariant (GIC) perturbations. For this purpose we recast the evolution of LTB models in terms of a covariant and gauge invariant formalism of local and non-local "exact fluctuations " on a Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) background defined by suitable averages of covariant scalars. We examine the properties of these fluctuations, which can be defined for a confined comoving domain or for an asymptotic domain extending to whole time slices. In particular, the non-local density fluctuation provides a covariant and precise definition for the notion of the "density contrast ". We show that in their linear regime these LTB exact fluctuations (local and non-local) are fully equivalent to the conventional cosmological perturbations in the synchronous-comoving gauge of CPT and to GIC perturbations. As an immediate consequence, we show the time-invariance of the spatial curvature perturbation in a simple form. The present work may provide important theoretical connections between the exact and perturbative (linear or no-linear) approach to the dynamics of dust sources in General Relativity.

We apply a recent proposal to define "gravitational entropy" to the expansion of cosmic voids within the framework of non-perturbative General Relativity. By considering CDM void configurations compatible with basic observational constraints, we show that this entropy grows from post-inflationary conditions towards a final asymptotic value in a late time fully non-linear regime described by the Lemaitre-Tolman-Bondi (LTB) dust models. A qualitatively analogous behavior occurs if we assume a positive cosmological constant consistent with a $\Lambda$-CDM background model. However, the $\Lambda$ term introduces a significant suppression of entropy growth with the terminal equilibrium value reached at a much faster rate.

We analyse Hybrid Natural Inflation in view of the recent results for the tensor index reported by BICEP2. We find that it predicts a large running of the scalar spectrum which is potentially detectable by large scale structure and $21\, \mathrm{cm} $ observations. The running of the running is also relatively large becoming close to $10^{-2}$. Along the way, we find general consistency relations at which observables are subject if the slow-roll approximation is imposed. Failure to satisfy these equations by the values obtained for the observables in surveys would be a failure of the slow-roll approximation itself.

We obtain covariant expressions that generalize the growing and decaying density modes of linear perturbation theory of dust sources by means of the exact density perturbation from the formalism of quasi--local scalars associated to weighed proper volume averages in LTB dust models. The relation between these density modes and theoretical properties of generic LTB models is thoroughly studied by looking at the evolution of the models through a dynamical system whose phase space is parametrized by variables directly related to the modes themselves. The conditions for absence of shell crossings, as well as sign conditions on the modes, become interrelated fluid flow preserved constraints that define phase space invariant subspaces. In the general case (both density modes being nonzero) the evolution of phase space trajectories exhibits the expected dominance of the decaying/growing in the early/late evolution times defined by past/future attractors characterized by asymptotic density inhomogeneity. In particular, the growing mode is also dominant for collapsing layers that terminate in a future attractor associated with a "Big Crunch" singularity, which is qualitatively different from the past attractor marking the "Big Bang". Suppression of the decaying mode modifies the early time evolution, with phase space trajectories emerging from an Einstein--de Sitter past attractor associated with homogeneous conditions. Suppression of the growing mode modifies the late time evolution as phase space trajectories terminate in future attractors associated with homogeneous states. General results are obtained relating the signs of the density modes and the type of asymptotic density profile (clump or void). A critical review is given of previous attempts in the literature to define these density modes for LTB models.

We introduce a weighed scalar average formalism ("q-average") for the study of the theoretical properties and the dynamics of spherically symmetric Lemaitre-Tolman-Bondi (LTB) dust models models. The "q-scalars" that emerge by applying the q-averages to the density, Hubble expansion and spatial curvature (which are common to FLRW models) are directly expressible in terms of curvature and kinematic invariants and identically satisfy FLRW evolution laws without the back-reaction terms that characterize Buchert's average. The local and non-local fluctuations and perturbations with respect to the q-average convey the effects of inhomogeneity through the ratio of curvature and kinematic invariants and the magnitude of radial gradients. All curvature and kinematic proper tensors that characterize the models are expressible as irreducible algebraic expansions on the metric and 4-velocity, whose coefficients are the q-scalars and their linear and quadratic local fluctuations. All invariant contractions of these tensors are quadratic fluctuations, whose q-averages are directly and exactly related to statistical correlation moments of the density and Hubble expansion scalar. We explore the application of this formalism to a definition of a gravitational entropy functional proposed by Hosoya et al (2004 Phys. Rev. Lett. 92 141302). We show that a positive entropy production follows from a negative correlation between fluctuations of the density and Hubble scalar, providing a brief outline on its fulfillment in various LTB models and regions. While the q-average formalism is specially suited for LTB and Szekeres models, it may provide a valuable theoretical insight on the properties of scalar averaging in inhomogeneous spacetimes in general.

Inhomogeneous cosmological models are able to fit cosmological observations without dark energy under the assumption that we live close to the "center" of a very large-scale under-dense region. Most studies fitting observations by means of inhomogeneities also assume spherical symmetry, and thus being at (or very near) the center may imply being located at a very special and unlikely observation point. We argue that such spherical voids should be treated only as a gross first approximation to configurations that follow from a suitable smoothing out of the non-spherical part of the inhomogeneities on angular scales. In this Letter we present a toy construction that supports the above statement. The construction uses parts of the Szekeres model, which is inhomogeneous and anisotropic thus it also addresses the limitations of spherical inhomogeneities. By using the thin-shell approximation (which means that the Israel-Darmois continuity conditions are not fulfilled between the shells) we construct a model of evolving cosmic structures, containing several elongated supercluster-like structures with underdense regions between them, which altogether provides a reasonable coarse-grained description of cosmic structures. While this configuration is not spherically symmetric, its proper volume average yields a spherical void profile of 250 Mpc that roughly agrees with observations. Also, by considering a non-spherical inhomogeneity, the definition of a "center" location becomes more nuanced, and thus the constraints placed by fitting observations on our position with respect to this location become less restrictive.

We consider spherically symmetric inhomogeneous dust models with a positive cosmological constant, $\Lambda$, given by the Lemaitre-Tolman-Bondi metric. These configurations provide a simple but useful generalization of the Lambda-CDM model describing cold dark matter (CDM) and a Lambda term, which seems to fit current cosmological observations. The dynamics of these models can be fully described by scalar evolution equations that can be given in the form of a proper dynamical system associated with a 4-dimensional phase space whose critical points and invariant subspaces are examined and classified. The phase space evolution of various configurations is studied in detail by means of two 2-dimensional subspaces: a projection into the invariant homogeneous subspace associated with Lambda-CDM solutions with FLRW metric, and a projection into a subspace generated by suitably defined fluctuations that convey the effects of inhomogeneity. We look at cases with perpetual expansion, bouncing and loitering behavior, as well as configurations with "mixed" kinematic patters, such as a collapsing region in an expanding background. In all cases, phase space trajectories emerge from and converge to stable past and future attractors in a qualitatively analogous way as in the case of the FLRW limit. However, we can identify in both projections of the phase space various qualitative features absent in the FLRW limit that can be useful in the construction of toy models of astrophysical and cosmological inhomogeneities.

We undertake a comprehensive and rigorous analytic study of the evolution of radial profiles of covariant scalars in regular Lemaitre-Tolman-Bondi dust models. We consider specifically the phenomenon of "profile inversions" in which an initial clump profile of density, spatial curvature or the expansion scalar, might evolve into a void profile (and vice versa). Previous work in the literature on models with density void profiles and/or allowing for density profile inversions is given full generalization, with some erroneous results corrected. We prove rigorously that if an evolution without shell crossings is assumed, then only the 'clump to void' inversion can occur in density profiles, and only in hyperbolic models or regions with negative spatial curvature. The profiles of spatial curvature follow similar patterns as those of the density, with 'clump to void' inversions only possible for hyperbolic models or regions. However, profiles of the expansion scalar are less restrictive, with profile inversions necessarily taking place in elliptic models. We also examine radial profiles in special LTB configurations: closed elliptic models, models with a simultaneous big bang singularity, as well as a locally collapsing elliptic region surrounded by an expanding hyperbolic background. The general analytic statements that we obtain allow for setting up the right initial conditions to construct fully regular LTB models with any specific qualitative requirements for the profiles of all scalars and their time evolution. The results presented can be very useful in guiding future numerical work on these models and in revising previous analytic work on all their applications.

We examine the radial asymptotic behavior of spherically symmetric Lemaitre-Tolman-Bondi dust models by looking at their covariant scalars along radial rays, which are spacelike geodesics parametrized by proper length $\ell$, orthogonal to the 4-velocity and to the orbits of SO(3). By introducing quasi-local scalars defined as integral functions along the rays, we obtain a complete and covariant representation of the models, leading to an initial value parametrization in which all scalars can be given by scaling laws depending on two metric scale factors and two basic initial value functions. Considering regular "open" LTB models whose space slices allow for a diverging $\ell$, we provide the conditions on the radial coordinate so that its asymptotic limit corresponds to the limit as $\ell\to\infty$. The "asymptotic state" is then defined as this limit, together with asymptotic series expansion around it, evaluated for all metric functions, covariant scalars (local and quasi-local) and their fluctuations. By looking at different sets of initial conditions, we examine and classify the asymptotic states of parabolic, hyperbolic and open elliptic models admitting a symmetry center. We show that in the radial direction the models can be asymptotic to any one of the following spacetimes: FLRW dust cosmologies with zero or negative spatial curvature, sections of Minkowski flat space (including Milne's space), sections of the Schwarzschild--Kruskal manifold or self--similar dust solutions.

We examine the dynamics of a self--gravitating magnetized neutron gas as a source of a Bianchi I spacetime described by the Kasner metric. The set of Einstein-Maxwell field equations can be expressed as a dynamical system in a 4-dimensional phase space. Numerical solutions of this system reveal the emergence of a point--like singularity as the final evolution state for a large class of physically motivated initial conditions. Besides the theoretical interest of studying this source in a fully general relativistic context, the resulting idealized model could be helpful in understanding the collapse of local volume elements of a neutron gas in the critical conditions that would prevail in the center of a compact object.

We introduce quasi-local integral scalar variables for the study of spherically symmetric Lemaitre-Tolman-Bondi (LTB) dust models. Besides providing a covariant, and theoretically appealing, interpretation for the parameters of these models, these variables allow us to study their dynamics (in their full generality) by means of fluid flow evolution equations that can be handled with simple numerical techniques and has a significant potential for astrophysical and cosmological applications. These evolution equations can also be understood in the framework of a gauge invariant and covariant formalism of spherical non-linear perturbations on a FLRW background. The covariant time splitting associated with the new variables leads, in a natural way, to rephrase the known analytic solutions within an initial value framework in which covariant scalars are given by simple scaling laws. By using this re-parametrization of the analytic solutions, we re-examine and provide an alternative outlook to various theoretical issues already treated in the literature: regularity conditions, an Omega parameter, as well as the fitting of a given LTB model to radial profiles of density or velocity at different cosmic times. Other theoretical issues and numeric applications will be examined in separate articles.

We introduce quasi--local (QL) scalar variables in spherically symmetric LTB models. If the QL scalars are defined as functionals, they become weighed averages that generalize the standard proper volume averages on space slices orthogonal to the 4-velocity. We examine the connection between QL functions and functionals and the "back-reaction" term $Q$ in the context of Buchert's scalar averaging formalism. With the help of the QL scalars we provide rigorous proof that back--reaction is positive for (i) all LTB models with negative and asymptotically negative spatial curvature, and (ii) models with positive curvature decaying to zero asymptotically in the radial direction. We show by means of qualitative, but robust, arguments that generic LTB models exist, either with clump or void profiles, for which an "effective" acceleration associated with Buchert's formalism can mimic the effects of dark energy.

We provide a comprehensive analytic study (rigorous and qualitative) of the conditions for the existence of a a positive kinematic back reaction term $\QQ>0$, in the context of Buchert's scalar averaging formalism applied to spherically symmetric Lema\^itre-Tolman-Bondi (LTB) dust solutions in which averaging domains are given as spherical comoving regions containing a symmetry center. We introduce proper volume and quasi-local average functionals and functions in order to examine the conditions for $\QQ\geq 0$, and in the process we also explore the relation between back reaction, spatial curvature and binding energy for a wide variety of LTB configurations. The back reaction term is positive for all "hyperbolic" regular domains with negative spatial curvature, either in the full radial range or in the radial asymptotic range. This result is also valid if these domains contain an inner "elliptic" region with positive curvature undergoing local collapse. For some cases in which positive spatial curvature decreases asymptotically, the conditions for a positive back reaction can still be met but seem to be more restrictive. Since $\QQ>0$ is a necessary condition for a positive "effective" acceleration that would mimic the effect of dark energy (in the context of Buchert's formalism), we examine this issue in LTB models in a follow up paper (part II).

We use the Szekeres inhomogeneous relativistic models in order to fit supernova combined data sets. We show that with a choice of the spatial curvature function that is guided by current observations, the models fit the supernova data almost as well as the LCDM model without requiring a dark energy component. The Szekeres models were originally derived as an exact solution to Einstein's equations with a general metric that has no symmetries and are regarded as good candidates to model the true lumpy universe that we observe. The null geodesics in these models are not radial. The best fit model found is also consistent with the requirement of spatial flatness at CMB scales. The first results presented here seem to encourage further investigations of apparent acceleration using various inhomogeneous models and other constraints from CMB and large structure need to be explored next.

A numerical approach is considered for spherically symmetric spacetimes that generalize Lemaitre-Tolman-Bondi dust solutions to nonzero pressure ("LTB spacetimes"). We introduce quasi-local (QL) variables that are covariant LTB objects satisfying evolution equations of Friedman-Lemaitre-Robertson-Walker (FLRW) cosmologies. We prove rigorously that relative deviations of the local covariant scalars from the QL scalars are non-linear, gauge invariant and covariant perturbations on a FLRW formal "background" given by the QL scalars. The dynamics of LTB spacetimes is completely determined by the QL scalars and these exact perturbations. Since LTB spacetimes are compatible with a wide variety of "equations of state", either single fluids or mixtures, a large number of known solutions with dark matter and dark energy sources in a FLRW framework (or with linear perturbations) can be readily examined under idealized but non-trivial inhomogeneous conditions. Coordinate choices and initial conditions are derived for a numerical treatment of the perturbation equations, allowing us to study non-linear effects in a variety of phenomena, such as gravitational collapse, non-local effects, void formation, dark matter and dark energy couplings and particle creation. In particular, the embedding of inhomogeneous regions can be performed by a smooth matching with a suitable FLRW solution, thus generalizing the Newtonian "top hat" models that are widely used in astrophysical literature. As examples of the application of the formalism, we examine numerically the formation of a black hole in an expanding Chaplygin gas FLRW universe, as well as the evolution of density clumps and voids in an interactive mixture of cold dark matter and dark energy.

We derive an expression for the entropy of a dark matter halo described by a Navarro-Frenk-White model with a core. The comparison of this entropy with the one of dark matter at the freeze-out era allows us to constraint the parameter space in mSUGRA models. Moreover, joining these constraints with the ones obtained from the usual abundance criteria and demanding both criteria to be consistent with the 2$\sigma$ bounds for the abundance of dark matter: $0.112\leq\Omega_{DM}h^2\leq0.122$, we are able to clearly discriminate validity regions among the values of $\tan \beta$, one of the parameters of the mSUGRA model. We found that for the explored regions of the parameter space, small values of tan$\beta$ are not favored; only for tan$\beta\simeq50$ are both criteria significantly consistent. In the region where both criteria are consistent we also found a lower bound for the neutralino mass, $m_{\chi}\geq 141$ GeV.

We consider a magnetized degenerate gas of fermions as the matter source of a homogeneous but anisotropic Bianchi I spacetime with a Kasner--like metric. We examine the dynamics of this system by means of a qualitative and numerical study of Einstein-Maxwell field equations which reduce to a non--linear autonomous system. For all initial conditions and combinations of free parameters the gas evolves from an initial anisotropic singularity into an asymptotic state that is completely determined by a stable attractor. Depending on the initial conditions the anisotropic singularity is of the ``cigar'' or ``plate'' types.

We use the Bianchi-I spacetime to study the local dynamics of a magnetized self-gravitating Fermi gas. The set of Einstein-Maxwell field equations for this gas becomes a dynamical system in a 4-dimensional phase space. We consider a qualitative study and examine numeric solutions for the degenerate zero temperature case. All dynamic quantities exhibit similar qualitative behavior in the 3-dimensional sections of the phase space, with all trajectories reaching a stable attractor whenever the initial expansion scalar H_{0} is negative. If H_{0} is positive, and depending on initial conditions, the trajectories end up in a curvature singularity that could be isotropic(singular "point") or anisotropic (singular "line"). In particular, for a sufficiently large initial value of the magnetic field it is always possible to obtain an anisotropic type of singularity in which the "line" points in the same direction of the field.

We examine numerically and qualitatively the Lema\^\i tre--Tolman--Bondi (LTB) inhomogeneous dust solutions as a 3--dimensional dynamical system characterized by six critical points. One of the coordinates of the phase space is an average density parameter, $<\Omega>$, which behaves as the ordinary $\Omega$ in Friedman-Lema\^\i tre--Robertson--Walker (FLRW) dust spacetimes. The other two coordinates, a shear parameter and a density contrast function, convey the effects of inhomogeneity. As long as shell crossing singularities are absent, this phase space is bounded or it can be trivially compactified. This space contains several invariant subspaces which define relevant particular cases, such as: ``parabolic'' evolution, FLRW dust and the Schwarzschild--Kruskal vacuum limit. We examine in detail the phase space evolution of several dust configurations: a low density void formation scenario, high density re--collapsing universes with open, closed and wormhole topologies, a structure formation scenario with a black hole surrounded by an expanding background, and the Schwarzschild--Kruskal vacuum case. Solution curves start expanding from a past attractor (source) in the plane $<\Omega>=1$, associated with self similar regime at an initial singularity. Depending on the initial conditions and specific configurations, the curves approach several saddle points as they evolve between this past attractor and other two possible future attractors: perpetually expanding curves terminate at a line of sinks at $<\Omega>=0$, while collapsing curves reach maximal expansion as $<\Omega>$ diverges and end up in sink that coincides with the past attractor and is also associated with self similar behavior.

We derive an expression for the entropy of a present dark matter halo described by a Navarro-Frenk-White modified model with a central core. The comparison of this entropy with the one of the halo at the freeze-out era allows us to obtain an expression for the relic abundance of neutralinos, which in turn is used to constrain the parameter space in mSUGRA models, when used with the WMAP observations. Moreover, by joning these results with the ones obtained from the usual abundance criteria, we are able to clearly discriminate validity regions among tan beta values of the mSUGRA model, by demanding both criteria to be consistent with the 2 sigma bounds of the WMAP observations for the relic density: 0.112 < Omega h^2 < 0.122. We found that for sign mu positive, small values of tan beta are not favored; only for tan beta ~ 50 are both criteria significantly consistent. The use of both criteria also allows us to put a lower bound on the neutralino mass of > 151 GeV.

We show that if the total internal energy of a black hole is constructed as the sum of $N$ photons all having a fixed wavelength chosen to scale with the Schwarzschild radius as $\lambda=\alpha R_{s}$, then $N$ will scale with $R_{s}^{2}$. A statistical mechanical calculation of the configuration proposed yields (\alpha = 4 \pi^2 / \ln(2)) and a total entropy of the system $S=k_{B}N \ln(2)$, in agreement with the Bekenstein entropy of a black hole . It is shown that the critical temperature for Bose-Einstein condensation for relativistic particles of $\lambda=\alpha R_{s}$ is always well below the Hawking temperature of a black hole, in support of the proposed internal configuration. We then examine our results from the point of view of recent loop quantum gravity ideas and find that a natural consistency of both approaches appears. We show that the Jeans criterion for gravitational instability can be generalised to the special and general relativistic regimes and holds for any type of mass--energy distribution.

Motivated by the possible conflict between the Navarro-Frenk-White(NFW) model predictions for the dark matter contents of galactic systems and its correlation with baryonic surface density, we will explore an alternative paradigm for the description of dark matter halos. Such an alternative emerges from Tsallis' non-extensive thermodynamics applied to self-gravitating systems and leads to the so-called ``stellar polytrope'' (SP) model. We consider that this could be a better approach to real structures rather than the isothermal model, given the fact that the first one takes into account the non-extensivity of energy and entropy present in these type of systems characterized by long-range interactions. We compare a halo based on the Navarro-Frenk-White (NFW) and one which follows the SP description. Analyzing the dark matter contents estimated by means of global physical parameters of galactic disks, obtained from a sample of actual galaxies, with the ones of the unobserved dark matter halos, we conclude that the SP model is favored over the NFW model in such a comparison.

We explore the use of two criteria to constraint the allowed parameter space in mSUGRA models. Both criteria are based in the calculation of the present density of neutralinos as dark matter in the Universe. The first one is the usual ``abundance'' criterion which is used to calculate the relic density after the ``freeze-out'' era. To compute the relic density we used the numerical public code micrOMEGAs. The second criterion applies the microcanonical definition of entropy to a weakly interacting and self-gravitating gas evaluating then the change in the entropy per particle of this gas between the ``freeze-out'' era and present day virialized structures. An ``entropy-consistency'' criterion emerges by comparing theoretical and empirical estimates of this entropy. The main objective of our work is to determine for which regions of the parameter space in the mSUGRA model are both criteria consistent with the 2$\sigma$ bounds according to WMAP for the relic density: $0.0945<\Omega_{CDM}h^2<0.1287$. As a first result, we found that for $A_0=0$, sgn$\mu=+$, small values of tan$\beta$ are not favored; only for tan$\beta\simeq50$ are both criteria significantly consistent.

We study a dark matter halo model from two points of view: the ``stellar polytrope'' (SP) model coming from Tsallis' thermodynamics, and the one coming from the Navarro-Frenk-White (NFW) paradigm. We make an appropriate comparison between both halo models and analyzing the relations between the global physical parameters of observed galactic disks, coming from a sample of actual galaxies, with the ones of the unobserved dark matter halos, we conclude that the SP model is favored over the NFW model in such a comparison.

We present an alternative for the description of galactic halos based on Tsallis' non--extensive entropy formalism; on this scheme, halos are stellar polytropes characterized by three parameters, the central density, $\rho_c$, the central velocity dispersion, $\sigma_c$ and the polytropic index, $n$. To evaluate these parameters we take the Navarro-Frenk-White paradigm as a comparative model and make the following assumptions: both halo models must have the same virial mass, the same total energy and the same maximal velocity. These three conditions fix all the parameters for a given stellar polytrope allowing us to compare both halo models. The halos studied have virial masses on the range $10^{12}-10^{15} M_\odot$, and it was found after the analysis that they are described, at all scales, by almost the same polytropic index, $n\approx 4.8$, implying an empirical estimation of Tsallis non--extensive parameter for this type of dynamical systems: $q\approx 1.3$.

We derive and analyze a class of spherically symmetric cosmological models whose source is an interactive mixture of inhomogeneous cold dark matter (DM) and a generic homogeneous dark energy (DE) fluid. If the DE fluid corresponds to a quintessense scalar field, the interaction term can be associated with a well motivated non--minimal coupling to the DM component. By constructing a suitable volume average of the DM component we obtain a Friedman evolution equation relating this average density with an average Hubble scalar, with the DE component playing the role of a repulsive and time-dependent $\Lambda$ term. Once we select an ``equation of state'' linking the energy density ($\mu$) and pressure ($p$) of the DE fluid, as well as a free function governing the radial dependence, the models become fully determinate and can be applied to known specific DE sources, such as quintessense scalar fields or tachyonic fluids. Considering the simple equation of state $p= (\gamma-1) \mu$ with $0\leq\gamma <2/3$, we show that the free parameters and boundary conditions can be selected for an adequate description of a local DM overdensity evolving in a suitable cosmic background that accurately fits current observational data. While a DE dominated scenario emerges in the asymptotic future, with total $\Omega$ and $q$ tending respectively to 1 and -1/2 for all cosmic observers, the effects of inhomogeneity and anisotropy yield different local behavior and evolution rates for these parameters in the local overdense region. We suggest that the models presented can be directly applied to explore the effects of various DE formalisms on local DM cosmological inhomogeneities.

In order to study their interplay with large scale cosmic evolution and with relativistic effects, such as gravitational lenses, quintessence sources or gravitational waves, we construct a post-Newtonian fluid framework for the "Navarro-Frenk-White'' (NFW) models of galactic halos that follow from N-body numerical simulations. Since these simulations are unable to resolve regions very near the halo center, the extrapolation of the fitting formula leads to a spherically averaged "universal'' density profile that diverges at the origin. We remove this inconvenient feature by replacing a small central region of the NFW halo with an interior Schwarzschild solution with constant density, continuously matched to the remaining NFW spacetime. A model of a single halo, as an isolated object with finite mass, follows by smoothly matching the NFW spacetime to a Schwarzschild vacuum exterior along the virial radius, the physical "cut-off'' customarily imposed, as the mass associated with NFW profiles diverges asymptotically. Numerical simulations assume weakly interacting collisionless particles, hence we suggest that NFW halos approximately satisfy an "ideal gas'' type of equation of state, where mass-density is the dominant rest--mass contribution to matter--energy, with the internal energy contribution associated with an anisotropic kinetic pressure. We show that, outside the central core, this pressure and the mass density roughly satisfy a polytropic relation. Since stellar polytropes are the equilibrium configurations in Tsallis' non-extensive formalism of Statistical Mechanics, we argue that NFW halos might provide a rough empirical estimate of the free parameter $q$ of Tsallis' formalism.

Using relativistic kinetic theory, we study spherically symmetric, static equilibrium configurations of a collisionless Maxwell-Boltzmann gas with non-standard self-interactions, modelled by an effective one--particle force. The resulting set of equilibrium conditions represents a generalization of the classical Tolman-Oppenheimer-Volkov equations. We specify these conditions for two types of Lorentz--like forces: one coupled to the 4-acceleration and the 4--velocity and the other one coupled to the Riemann tensor. We investigate the weak field limits in each case and show that they lead to various Newtonian type configurations that are different from the usual isothermal sphere characterizing the conventional Newtonian Maxwell--Boltzmann gas. These configurations could provide a plausible phenomenological and theoretical description of galactic dark matter halo structures. We show how the self--interaction may act phenomenologically as an effective cosmological constant and discuss possible connections with Modified Newtonian Dynamics(MOND).

We examine the possibility that galactic halos are collisionless ensembles of scalar field ``massive compact halo objects'' (MACHOs). Using mass constraints from MACHO microlensing and from theoretical arguments on halos made up of massive black holes, as well as demanding also that scalar MACHO ensambles of all scales do not exhibit gravothermal instability (as required by consistency with observations of LSB galaxies), we obtain the range: $m\alt 10^{-7} M_\odot$ or $30 M_\odot\alt m\alt 100 M_\odot$. The rather narrow mass range of large MACHOs seems to indicate that the ensambles we are suggesting should be probably made up of scalar MACHOs in the low mass range (``mini--MACHOs''). The proposed model allows one to consider a non--baryonic and non--thermal fundamental nature of dark matter, while at the same time keeping the same phenomenology of the CDM paradigm.

Applying the microcanonical definition of entropy to a weakly interacting and self--gravitating neutralino gas, we evaluate the change in the local entropy per particle of this gas between the freeze out era and present day virialized halo structures. An ``entropy consistency'' criterion emerges by comparing theoretical and empirical estimates of this entropy. We apply this criterion to the cases when neutralinos are mostly B-inos and mostly Higgsinos, in conjunction with the usual ``abundance'' criterion requiring that present neutralino relic density complies with $0.2<\Omega_{{\tilde\chi^1_0}} < 0.4$ for $h\simeq 0.65$. The joint application of both criteria reveals that a much better fitting occurs for the B-ino than for the Higgsino channels, so that the former seems to be a favored channel along the mass range of $155 \hbox{GeV} < m_{{\tilde\chi^1_0}} < 230 \hbox{GeV}$. These results are consistent with neutralino annihilation patterns that emerge from recent theoretical analysis on cosmic ray positron excess data reported by the HEAT collaboration. The suggested methodology can be applied to test other annihilation channels of the neutralino, as well as other particle candidates of thermal WIMP gas relics.

We show how an adequate post--Newtonian generalization can be obtained for Newtonian dark matter halos associated with an empiric density profile. Applying this approach to halos that follow from the well known numerical simulations of Navarro, Frenk and White (NFW), we derive all dynamical variables and show that NFW halos approximatelly follow an ideal gas type of equation of state which fits very well to a polytropic relation in the region outside the core. This fact suggests that ``outer'' regions of NFW halos might be related to equilibrium states in the non--extensive Statistical Mechanics formalism proposed by Tsallis.

This paper was withdrawn by the authors due to serious flaws in the physics that render the results meaningless.

We examine isothermal dark matter halos in hydrostatic equilibrium with a cosmological constant Lambda =Omega_\Lambda rho_{crit}c^2, where Omega_\Lambda=0.7, and rho_{crit} is the present value of the critical density with h=0.65. The Newtonian limit of General Relativity yields equilibrium equations that are different from those arising by merely coupling an ``isothermal sphere'' to the Lambda-field within a Newtonian framework. The conditions for the existence and stability of circular geodesic orbits show the existence of (I) an ``isothermal region'' (0<r<r_2), circular orbits are stable and all variables behave almost identically to those of an isothermal sphere; (II) an ``asymptotic region'' (r>r_1) dominated by the Lambda-field, where the Newtonian potential oscillates and circular orbits exist in disconnected patches of the domain of r; (III) a ``transition region'' (r_2<r<r_1), circular orbits exist but are unstable. We also find that no stable configuration exists with central density, rho_c, smaller than 2 Lambda, hence any galactic haloes which virialized at $z< 30$ in must have rho_c >0.008 M_\odot {pc}^3, in agreement with rotation curve studies of dwarf galaxies. Since r_2 marks the largest radius of a stable circular orbit, it provides a ``cut off'' radius. For current estimates of rho_c and velocity dispersion of galactic structures, this is around five times larger than the virialization radius. The effects of the Lambda$--field can hence be ignored in structure formation models, but could be significant in the dynamics of superclusters in the linear regime or of gravitational clustering at large scales.

Applying the microcanonical definition of entropy to a weakly interacting and self--gravitating neutralino gas, we evaluate the change in the local entropy per particle of this gas between the freeze out era and present day virialized halo structures. An ``entropy consistency'' criterion emerges by comparing the obtained theoretical entropy per particle of the virialized halos with an empirical entropy per particle given in terms of dynamical halo variables of actual galactic structures. We apply this criterion to the cases when neutralinos are mostly B-inos and mostly Higgsinos, in conjunction with the usual ``abundance'' criterion requiring that present neutralino relic density complies with $0.2<\Omega_{{\tilde\chi^1_0}} < 0.4$ for $h\simeq 0.65$. The joint application of both criteria reveals that a much better fitting occurs for the B-ino than for the Higgsino channels, so that the former seems to be a favored channel along the mass range of $150 \hbox{GeV} < m_{{\tilde\chi^1_0}} < 250 \hbox{GeV}$. These results are consistent with neutralino annihilation patterns that emerge from recent theoretical analysis on cosmic ray positron excess data reported by the HEAT collaboration. The suggested methodology can be applied to test other annihilation channels of the neutralino, as well as other particle candidates of thermal WIMP gas relics.

We examine, from both a qualitative and a numerical point of view, the evolution of Kantowski--Sachs cosmological models whose source is a mixture of a gas of weakly interacting massive particles (WIMP's) and a radiative gas made up of a ``tightly coupled'' mixture of electrons, baryons and photons from the end of nucleosynthesis up to the duration of radiative interactions. In this cosmic era annihilation processes are negligible, while the WIMP's only interact gravitationally with the radiative gas and the latter behaves as a single dissipative fluid that can be studied within a hydrodynamical framework. Applying the full transport equations of Extended Irreversible Thermodynamics, coupled with the field and balance equations, we obtain a set of governing equations that becomes an autonomous system of ordinary differential equations once the shear viscosity relaxation time is specified; assuming the latter to be proportional to the Hubble time, the qualitative analysis indicates that models begin in the radiation dominated epoch close to an isotropic equilibrium point. We show how the relaxation time governs the relaxation of the models towards an equilibrium photon entropy, leading to ``near-Eckart'' and transient regimes associated with `abrupt' and `smooth' relaxation processes, respectively. Assuming the WIMP particle to be a super-symmetric neutralino of a 100$ GeV, the numerical analysis reveals that a physically plausible evolution, compatible with a stable equilibrium state and with observed bounds on CMB anisotropies and neutralino abundance, is only possible for models characterized by initial conditions associated with nearly zero spatial curvature and total initial energy density very close to unity.

We consider simple hydrodynamical models of galactic dark matter in which the galactic halo is a self-gravitating and self-interacting gas that dominates the dynamics of the galaxy. Modeling this halo as a sphericaly symmetric and static perfect fluid satisfying the field equations of General Relativity, visible barionic matter can be treated as ``test particles'' in the geometry of this field. We show that the assumption of an empirical ``universal rotation curve'' that fits a wide variety of galaxies is compatible, under suitable approximations, with state variables characteristic of a non-relativistic Maxwell-Boltzmann gas that becomes an isothermal sphere in the Newtonian limit. Consistency criteria lead to a minimal bound for particle masses in the range $30 \hbox{eV} \leq m \leq 60 \hbox{eV}$ and to a constraint between the central temperature and the particles mass. The allowed mass range includes popular supersymmetric particle candidates, such as the neutralino, axino and gravitino, as well as lighter particles ($m\approx$ keV) proposed by numerical N-body simulations associated with self-interactive CDM and WDM structure formation theories.

We re-examine the Lem\^aitre-Tolman-Bondi (LTB) solutions with a dust source admitting symmetry centers. We consider as free parameters of the solutions the initial value functions: $Y_i$, $\rho_i$ and $\Ri$, obtained by restricting the curvature radius, $Y\equiv \sqrt{g_{\theta\theta}}$, the rest mass density, $\rho$, and the 3-dimensional Ricci scalar of the rest frames, $\R$, to an arbitrary regular Cauchy hypersurface, $\Ti$, marked by constant cosmic time ($t=t_i$). Using $Y_i$ to fix the radial coordinate and the topology (homeomorphic class) of $\Ti$, and scaling the time evolution in terms of an adimensional scale factor $y=Y/Y_i$, we show that the dynamics, regularity conditions and geometric features of the models are determined by $\rho_i$, $\Ri$ and by suitably constructed volume averages and contrast functions expressible in terms of invariant scalars defined in $\Ti$. These quantities lead to a straightforward characterization of initial conditions in terms of the nature of the inhomogeneity of $\Ti$, as density and/or curvature overdensities (``lumps'') and underdensities (''voids'') around a symmetry center. In general, only models with initial density and curvature lumps evolve without shell crossing singularities, though special classes of initial conditions, associated with a simmultaneous big bang, allow for a regular evolution for initial density and curvature voids. Specific restrictions are found so that a regular evolution for $t>t_i$ is possible for initial voids. A step-by-step guideline is provided for using the new variables in the construction of LTB models and for plotting all relevant quantities.

A new class of solutions of Einstein's equations is derived providing a physically plausible hydrodynamic description of cosmological matter in the radiation era of the Universe expansion. The solutions are characterized by the LTB metric with dissipative fluid sources, subjected to the following conditions: (i) tthe nonequilibrium state variables satisfy the equations of state of a mixture of relativistic and non-relativistic ideal gases, where the internal energy of the latter has been neglected, (ii) the particle numbers of the mixture components are independently conserved, (iii) the viscous stress is consistent with the transport equation and entropy balance law of extended irreversible thermodynamics with the coefficients of shear viscosity provided by kinetic theory. The Jeans mass at decoupling is of the same order of magnitude as that of baryons dominated models.

We consider the Lema\^{\i}tre-Tolman-Bondi metric with an inhomogeneous viscous fluid source satisfying the equation of state of an interactive mixture of radiation and matter. Assuming conditions prior to the decoupling era, we apply Extended Irreversible Thermodynamcs (EIT) to this mixture. Using the full transport equation of EIT we show that the relaxation time of shear viscosity can be several orders of magnitude larger than the Thomson collision time between photons and electrons. A comparison with the ``truncated'' transport equation for these models reveals that the latter cannot describe properly the decoupling of matter and radiation

A physicaly reasonable interpretation is provided for the perfect fluid, sphericaly symmetric, conformally flat ``Stephani Universes''. The free parameters of this class of exact solutions are determined so that the ideal gas relation $p=n k T$ is identicaly fulfiled, while the full equation of state of a classical monatomic ideal gas and a matter-radiation mixture holds up to a good approximation in a near dust, matter dominated regime. Only the models having spacelike slices with positive curvature admit a regular evolution domain that avoids an unphysical singularity. In the matter dominated regime these models are dynamicaly and observationaly indistinguishable from ``standard'' FLRW cosmology with a dust source.

We derive a new class of exact solutions of Einstein's equations providing a physically plausible hydrodynamical description of cosmological matter in the radiative era ($10^6 K > T > 10^3 K$), between nucleosynthesis and decoupling. The solutions are characterized by the Lema\^{\i}tre-Tolman -Bondi metric with a viscous fluid source, subjected to the following conditions: (a) the equilibrium state variables satisfy the equation of state of a mixture of an ultra-relativistic and a non-relativistic ideal gases, where the internal energy of the latter has been neglected, (b) the particle numbers of the mixture components are independently conserved, (c) the viscous stress is consistent with the transport equation and entropy balance law of Extended Irreversible Thermodynamics, with the coefficient of shear viscosity provided by Kinetic Theory for the `radiative gas' model. The fulfilment of (a), (b) and (c) restricts initial conditions in terms of an initial value function, $\Delta_i^{(s)}$, related to the average of spatial gradients of the fluctuations of photon entropy per baryon in the initial hypersurface. Constraints on the observed anisotropy of the microwave cosmic radiation and the condition that decoupling occurs at $T=T_{_D}\approx 4\times 10^3$ K yield an estimated value: $|\Delta_i^{(s)}|\approx 10^{-8}$ which can be associated with a bound on promordial entropy fluctuations. The Jeans mass at decoupling is of the same order of magnitude as that of baryon dominated perturbation models ($\approx 10^{16} M_\odot$)

We derive a new class of exact solutions characterized by the Szekeres-Szafron metrics (of class I), admitting in general no isometries. The source is a fluid with viscosity but zero heat flux (adiabatic but irreversible evolution) whose equilibrium state variables satisfy the equations of state of: (a) ultra-relativistic ideal gas, (b) non-relativistic ideal gas, (c) a mixture of (a) and (b). Einstein's field equations reduce to a quadrature that is integrable in terms of elementary functions (cases (a) and (c)) and elliptic integrals (case (b)). Necessary and sufficient conditions are provided for the viscous dissipative stress and equilibrium variables to be consistent with the theoretical framework of Extended Irreversible Thermodynamics and Kinetic Theory of the Maxwell-Boltzmann and radiative gases. Energy and regularity conditions are discussed. We prove that a smooth matching can be performed along a spherical boundary with a FLRW cosmology or with a Vaidya exterior solution. Possible applications are briefly outlined.

New exact solutions emerge by replacing the dust source of the Lem\^aitre-Tolman-Bondi metrics with a viscous fluid satisfying the monatomic gas equation of state. The solutions have a consistent thermodynamical interpretation. The most general transport equation of Extended Irreversible Thermodynamics is satisfied, with phenomenological coefficients bearing a close resemblance to those characterizing a non relativistic Maxwell-Bolzmann gas.