We analyze the global behaviour of the growth index of cosmic inhomogeneities in an isotropic homogeneous universe filled by cold non-relativistic matter and dark energy (DE) with an arbitrary equation of state. Using a dynamical system approach, we find the critical points of the system. That unique trajectory for which the growth index $\gamma$ is finite from the asymptotic past to the asymptotic future is identified as the so-called heteroclinic orbit connecting the critical points $(\Omega_m=0,~\gamma_{\infty})$ in the future and $(\Omega_m=1,~\gamma_{-\infty})$ in the past. The first is an attractor while the second is a saddle point, confirming our earlier results. Further, in the case when a fraction of matter (or DE tracking matter) $\varepsilon \Omega^{\rm tot}_m$ remains unclustered, we find that the limit of the growth index in the past $\gamma_{-\infty}^{\varepsilon}$ does not depend on the equation of state of DE, in sharp contrast with the case $\varepsilon=0$ (for which $\gamma_{-\infty}$ is obtained). We show indeed that there is a mathematical discontinuity: one cannot obtain $\gamma_{-\infty}$ by taking $\lim_{\varepsilon \to 0} \gamma^{\varepsilon}_{-\infty}$ (i.e. the limits $\varepsilon\to 0$ and $\Omega^{\rm tot}_m\to 1$ do not commute). We recover in our analysis that the value $\gamma_{-\infty}^{\varepsilon}$ corresponds to tracking DE in the asymptotic past with constant $\gamma=\gamma_{-\infty}^{\varepsilon}$ found earlier.

We study gravitational collapse in K-essence model with shift symmetry. For these models, we have the formation of two types of horizons, event and sonic. For the particular case $K(X)=X+\beta X^2$ we found three different regimes. In the weak field regime the scalar field disperses to infinity, in the very strong regime both horizons form at the same time and finally for the intermediate regime, the sonic horizon could form first or both horizons form at the same time. The threshold of formation of black hole is found in the regime where the sonic horizon forms first. We observe a universal behavior with scaling parameter $\gamma\simeq 0.51$. Interestingly this universal behavior is already encoded in the sonic horizon much before the black hole forms and therefore the emergence of the event horizon.

In the present paper, we explore various gravitational aspects such as energy extraction (via the Penrose process and Superradiance), particle collisions around a $\mathcal{N}=2$, $U(1)^2$ dyonic rotating black hole (BH) in the gauged supergravity model. The impact of the rotation parameter ($a$) and the gauge coupling constant ($g$) on the behaviour of horizon and ergoregion of the BH is studied. It is of interest to note that, compared with the extremal Kerr BH, the gauge coupling constant, under certain constraints, can enhance the maximum efficiency of energy extraction by the Penrose process almost double. Under the same constraints, we can extract approximately 60.75\% of the initial mass energy from the BH which is noticeably higher in contrast to the extremal Kerr BH. The limit of energy extraction in terms of the local speeds of the fragments is also examined with the help of the Wald inequality. We identify an upper limit on the gauge coupling constant up to which the phenomenon of Superradiance is likely to occur. Finally, we computed the center-of-mass energy ($E_{CM}$) of two particles with the same rest masses moving in the equatorial plane of the BH. Our study also aims to sensitize $E_{CM}$ to the rotation parameter and the gauge coupling constant for extremal and nonextremal spacetime as well. Especially, for the extremal case, an infinitely large amount of $E_{CM}$ can be achieved closer to the horizon which allows the BH to serve as a more powerful Planck-energy-scale collider as compared to Kerr and any other generalized BHs in the Kerr family explored so far in general relativity. However, $E_{CM}$ for the nonextremal spacetime is shown to be finite and has an upper bound.

We perform here a global analysis of the growth index $\gamma$ behaviour from deep in the matter era till the far future. For a given cosmological model in GR or in modified gravity, the value of $\gamma(\Omega_{m})$ is unique when the decaying mode of scalar perturbations is negligible. However, $\gamma_{\infty}$, the value of $\gamma$ in the asymptotic future, is unique even in the presence of a nonnegligible decaying mode today. Moreover $\gamma$ becomes arbitrarily large deep in the matter era. Only in the limit of a vanishing decaying mode do we get a finite $\gamma$, from the past to the future in this case. We find further a condition for $\gamma(\Omega_{m})$ to be monotonically decreasing (or increasing). This condition can be violated inside general relativity (GR) for varying $w_{DE}$ though generically $\gamma(\Omega_{m})$ will be monotonically decreasing (like $\Lambda$CDM), except in the far future and past. A bump or a dip in $G_{\rm eff}$ can also lead to a significant and rapid change in the slope $\frac{d\gamma}{d\Omega_{m}}$. On a $\Lambda$CDM background, a $\gamma$ substantially lower (higher) than $0.55$ with a negative (positive) slope reflects the opposite evolution of $G_{\rm eff}$. In DGP models, $\gamma(\Omega_{m})$ is monotonically increasing except in the far future. While DGP gravity becomes weaker than GR in the future and $w^{DGP}\to -1$, we still get $\gamma_{\infty}^{DGP}= \gamma_{\infty}^{\Lambda CDM}=\frac{2}{3}$. In contrast, despite $G^{DGP}_{\rm eff}\to G$ in the past, $\gamma$ does not tend to its value in GR because $\frac{dG^{DGP}_{\rm eff}}{d\Omega_{m}}\Big|_{-\infty}\ne 0$.

The large scale structure bispectrum in the squeezed limit couples large with small scales. Since relativity is important at large scales and non-linear loop corrections are important at small scales, the proper calculation of the observed bispectrum in this limit requires a non-linear relativistic calculation. We compute the matter bispectrum in general relativity in the weak field approximation. The calculation is as involved as existing second-order results. We find several differences with the Newtonian calculation such as the non-cancellation of IR divergences, the need to renormalize the background, and the fact that initial conditions must be set at second order in perturbation theory. For the bispectrum, we find relativistic corrections to be as large as the newtonian result in the squeezed limit. In that limit relativistic one-loop contributions, which we compute for the first time, can be as large as tree level results and have the same $1/k^2$ dependence as a primordial local non-Gaussianity signal where $k$ is the momentum approaching zero. Moreover, we find the time dependence of the relavistic corrections to the bispectrum to be the same as that of a primordial non-Gaussianity signal.

Assuming a simple form for the growth index gamma(z) depending on two parameters gamma_0 = gamma(z=0) and gamma_1 = gamma'(z=0), we show that these parameters can be constrained using background expansion data. We explore systematically the preferred region in this parameter space. Inside General Relativity we obtain that models with a quasi-static growth index and gamma_1 = -0.02 are favoured. We find further the lower bounds gamma_0 > 0.53 and gamma_1 > -0.15 for models inside GR. Models outside GR having the same background expansion as LCDM and arbitrary gamma(z) with gamma_0 = gamma_0^{LCDM}, satisfy G_{eff,0}>G for gamma_1 > gamma_1^{LCDM}, and G_{eff,0}<G for gamma_1 < gamma_1^{LCDM}. The first models will cross downwards the value G_{eff}=G on very low redshifts z<0.3, while the second models will cross upwards G_{eff}=G in the same redshift range. This makes the realization of such modified gravity models even more problematic.

Recent analyses \cite{Nesseris:2017vor,Kazantzidis:2018rnb} have indicated that an effective Newton's constant $G_{\rm eff}(z)$ decreasing with redshift may relieve the observed tension between the Planck15 best fit $\Lambda$CDM cosmological background ({\it i.e.} Planck15/$\Lambda$CDM) and the corresponding $\Lambda$CDM background favored by growth $f\sigma_8$ and weak lensing data. We investigate the consistency of such a decreasing $G_{\rm eff}(z)$ with some viable scalar-tensor models and $f(R)$ theories. We stress that $f(R)$ theories generically can not lead to a decreasing $G_{\rm eff}(z)$ for any cosmological background. For scalar-tensor models we deduce that in the context of a $\Lambda$CDM cosmological background, a decreasing $G_{\rm eff}(z)$ is not consistent with a large Brans-Dicke parameter $\omega_{BD,0}$ today. This inconsistency remains and amplifies in the presence of a phantom dark energy equation of state parameter ($w < -1$). However it can be avoided for $w >-1$. We also find that any modified gravity model with the required decreasing $G_{\rm eff}(z)$ and $G_{{\rm eff},0}=G$, would have a characteristic signature in its growth index $\gamma$ with $0.61\lesssim \gamma_0\lesssim 0.69$ and large slopes $\gamma_0'$, $0.16\lesssim \gamma_0'\lesssim 0.4$, which is a characteristic signature of a decreasing (with $z$) $G_{\rm eff}(z)<G$ on small redshifts. This is a substantial departure today from the quasi-static behaviour in $\Lambda$CDM with $(\gamma_0,\gamma_0')\approx (0.55,-0.02)$.

We consider an Einstein-Scalar-Gauss-Bonnet gravitational theory, and argue that at early times the Ricci scalar can be safely ignored. We then demonstrate that the pure scalar-Gauss-Bonnet theory, with a quadratic coupling function, naturally supports inflationary -- de Sitter -- solutions. During inflation, the scalar field decays exponentially and its effective potential remains always bounded. The theory contains also solutions where these de Sitter phases possess a natural exit mechanism and are replaced by linearly expanding -- Milne -- phases.

In this work, we consider a generalised gravitational theory that contains the Einstein term, a scalar field and the quadratic Gauss-Bonnet term. We focus on the early-universe dynamics, and demonstrate that a simple choice of the coupling function between the scalar field and the Gauss-Bonnet term and a simplifying assumption regarding the role of the Ricci scalar can lead to new, analytical, elegant solutions with interesting characteristics. We first argue, and demonstrate in the context of two different models, that the presence of the Ricci scalar in the theory at early times, when the curvature is strong, does not affect the actual cosmological solutions. By considering therefore a pure scalar-GB theory with a quadratic coupling function we derive a plethora of interesting, analytic solutions: for a negative coupling parameter, we obtain inflationary, de Sitter-type solutions or expanding solutions with a de Sitter phase in their past and a natural exit mechanism at later times; for a positive coupling function, we find instead singularity-free solutions with no Big-Bang singularity. We show that the aforementioned solutions arise only for this particular choice of coupling function, a result that may hint to some fundamental role that this coupling function may hold in the context of an ultimate theory.

We study the structure of neutron stars in $f(R)=R+\alpha R^{2}$ theory of gravity (Starobinsky model), in an exact and non-perturbative approach. In this model, apart from the standard General Relativistic junction conditions, two extra conditions, namely the continuity of the curvature scalar and its first derivative needs to be satisfied. For an exterior Schwarzschild solution, the curvature scalar and its derivative has to be zero at the stellar surface. We show that for some equation of state (EoS) of matter, matching all conditions at the surface of the star is impossible. Hence the model brings two major fine-tuning problems: (i) only some particular classes of EoS are consistent with Schwarzschild at the surface and (ii) given that EoS, only a very particular set of boundary conditions at the centre of the star will satisfy the given boundary conditions at the surface. Hence we show that this model (and subsequently many other f(R) models where uniqueness theorem is valid) is highly unnatural, for the existence of compact astrophysical objects. This is because the EoS of a compact star should be completely determined by the physics of nuclear matter at high density and not the theory of gravity.

We consider symmetron model in a generalized background with a hope to make it compatible with dark energy. We observe a "no go" theorem at least in case of a conformal coupling. Being convinced of symmetron incapability to be dark energy, we try to retain its role for spontaneous symmetry breaking and assign the role of dark energy either to standard quintessence or $F(R)$ theory which are switched on by symmetron field in the symmetry broken phase. The scenario reduces to standard Einstein gravity in the high density region. After the phase transition generated by symmetron field, either the $F(R)$ gravity or the standard quintessence are induced in the low density region. we demonstrate that local gravity constraints and other requirements are satisfied although the model could generate the late-time acceleration of Universe.

We investigate cosmological behavior in the quasi-Dilaton non-linear massive gravity. We perform a detailed dynamical analysis and examine the stable late-time solutions relevant to late time cosmic acceleration. We demonstrate that a dark-energy dominated, cosmological-constant-like solution is a late attractor of the dynamics. We also analyze the evolution of the universe at intermediate times, showing that the observed epoch sequence can be easily obtained in the model under consideration. Furthermore, we study the non-singular bounce and the cosmological turnaround which can be realized for a region of the parameter space. Last but not least, we use observational data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO) and Cosmic Microwave Background (CMB), in order to constrain the parameters of the theory.

In this paper, we examine the cosmological viability of a light mass galileon field consistent with local gravity constraints. The minimal, L_3=\Box\phi(\partial_\mu \phi)^2, massless galileon field requires an additional term in order to give rise to a viable ghost free late time acceleration of Universe. The desired cosmological dynamics can either be achieved by incorporating an additional terms in the action such as (L_4,L_5) - the higher order galileon Lagrangians or by considering a light mass field a la galileon field potential. We analyse the second possibility and find that: (1) The model produces a viable cosmology in the regime where the non-linear galileon field is subdominant, (2) The Vainshtein mechanism operates at small scales where the non-linear effects become important and contribution of the field potential ceases to be significant. Also the small mass of the field under consideration is protected against strong quantum corrections thereby providing quantum stability to the system.

We derive field equations of Gauss-Bonnet gravity in 4 dimensions after dimensional reduction of the action and demonstrate that in this scenario Vainshtein mechanism operates in the flat spherically symmetric background. We show that inside this Vainshtein sphere the fifth force is negligibly small compared to the gravitational force. We also investigate stability of the spherically symmetric solution, clarify the vocabulary used in the literature about the hyperbolicity of the equation and the ghost-Laplacian stability conditions. We find superluminal behavior of the perturbation of the field in the radial direction. However, because of the presence of the non linear terms, the structure of the space-time is modified and as a result the field does not propagate in the Minkowski metric but rather in an "aether" composed by the scalar field $\pi(r)$. We thereby demonstrate that the superluminal behavior does not create time paradoxes thank to the absence of Causal Closed Curves. We also derive the stability conditions for Friedmann Universe in context with scalar and tensor perturbations.

We investigate the Weyl space-time extension of general relativity (GR) for studying the FLRW cosmology through focusing and defocusing of the geodesic congruences. We have derived the equations of evolution for expansion, shear and rotation in the Weyl space-time. In particular we consider the Starobinsky modification, $f(R)=R+\beta R^2-2\Lambda$, of gravity in the Einstein-Palatini formalism, which turns out to reduce to the Weyl integrable space-time (WIST) with the Weyl vector being a gradient. The modified Raychaudhuri equation takes the form of the Hill-type equation which is then analysed for studying formation of the caustics. In this model, it is possible to have a Big Bang singularity free cyclic Universe but unfortunately the periodicity turns out to be extremely short.

In this paper we examine the cosmological consequences of fourth order Galileon gravity. We carry out detailed investigations of the underlying dynamics and demonstrate the stability of one de Sitter phase. The stable de Sitter phase contains a Galileon field $\pi$ which is an increasing function of time (\dot{\pi}>0). Using the required suppression of the fifth force, supernovae, BAO and CMB data, we constrain parameters of the model. We find that the $\pi$ matter coupling parameter $\beta$ is constrained to small numerical values such that $\beta$<0.02. We also show that the parameters of the third and fourth order in the action (c_3,c_4) are not independent and with reasonable assumptions, we obtain constraints on them. We investigate the growth history of the model and find that the sub-horizon approximation is not allowed for this model. We demonstrate strong scale dependence of linear perturbations in the fourth order Galileon gravity.

In chameleon dark energy models, local gravity constraints tend to rule out parameters in which observable cosmological signatures can be found. We study viable chameleon potentials consistent with a number of recent observational and experimental bounds. A novel chameleon field potential, motivated by f(R) gravity, is constructed where observable cosmological signatures are present both at the background evolution and in the growth-rate of the perturbations. We study the evolution of matter density perturbations on low redshifts for this potential and show that the growth index today gamma_0 can have significant dispersion on scales relevant for large scale structures. The values of gamma_0 can be even smaller than 0.2 with large variations of gamma on very low redshifts for the model parameters constrained by local gravity tests. This gives a possibility to clearly distinguish these chameleon models from the Lambda-Cold-Dark-Matter model in future high-precision observations.

We consider the covariant galileon gravity taking into account the third order and fourth order scalar field Lagrangians L_3(\pi) and L_4(\pi) consisting of three and four $\pi$'s with four and five derivatives acting on them respectively. The background dynamical equations are set up for the system under consideration and the stability of the self accelerating solution is demonstrated in general setting. We extended this study to the general case of the fifth order theory. For spherically symmetric static background, we spell out conditions for suppression of fifth force effects mediated by the galileon field $\pi$. We study the field perturbations in the fixed background and investigate conditions for their causal propagation. We also briefly discuss metric fluctuations and derive evolution equation for matter perturbations in galileon gravity.

In this paper, we carry out a study of viable cosmological models in $f(R)$-gravity at the background level. We use observable parameters like $\Omega$ and $\gamma$ to form autonomous system of equations and show that the models under consideration exhibit two different regimes in their time evolution, namely, a phantom phase followed by a quintessence like behavior. We employ statefinder parameters to emphasize a characteristic discriminative signature of these models.

We study the growth of matter density perturbations delta_m for a number of viable f(R) gravity models that satisfy both cosmological and local gravity constraints, where the Lagrangian density f is a function of the Ricci scalar R. If the parameter m=Rf_{,RR}/f_{,R} today is larger than the order of 10^{-6}, linear perturbations relevant to the matter power spectrum evolve with a growth rate s=d (ln delta_m)/d (ln a) (a is the scale factor) that is larger than in the LCDM model. We find the window in the free parameter space of our models for which spatial dispersion of the growth index gamma_0= gamma(z=0) (z is the redshift) appears in the range of values 0.40< gamma_0<0.55, as well as the region in parameter space for which there is essentially no dispersion and gamma_0 converges to values around 0.40<gamma_0<0.43. These latter values are much lower than in the LCDM model. We show that these unusual dispersed or converged spectra are present in most of the viable f(R) models with m(z=0) larger than the order of 10^{-6}. These properties will be essential in the quest for f(R) modified gravity models using future high-precision observations and they confirm the possibility to distinguish clearly most of these models from the LCDM model.

The $f(R)$ gravity models proposed by Hu-Sawicki and Starobinsky are generic for local gravity constraints to be evaded. The large deviations from these models either result into violation of local gravity constraints or the modifications are not distinguishable from cosmological constant. The curvature singularity in these models is generic but can be avoided provided that proper fine tuning is imposed on the evolution of scalaron in the high curvature regime. In principle, the problem can be circumvented by incorporating quadratic curvature correction in the Lagrangian though it might be quite challenging to probe the relevant region numerically.

We consider the linear growth of matter perturbations on low redshifts in modified gravity Dark Energy (DE) models where G_eff(z,k) is explicitly scale-dependent. Dispersion in the growth today will only appear for scales of the order the critical scale ~ \lambda_{c,0}, the range of the fifth-force today. We generalize the constraint equation satisfied by the parameters \gamma_0(k) and \gamma'_0(k) \equiv \frac{d\gamma(z,k)}{dz}(z=0) to models with G_{eff,0}(k) \ne G. Measurement of \gamma_0(k) and \gamma'_0(k) on several scales can provide information about \lambda_{c,0}. In the absence of dispersion when \lambda_{c,0} is large compared to the probed scales, measurement of \gamma_0 and \gamma'_0 provides a consistency check independent of \lambda_{c,0}. This applies in particular to results obtained earlier for a viable f(R) model.

We consider the linear growth of matter perturbations on low redshifts in some $f(R)$ dark energy (DE) models. We discuss the definition of dark energy (DE) in these models and show the differences with scalar-tensor DE models. For the $f(R)$ model recently proposed by Starobinsky we show that the growth parameter $\gamma_0\equiv \gamma(z=0)$ takes the value $\gamma_0\simeq 0.4$ for $\Omega_{m,0}=0.32$ and $\gamma_0\simeq 0.43$ for $\Omega_{m,0}=0.23$, allowing for a clear distinction from $\Lambda$CDM. Though a scale-dependence appears in the growth of perturbations on higher redshifts, we find no dispersion for $\gamma(z)$ on low redshifts up to $z\sim 0.3$, $\gamma(z)$ is also quasi-linear in this interval. At redshift $z=0.5$, the dispersion is still small with $\Delta \gamma\simeq 0.01$. As for some scalar-tensor models, we find here too a large value for $\gamma'_0\equiv \frac{d\gamma}{dz}(z=0)$, $\gamma'_0\simeq -0.25$ for $\Omega_{m,0}=0.32$ and $\gamma'_0\simeq -0.18$ for $\Omega_{m,0}=0.23$. These values are largely outside the range found for DE models in General Relativity (GR). This clear signature provides a powerful constraint on these models.

We consider asymptotically stable scalar-tensor dark energy (DE) models for which the equation of state parameter $w_{DE}$ tends to zero in the past. The viable models are of the phantom type today, however this phantomness is milder than in General Relativity if we take into account the varying gravitational constant when dealing with the SNIa data. We study further the growth of matter perturbations and we find a scaling behaviour on large redshifts which could provide an important constraint. In particular the growth of matter perturbations on large redshifts in our scalar-tensor models is close to the standard behaviour $\delta_m \propto a$, while it is substantially different for the best-fit model in General Relativity for the same parametrization of the background expansion. As for the growth of matter perturbations on small redshifts, we show that in these models the parameter $\gamma'_0\equiv \gamma'(z=0)$ can take absolute values much larger than in models inside General Relativity. Assuming a constant $\gamma$ when $\gamma'_0$ is large would lead to a poor fit of the growth function $f$. This provides another characteristic discriminative signature for these models.

We consider the linear growth of matter perturbations in various dark energy (DE) models. We show the existence of a constraint valid at $z=0$ between the background and dark energy parameters and the matter perturbations growth parameters. For $\Lambda$CDM $\gamma'_0\equiv \frac{d\gamma}{dz}_0$ lies in a very narrow interval $-0.0195 \le \gamma'_0 \le -0.0157$ for $0.2 \le \Omega_{m,0}\le 0.35$. Models with a constant equation of state inside General Relativity (GR) are characterized by a quasi-constant $\gamma'_0$, for $\Omega_{m,0}=0.3$ for example we have $\gamma'_0\approx -0.02$ while $\gamma_0$ can have a nonnegligible variation. A smoothly varying equation of state inside GR does not produce either $|\gamma'_0|>0.02$. A measurement of $\gamma(z)$ on small redshifts could help discriminate between various DE models even if their $\gamma_0$ is close, a possibility interesting for DE models outside GR for which a significant $\gamma'_0$ can be obtained.

We derive the conditions under which dark energy models whose Lagrangian densities f are written in terms of the Ricci scalar R are cosmologically viable. We show that the cosmological behavior of f(R) models can be understood by a geometrical approach consisting in studying the m(r) curve on the (r, m) plane, where m=Rf_{,RR}/f_{,R} and r=-Rf_{,R}/f with f_{,R}=df/dR. This allows us to classify the f(R) models into four general classes, depending on the existence of a standard matter epoch and on the final accelerated stage. The existence of a viable matter dominated epoch prior to a late-time acceleration requires that the variable m satisfies the conditions m(r) approx+0 and dm/dr>-1 at r approx-1. For the existence of a viable late-time acceleration we require instead either (i) m=-r-1, (sqrt{3}-1)/2<m<1 and dm/dr<-1 or (ii) 0<m<1 at r=-2. These conditions identify two regions in the (r, m) space, one for the matter era and the other for the acceleration. Only models with a m(r) curve that connects these regions and satisfy the requirements above lead to an acceptable cosmology. The models of the type f(R)=alpha R^{-n} and f=R+alpha R^{-n} do not satisfy these conditions for any n>0 and n<-1 and are thus cosmologically unacceptable. Similar conclusions can be reached for many other examples discussed in the text. In most cases the standard matter era is replaced by a cosmic expansion with scale factor a propto t^{1/2}. We also find that f(R) models can have a strongly phantom attractor but in this case there is no acceptable matter era.

We present here some recent results concerning scalar-tensor Dark Energy models. These models are very interesting in many respects: they allow for a consistent phantom phase, the growth of matter perturbations is modified. Using a systematic expansion of the theory at low redshifts, we relate the possibility to have phantom like DE to solar system constraints.

We consider the viability of dark energy (DE) models in the framework of the scalar-tensor theory of gravity, including the possibility to have a phantom DE at small redshifts $z$ as admitted by supernova luminosity-distance data. For small $z$, the generic solution for these models is constructed in the form of a power series in $z$ without any approximation. Necessary constraints for DE to be phantom today and to cross the phantom divide line $p=-\rho$ at small $z$ are presented. Considering the Solar System constraints, we find for the post-Newtonian parameters that $\gamma_{PN}<1$ and $\gamma_{PN,0}\approx 1$ for the model to be viable, and $\beta_{PN,0}>1$ (but very close to 1) if the model has a significantly phantom DE today. However, prospects to establish the phantom behaviour of DE are much better with cosmological data than with Solar System experiments. Earlier obtained results for a $\Lambda$-dominated universe with the vanishing scalar field potential are extended to a more general DE equation of state confirming that the cosmological evolution of these models rule them out. Models of currently fantom DE which are viable for small $z$ can be easily constructed with a constant potential; however, they generically become singular at some higher $z$. With a growing potential, viable models exist up to an arbitrary high redshift.