In this paper we investigate the cosmological dynamics of geometric inflation by means of the tools of the dynamical systems theory. We focus in the study of two explicit models where it is possible to sum the infinite series of higher curvature corrections that arises in the formalism. These would be very interesting possibilities since, if regard gravity as a quantum effective theory, a key feature is that higher powers of the curvature invariants are involved at higher loops. Hence, naively, consideration of the whole infinite tower of curvature invariants amounts to consideration of all of the higher order loops. The global dynamics of these toy models in the phase space is discussed and the quantum origin of primordial inflation is exposed.

In this paper we investigate the asymptotic dynamics of inflationary cosmological models that are based in scalar-tensor theories of gravity. Our main aim is to explore the global structure of the phase space in the framework of single-field inflation models. For this purpose we make emphasis in the adequate choice of the variables of the phase space. Our results indicate that, although single-field inflation is generic in the sense that the corresponding critical point in the phase space exists for a wide class of potentials, along given phase space orbits -- representing potential cosmic histories -- the occurrence of the inflationary stage is rather dependent on the initial conditions. We have been able to give quantitative estimates of the relative probability (RP) for initial conditions leading to slow-roll inflation. For the non-minimal coupling model with the $\phi^2$-potential our rough estimates yield to an almost vanishing relative probability: $10^{-13}\,\%\lesssim RP\ll 10^{-8}\,\%$. These bonds are greatly improved in the scalar-tensor models, including the Brans-Dicke theory, where the relative probability $1\,\%\lesssim RP\leq 100\,\%$. Hence slow-roll inflation is indeed a natural stage of the cosmic expansion in Brans-Dicke models of inflation. It is confirmed as well that the dynamics of vacuum Brans-Dicke theories with arbitrary potentials are non-chaotic.

In this paper we investigate the cosmological dynamics of an up to cubic curvature correction to General Relativity (GR) known as Cosmological Einsteinian Cubic Gravity (CECG), whose vacuum spectrum consists of the graviton exclusively and its cosmology is well-posed as an initial value problem. We are able to uncover the global asymptotic structure of the phase space of this theory. It is revealed that an inflationary matter-dominated bigbang is the global past attractor which means that inflation is the starting point of any physically meaningful cosmic history. Given that higher order curvature corrections to GR are assumed to influence the cosmological dynamics at early times -- high energies/large curvature limit -- the late-time inflation can not be a consequence of the up to cubic order curvature modifications. We confirm this assumption by showing that late-time acceleration of the expansion in the CECG model is possible only if add a cosmological constant term.

In the bibliography a certain confusion arises in what regards to the classification of the gravitational theories into scalar-tensor theories and general relativity with a scalar field either minimally or non-minimally coupled to matter. Higher-derivatives Horndeski and beyond Horndeski theories that at first sight do not look like scalar-tensor theories only add to the confusion. To further complicate things, the discussion on the physical equivalence of the different conformal frames in which a given scalar-tensor theory may be formulated, makes even harder to achieve a correct classification. In this paper we propose a specific criterion for an unambiguous identification of scalar-tensor theories and discuss its impact on the conformal transformations issue. The present discussion carries not only pedagogical but also scientific interest since an incorrect classification of a given theory as a scalar-tensor theory of gravity may lead to conceptual issues and to the consequent misunderstanding of its physical implications.

In this paper we show that an equivalence between Horndeski and beyond Horndeski theories and general relativity with an effective imperfect fluid can be formally established. The formal equivalence is discussed for several particular cases of interest. Working in the cosmological framework, it is shown that, while the effective stress-energy tensor of viable Horndeski theories is formally equivalent to that of an imperfect fluid with anisotropic stresses and vanishing heat flux vector, the effective stress-energy tensor of beyond Horndeski theories is equivalent to the one of a perfect fluid instead.

Scalar fields have played an important role in the development of the fundamental theories of physics as well as in other branches of physics such as gravitation and cosmology. For a long time these escaped detection until 2012 year when the Higgs boson was observed for the first time. Since then alternatives to the general theory of relativity like the Brans-Dicke theory, scalar-tensor theories of gravity and their higher derivative generalizations -- collectively known as Horndeski theories -- have acquired renewed interest. In the present review we discuss on several selected topics regarding these theories, mainly from the theoretical perspective but with due mention of the observational aspect. Among the topics covered in this review we pay special attention to the following: 1) the asymptotic dynamics of cosmological models based in the Brans-Dicke, scalar-tensor and Horndeski theories, 2) inflationary models, extended quintessence and the Galileons, with emphasis in causality and stability issues, 3) the chameleon and Vainshtein screening mechanisms that may allow the elusive scalar field to evade the tight observational constraints implied by the solar system experiments, 4) the conformal frames conundrum with a brief discussion on the disformal transformations and 5) the role of Weyl symmetry and scale invariance in the gravitation theories. The review is aimed at specialists as well as at non-specialists in the subject, including postgraduate students.

Here we follow the mainstream of thinking about physical equivalence of different representations of a theory, regarded as the consequence of invariance of the laws of physics -- represented by an action principle and the derived motion equations -- under given transformations; be it coordinate, gauge or conformal transformations. Accordingly the conformal transformations' issue is discussed by invoking the assumed invariance of the laws of physics -- in particular the laws of gravity -- under conformal transformations of the metric. It is shown that Brans-Dicke and scalar-tensor theories are not well-suited to address physical equivalence of the conformal frames since the corresponding laws of gravity are not invariant under the conformal transformations or Weyl rescalings. The search for conformal symmetry leads us to explore the physical consequences of Weyl-invariant theories of gravity, that represent a natural arena where to discuss on physical equivalence of the conformally related representations. We show that conformal invariance of the action of a (supposedly conformal invariant) theory and of the derived motion equations is not enough to ensure actual Weyl invariance. It is required, also, that the underlying geometrical structure of the background spacetime be, at least, Weyl-integrable. Otherwise, if assume (as usual) spacetimes of (pseudo)Riemannian geometrical structure, the resulting -- apparently conformal invariant -- theory is anomalous in that, only massless matter fields can be consistently coupled. Gauge freedom, a distinctive feature of actually Weyl-invariant theories of gravity, leads to very unusual consequences.

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.

In this paper we investigate the so called "phantom barrier crossing" issue in a cosmological model based in the scalar-tensor theory with non-minimal derivative coupling to the Einstein's tensor. Special attention will be paid to the physical bounds on the squared sound speed. The numeric results are geometrically illustrated by means of a qualitative procedure of analysis that is based on the mapping of the orbits in the phase plane onto the surfaces that represent physical quantities in the extended phase space, that is: the phase plane complemented with an additional dimension relative to the given physical parameter. We find that the cosmological model based in the non-minimal derivative coupling theory -- this includes both the quintessence and the pure derivative coupling cases -- has serious causality problems related with superluminal propagation of the scalar and tensor perturbations. Even more disturbing is the finding that, despite that the underlying theory is free of the Ostrogradsky instability, the corresponding cosmological model is plagued by the Laplacian (classical) instability related with negative squared sound speed. This instability leads to an uncontrollable growth of the energy density of the perturbations that is inversely proportional to their wavelength. We show that independent of the self-interaction potential, for the positive coupling the tensor perturbations propagate superluminally, while for the negative coupling a Laplacian instability arises. This latter instability invalidates the possibility for the model to describe the primordial inflation.

In this paper we apply the tools of the dynamical systems theory in order to uncover the whole asymptotic structure of the vacuum interactions of a galileon model with a cubic derivative interaction term. It is shown that, contrary to what occurs in the presence of background matter, the galileon interactions of vacuum appreciably modify the late-time cosmic dynamics. In particular, a local late-time attractor representing phantom behavior arises which is inevitably associated with a big rip singularity. It seems that the gravitational interactions of the background matter with the galileon screen the effects of the gravitational self-interactions of the galileon, thus erasing any potential modification of the late-time dynamics by the galileon vacuum processes. Unlike other galileon models inspired in the DGP scenario, self-accelerating solutions do not arise in this model.

In this paper we investigate the chameleon effect in the different conformal frames of the Brans--Dicke theory. Given that, in the standard literature on the subject, the chameleon is described in the Einstein frame almost exclusively, here we pay special attention to the description of this effect in the Jordan and in the string frames. It is shown that, in general, terrestrial and solar system bounds on the mass of the BD scalar field, and bounds of cosmological origin, are difficult to reconcile at once through a single chameleon potential. We point out that, in a cosmological context, provided that the effective chameleon potential has a minimum within a region of constant density of matter, the Brans--Dicke theory transmutes into general relativity with a cosmological constant, in that region. This result, however, can be only locally valid. In cosmological settings de Sitter--general relativity is a global attractor of the Brans--Dicke theory only for the quadratic potential $V(\phi)=M^2\phi^2$, or for potentials that asymptote to $M^2\phi^2$.

In this paper we seek for relevant information on the asymptotic cosmological dynamics of the Brans--Dicke theory of gravity for several self-interaction potentials. By means of the simplest tools of the dynamical systems theory, it is shown that the general relativity de Sitter solution is an attractor of the Jordan frame (dilatonic) Brans--Dicke theory only for the exponential potential $U(\vphi)\propto\exp\vphi$, which corresponds to the quadratic potential $V(\phi)\propto\phi^2$ in terms of the original Brans--Dicke field $\phi=\exp\vphi$, or for potentials which asymptote to $\exp\vphi$. At the stable de Sitter critical point, as well as at the stiff-matter equilibrium configurations, the dilaton is necessarily massless. We find bounds on the Brans--Dicke coupling constant $\omega_\textsc{bd}$, which are consistent with well-known results.

The theory of the dynamical systems is a very complex subject which has brought several surprises in the recent past in connection with the theory of chaos and fractals. The application of the tools of the dynamical systems in cosmological settings is less known in spite of the amount of published scientific papers on this subject. In this paper a -- mostly pedagogical -- introduction to the application in cosmology of the basic tools of the dynamical systems theory is presented. It is shown that, in spite of their amazing simplicity, these allow to extract essential information on the asymptotic dynamics of a wide variety of cosmological models. The power of these tools is illustrated within the context of the so called $\Lambda$CDM and scalar field models of dark energy. This paper is suitable for teachers, undergraduate and postgraduate students from physics and mathematics disciplines.

Models of the universe with arbitrary (non gravitational) interaction between the components of the cosmic fluid: the phantom energy and the background, are investigated. A general form of the interaction that is inspired in scalar-tensor theories of gravity is considered. No specific model for the phantom fluid is assumed. We concentrate our investigation on solutions that are free of the coincidence problem. We found a wide region in the parameter space where the solutions are free of the big rip singularity also. Physical arguments, together with arguments based on the analysis of the observational evidence, suggest that phantom models without big rip singularity might be preferred by Nature.

Here we investigate the cosmic dynamics of Friedmann-Robertson-Walker universes -- flat spatial sections -- which are driven by nonlinear electrodynamics (NLED) Lagrangians. We pay special attention to the check of the sign of the square sound speed since, whenever the latter quantity is negative, the corresponding cosmological model is classically unstable against small perturbations of the background energy density. Besides, based on causality arguments, one has to require that the mentioned small perturbations of the background should propagate at most at the local speed of light. We also look for the occurrence of curvature singularities. Our results indicate that several cosmological models which are based in known NLED Lagrangians, either are plagued by curvature singularities of the sudden and/or big rip type, or are violently unstable against small perturbations of the cosmological background -- due to negative sign of the square sound speed -- or both. In addition, causality issues associated with superluminal propagation of the background perturbations may also arise.

In a recent paper [Phys. Rev. D{\bf 85}, 023528 (2012)] the authors proposed a generalized Born-Infeld electrodynamics coupled to general relativity which produces a nonsingular bouncing universe. For a magnetic universe the resulting cosmic evolution inevitably interpolates between asymptotic de Sitter states. Here we shall show that (i) the above theory does not have the standard weak field Maxwell limit, (ii) a sudden curvature singularity -- not better than the big bang -- arises, (iii) the speed of sound squared is a negative quantity signaling instability against small perturbations of the background energy density, and that (iv) the conclusion about the inevitability of the asymptotic vacuum regime in a magnetic universe is wrong.

By means of a combined use of the type Ia supernovae and H(z) data tests, together with the study of the asymptotic properties in the equivalent phase space - through the use of the dynamical systems tools - we demonstrate that the bulk viscous matter-dominated scenario is not a good model to explain the accepted cosmological paradigm, at least, under the parametrization of bulk viscosity considered in this paper. The main objection against such scenarios is the absence of conventional radiation and matter-dominated critical points in the phase space of the model. This entails that radiation and matter dominance are not generic solutions of the cosmological equations, so that these stages can be implemented only by means of unique and very specific initial conditions, i. e., of very unstable particular solutions. Such a behavior is in marked contradiction with the accepted cosmological paradigm which requires of an earlier stage dominated by relativistic species, followed by a period of conventional non-relativistic matter domination, during which the cosmic structure we see was formed. Also, we found that the bulk viscosity is positive just until very late times in the cosmic evolution, around z < 1. For earlier epochs it is negative, been in tension with the local second law of thermodynamics.

We study the phase space asymptotics of the so called Veneziano ghost dark energy models. Models where the ghost field's energy density: i) $\rho_{ghost}\propto H$, and ii) $\rho_{ghost}\propto H+H^2$, are investigated. Both, cases with and without additional non-gravitational interaction between cold dark matter and ghost dark energy, are subject to scrutiny. We pay special attention to the choice of phase space variables leading to bounded and compact phase space so that no critical point of physical interest is missing. A rich asymptotic structure is revealed: depending on the kind of non-minimal coupling critical points associated with radiation dominance, matter dominance, cold dark matter/ghost dark energy scaling, and ghost dark energy dominance, are found. Past and future attractors, as well as saddle equilibrium points, are identified in the corresponding phase spaces.

In this paper we show, through the study of concrete examples, that, depending on the cosmic dynamics of the energy density of matter degrees of freedom living in the brane, Randall-Sundrum (RS) brane effects can be important not only at short distances (UV regime), but also at large cosmological scales (IR regime). Our first example relies on the study, by means of the dynamical systems tools, of a toy model based in a non-linear electrodynamics (NLED) Lagrangian. Then we show that other, less elaborated models, such as the inclusion of a scalar phantom field, and of a tachyon phantom field -- trapped in the brane -- produce similar results. The form of the RS correction seems to convert what would have been future attractors in general relativity into saddle points. The above "mixing of scales" effect is distinctive only of theories that modify the right-hand-side (matter part) of the Friedmann equation, so that, for instance, DGP-brane models do not show this feature.

In this paper we investigate to which extent noncommutativity, a intrinsically quantum property, may influence the Friedmann-Robertson-Walker cosmological dynamics at late times/large scales. To our purpose it will be enough to explore the asymptotic properties of the cosmological model in the phase space. Our recipe to build noncommutativity into our model is based in the approach of reference [Phys. Rev. Lett. {\bf 88} (2002) 161301], and can be summarized in the following steps: i) the Hamiltonian is derived from the Einstein-Hilbert action (plus a self-interacting scalar field action) for a Friedmann-Robertson-Walker spacetime with flat spatial sections, ii) canonical quantization recipe is applied, i. e., the minisuperspace variables are promoted to operators, and the WDW equation is written in terms of these variables, iii) noncommutativity in the minisuperspace is achieved through the replacement of the standard product of functions by the Moyal star product in the WDW equation, and, finally, iv) semi-classical cosmological equations are obtained by means of the WKB approximation applied to the (equivalent) modified Hamilton-Jacobi equation. We demonstrate, indeed, that noncommutative effects of the kind considered here, can be those responsible for the present speed up of the cosmic expansion.

We investigate in detail the asymptotic properties of tachyon cosmology for a broad class of self-interaction potentials. The present approach relies in an appropriate re-definition of the tachyon field, which, in conjunction with a method formerly applied in the bibliography in a different context, allows to generalize the dynamical systems study of tachyon cosmology to a wider class of self-interaction potentials beyond the (inverse) square-law one. It is revealed that independent of the functional form of the potential, the matter-dominated solution and the ultra-relativistic (also matter-dominated) solution, are always associated with equilibrium points in the phase space of the tachyon models. The latter is always the past attractor, while the former is a saddle critical point. For inverse power-law potentials $V\propto\phi^{-2\lambda}$ the late-time attractor is always the de Sitter solution, while for sinh-like potentials $V\propto\sinh^{-\alpha}(\lambda\phi)$, depending on the region of parameter space, the late-time attractor can be either the inflationary tachyon-dominated solution or the matter-scaling (also inflationary) phase. In general, for most part of known quintessential potentials, the late-time dynamics will be associated either with de Sitter inflation, or with matter-scaling, or with scalar field-dominated solutions.

We study the consequences of further modification of $f(R,R_{\mu\nu} R^{\mu\nu}, R_{\mu\nu\sigma\rho} R^{\mu\nu\sigma\rho}) / f(R)$-theories by means of the Dirac-Born-Infeld deformation procedure, which amounts to the replacement of $f$ by $\lambda(\sqrt{1+2f/\lambda}-1)$ (the free parameter $\lambda$ fixes an additional energy scale). We pay special attention to the definition of masses of the linearized propagating degrees of freedom since these are important to judge about the stability of the linearization around vacuum background spaces. In this context we discuss the subtleties associated with expanding $f(R,R_{\mu\nu} R^{\mu\nu},R_{\mu\nu\sigma\rho} R^{\mu\nu\sigma\rho})$-Lagrangians around maximally symmetric spaces of constant curvature, as well as with equivalence of the linearized Lagrangian to a scalar-tensor theory. Investigation of the consequences of applying the Dirac-Born-Infeld strategy to further modify quadratic theories, on the stability of de Sitter vacuum, as well as its impact on the cosmological dynamics, is the main concern of this paper.

The dynamics of a cosmological model fueled by scalar field dark matter with a cosh-like potential plus a cosmological constant is investigated in detail. It is revealed that the late-time attractor is always the de Sitter solution, and that, depending on the values of the free parameters, the oscillating solution of the scalar field -- modeling cold dark matter -- mediates between some early stage (say, the radiation-dominated solution) and the accelerating de Sitter attractor.

In this work we investigate the evolution of matter density perturbations for quintessence models with a self-interaction potential that is a combination of exponentials. One of the models is based on the Einstein theory of gravity, while the other is based on the Brans-Dicke scalar tensor theory. We constrain the parameter space of the models using the determinations for the growth rate of perturbations derived from data of the 2-degree Field Galaxy Redshift Survey.

In this paper we investigate quintom cosmologies with arbitrary potentials from the dynamical systems perspective. The dynamical systems analysis is complete in the sense that it includes the asymptotic regime where both scalar fields diverge, which proves to be particularly relevant in connection with the existence of tracking phases. The results of the present study indicate that the existence of phantom attractors is not generic: for quintom models there may exist either de Sitter attractors associated with the saddle points of the potential, or tracking attractors in the asymptotic regime where the scalar fields diverge.

We explore quintessence models of dark energy which exhibit non-minimal coupling between the dark matter and the dark energy components of the cosmic fluid. The kind of coupling chosen is inspired in scalar-tensor theories of gravity. We impose a suitable dynamics of the expansion allowing to derive exact Friedmann-Robertson-Walker solutions once the coupling function is given as input. Self-interaction potentials of single and double exponential types emerge as result of our choice of the coupling function. The stability and existence of the solutions is discussed in some detail. Although, in general, models with appropriated interaction between the components of the cosmic mixture are useful to handle the coincidence problem, in the present study the coincidence can not be evaded due to the choice of the solution generating ansatz.

A class of braneworld models can lead to phantom-like acceleration of the late universe, but without the need for any phantom matter. In the simplest models, the universe contains only cold dark matter and a cosmological constant. We generalize these models by introducing a quintessence field. The new feature in our models is that quintessence leads to a crossing of the phantom divide, $w=-1$. This is a purely gravitational effect, and there is no phantom instability. Furthermore, the Hubble parameter is always decreasing, and there is no big rip singularity in the future.

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.

We explore quintessence models of dark energy which exhibit non-minimal coupling between the dark matter and the dark energy components of the cosmic fluid. The kind of coupling chosen is inspired in scalar-tensor theories of gravity. We impose a suitable dynamics of the expansion allowing to derive exact Friedmann-Robertson-Walker solutions once the coupling function is given as input. Coupling functions that lead to self-interaction potentials of the single and double exponential types are the target of the present investigation. The stability and existence of the solutions is discussed in some detail. The models are tested against the observational evidence.

A linear relationship between the Hubble expansion parameter and the time derivative of the scalar field is assumed in order to derive exact analytic cosmological solutions to Einstein's gravity with two fluids: a barotropic perfect fluid of ordinary matter, together with a self-interacting scalar field fluid accounting for the dark energy in the universe. A priori assumptions about the functional form of the self-interaction potential or about the scale factor behavior are not neccessary. These are obtained as outputs of the assumed linear relationship between the Hubble expansion parameter and the time derivative of the scalar field. As a consequence only a class of exponential potentials and their combinations can be treated. The relevance of the solutions found for the description of the cosmic evolution are discussed in some detail. The possibility to have superaccelerated expansion within the context of normal quintessence models is also discussed.

In this work we present a model of the universe in which dark energy is modelled explicitely with both a dynamical quintessence field and a cosmological constant. Our results confirm the possibility of a future collapsing universe (for a given region of the parameter space), which is necessary for a consistent formulation of string theory and quantum field theory. We have also reproduced the measurements of modulus distance from supernovae with good accuracy.

In this work we explore a model of the universe in which dark energy is modelled explicitely with both a dynamical quintessence field (with a double exponential self-interaction potential) and a cosmological constant. For a given region of the parameter space, our results confirm the possibility of a collapsing universe, which is necessary for an adequate definition of both perturbative quantum field and string theories. We have also reproduced the measurements of modulus distance from supernovae with good accuracy.

In this talk we present a model of the universe in which dark energy is modelled explicitely with both a dynamical quintessence field and a cosmological constant. Our results confirm the possibility of a collapsing universe (for a given region of the parameter space), which is advantageous for an adequate formulation of both perturbative quantum field and string theories. We have also reproduced the measurements of modulus distance from supernovae with good accuracy.

We briefly comment on a paper by Rubano and Scudellaro [astro-ph/0103335] where they found general exact solutions for two classes of exponential potentials in a scalar field model for quintessence. In that paper the authors were led to some interesting conclusions after a proper choice of the integration constants. By using dimensionless variables we show that the integration constants can be found explicitly without additional assumptions. In consequence we revise some results and conclusions in that paper. We also reproduce observations for Type-1a supernovae with good accuracy.