In this paper, two inhomogeneous single fluid models for the Universe, which are able to naturally solve the $H_{0}$ tension problem, are discussed. The analysis is based on a Bayesian Machine Learning approach that uses a generative process. The adopted method allows to constrain the free parameters of each model by using the model itself, only. The observable is taken to be the Hubble parameter, obtained from the generative process. Using the full advantages of our method, the models are constrained for two redshift ranges. Namely, first this is done with mock $H(z)$ data over $z\in [0,2.5]$, thus covering known $H(z)$ observational data, which are most helpful to validate the fit results. Then, aiming to extend to redshift ranges to be covered by the most recent ongoing and future planned missions, the models are constrained for the range $z\in[0,5]$, too. Full validation of the results for this extended redshift range will have to wait for the near future, when higher redshift $H(z)$ data become available. This makes our models fully falsifiable. Finally, our second model here is able to explain the BOSS reported value for $H(z)$ at $z=2.34$.

In this work, we revisit the dynamics of pre-inflationary universe with a family of $\alpha-$attractor potentials, in the framework of loop quantum cosmology, in which the big bang singularity is generically resolved purely with quantum geometric effects, and replaced by a quantum bounce. At the bounce, the background evolution is divided into two distinct classes, the first is dominated by the kinetic energy of the inflaton field and the second by the potential energy. In both classes, we find the physically viable initial conditions numerically that provide not only the slow-roll inflation, but also sufficient e-folds to be compatible with observations. In the entire range of kinetic energy dominated initial conditions (except some subsets of Models 2 and 4), the background evolution prior to reheating is always split into three different phases: bouncing, transition and slow-roll inflation. In the bouncing phase, the numerical evolution of the scale factor is independent not only of the initial data, but also the inflationary potentials, as long as it is dominated by the kinetic energy, and can be well approximated by an analytical solution, whereas in the potential energy dominated case, such approximated results do not exist. Moreover, we study the phase space analysis for a class of $\alpha-$attractor potentials, and discuss the phase space trajectories for physically viable initial conditions of the inflaton field.

We show that there is a correspondence between cosmology from non-extensive thermodynamics and cosmology with fluids of redefined and generalized equation of state. We first establish the correspondence in the case of basic non-extensive thermodynamics, and then we proceed by investigating the more consistent case, from the quantum field theoretical point of view, of varying exponent, namely depending on the scale. The obtained duality provides a way of explaining the complicated phenomenological forms of the effective fluid equation-of-state parameters that are being broadly used in the literature, since their microphysical origin may indeed lie in the non-extensive thermodynamics of spacetime. Finally, concerning the cosmological behavior, we show that at late times the effective fluid may drive the universe acceleration even in the absence of an explicit cosmological constant, and even if the initial fluid is the standard dust matter one. Similarly, at early times we obtain an effective cosmological constant which is enhanced through screening, and hence it can drive a a successful inflation without spoiling the correct late-time acceleration.

In this paper, we investigate a scenario in which late time cosmic acceleration might arise due to coupling between dark matter and baryonic matter without resorting to dark energy or large scale modification of gravity associated with extra degrees of freedom. The scenario can give rise to late time acceleration in Jordan frame and no acceleration in Einstein frame - \textit{generic modification of gravity} caused by disformal coupling. Using a simple parametrization of the coupling function, in maximally disformal case, we constrain the model parameters by using the age constraints due to globular cluster data. We also obtain observational constraints on the parameters using $H(z)+SNIa+BAO$ data sets. In this case, we distinguish between phantom and non phantom acceleration and show that the model can give rise to phantom behavior in a narrow region of parameter space.

In order to apply holography and entropy relations to the whole universe, which is a gravitational and thus nonextensive system, for consistency one should use the generalized definition for the universe horizon entropy, namely Tsallis nonextensive entropy. We formulate Tsallis holographic dark energy, which is a generalization of standard holographic dark energy quantified by a new dimensionless parameter $\delta$, possessing the latter as a particular sub-case. We provide a simple differential equation for the dark energy density parameter, as well as an analytical expression for its equation-of-state parameter. In this scenario the universe exhibits the usual thermal history, namely the successive sequence of matter and dark-energy epochs, before resulting in a complete dark energy domination in the far future. Additionally, the dark energy equation-of-state parameter presents a rich behavior and, according to the value of $\delta$, it can be quintessence-like, phantom-like, or experience the phantom-divide crossing before or after the present time. Finally, we confront the scenario with Supernovae type Ia and Hubble parameter observational data, and we show that the agreement is very good, with $\delta$ preferring a value slightly larger than its standard value 1.

We investigate the cosmological applications of fluids having an equation of state which is the analog to the one related to the isotropic deformation of crystalline solids, that is containing logarithmic terms of the energy density, allowing additionally for a bulk viscosity. We consider two classes of scenarios and we show that they are both capable of triggering the transition from deceleration to acceleration at late times. Furthermore, we confront the scenarios with data from Supernovae type Ia (SN Ia) and Hubble function observations, showing that the agreement is excellent. Moreover, we perform a dynamical system analysis and we show that there exist asymptotic accelerating attractors, arisen from the logarithmic terms as well as from the viscosity, which in most cases correspond to a phantom late-time evolution. Finally, for some parameter regions we obtain a nearly de Sitter late-time attractor, which is a significant capability of the scenario since the dark energy, although dynamical, stabilizes at the cosmological constant value.

We study the bounce and cyclicity realization in the framework of new gravitational scalar-tensor theories. In these theories the Lagrangian contains the Ricci scalar and its first and second derivatives, in a specific combination that makes them free of ghosts, and transformed into the Einstein frame they are proved to be a subclass of bi-scalar extensions of general relativity. We present analytical expressions for the bounce requirements, and we examine the necessary qualitative behavior of the involved functions that can give rise to a given scale factor. Having in mind these qualitative forms, we reverse the procedure and we construct suitable simple Lagrangian functions that can give rise to a bounce or cyclic scale factor.

We study the generalized $\alpha$ attractor model in the context of the late time cosmic acceleration. The model interpolates between the scaling freezing and thawing dark energy models. In the slow roll region, the original potential is modified whereas the modification ceases in the asymptotic region and the effective potential behaves as the quadratic one. In our setting, the field rolls slowly around the present epoch and mimics the scaling behavior in the future. We obtain observational constraints on the model parameters by using an integrated data base (SN+Hubble+BAO+CMB).

We study relic gravitational waves in the paradigm of quintessential inflation. In this framework, irrespective of the underlying model, inflation is followed by the kinetic regime. Thereafter, the field energy density remains sub-dominant before the onset of acceleration. We carry out model independent analysis to obtain the temperature at the end of inflation and the estimate for upper bound on the Hubble parameter to circumvent the problem due to relic gravitational waves. In this process, we used Planck 2015 data to constrain the inflationary phase. We demonstrate that the required temperature can be produced by the mechanism of instant preheating. The generic feature of the scenario includes the presence of kinetic regime after inflation which results into blue spectrum of gravitational wave background at high frequencies. We discuss the prospects of detection of relic gravitational wave background in the advanced LIGO and LISA space-born gravitational wave missions. Finally we consider a concrete model to realize the paradigm of quintessential inflation and show that inflationary as well as post-inflationary evolution can successfully be described by the inflaton potential, $V(\phi) \propto Exp(-\lambda \phi^n/\Mpl^n)(n>1)$, by suitably constraining the parameters of the model.

In this paper, we will analyse the cosmological models in Einstein-aether gravity, which is a modified theory of gravity in which a time-like vector field breaks the Lorentz symmetry. We will use this formalism to analyse different cosmological models with different behavior of the scale factor. In this analysis, we will use a certain functional dependence of the dark energy on the Hubble parameter. It will be demonstrated that the aether vector field has a non-trivial effect on these cosmological models. We will also perform the \emph{Om} diagnostic in Einstein-aether gravity. Thus, we will fit parameters of the cosmological models using recent observational data.

In this paper, we consider rolling tachyon, with steep run-away type of potentials non-minimally coupled to massive neutrino matter. The coupling dynamically builds up at late times as neutrino matter turns non-relativistic. In case of scaling and string inspired potentials, we have shown that non-minimal coupling leads to minimum in the field potential. Given a suitable choice of model parameters, it is shown to give rise to late-time acceleration with the desired equation of state.

In this paper, we study the Power Law Entropy Corrected Holographic Dark Energy (PLECHDE) model in the framework of a non-flat Universe and of Ho\v{r}ava-Lifshitz cosmology with infrared cut-off given by recently proposed Granda-Oliveros cut-off, which contains one term proportional to the Hubble parameter squared $H^2$ and one term proportional to the first time derivative of the Hubble parameter $\dot{H}$. Moreover, this cut-off is characterized by two constant parameters, $\alpha$ and $\beta$. For the two cases corresponding to non-interacting and interacting DE and Dark Matter (DM), we derive the evolutionary form of the energy density of DE $\Omega_D'$, the Equation of State (EoS) parameter of DE $\omega_D$ and the deceleration parameter $q$. Using the parametrization of the EoS parameter $\omega_D\left(z\right)=\omega_0+\omega_1 z$, we obtain the expressions of the two parameters $\omega_0$ and $\omega_1$. We also study the statefinder parameters $\left\{ r,s \right\}$, the snap and lerk cosmographic parameters and the squared speed of the sound $v_s^2$. We also calculate the values of the quantities we study for different values of the running parameter $\lambda$ and for different set of values of $\alpha$ and $\beta$.

In the present work, we study the cosmological model with fermionic field and with the non-canonical kinetic term (f-essence). We also present some important reductions of the model as well as its some generalizations. We also found the exact solution of the model and examine the influence of such gravity-fermion interaction on the observed accelerated expansion of our universe. Several fermionic and bosonic-fermionic DBI models were constructed as some particular reductions of f-essence and g-essence, respectively and their Chaplygin gas counterparts are found. Finally, some nonlinear models of bouncing and cyclic universes are proposed.

We consider generic models of quintessence and we investigate the influence of massive neutrino matter with field-dependent masses on the matter power spectrum. In case of minimally coupled neutrino matter, we examine the effect in tracker models with inverse power-law and double exponential potentials. We present detailed investigations for the scaling field with a steep exponential potential, non-minimally coupled to massive neutrino matter, and we derive constraints on field-dependent neutrino masses from the observational data.

In literature there is a model of modified gravity in which the matter Lagrangian is coupled to the geometry via trace of the stress-energy momentum tensor $T=T_{\mu}^{\mu}$. This type of modified gravity is called as $f(R,T)$ in which $R$ is Ricci scalar $R=R_{\mu}^{\mu}$. We extend manifestly this model to include the higher derivative term $\Box R$. We derived equation of motion (EOM) for the model by starting from the basic variational principle. Later we investigate FLRW cosmology for our model. We show that de Sitter solution is unstable for a generic type of $f(R,\Box R, T)$ model. Furthermore we investigate an inflationary scenario based on this model. A graceful exit from inflation is guaranteed in this type of modified gravity.

We briefly review the problems and prospects of the standard lore of dark energy and notice that the cosmological constant problems are often over emphasized. In order to keep the discussion pedagogical aimed at a wider audience, we have avoided technical complications in several places and resorted to heuristic arguments based on physical perceptions. We presented underlying ideas of modified theories based upon chameleon mechanism and Vainshtein screening. We have given a lucid illustration of recently investigated ghost free non linear massive gravity. Again we have sacrificed rigor and confined to the basic ideas that led to the formulation of $dRGT$. The review ends with a brief discussion on the difficulties of the theory applied to cosmology with an optimism that a consistent solution of the problem is round the corner.

We investigate two classes of models of quintessential inflation, based upon canonical as well as noncanonical scalar fields. In particular, introducing potentials steeper than the standard exponential, we construct models that can give rise to a successful inflationary phase, with signatures consistent with Planck 2015 results. Additionally, using nonminimal coupling of the scalar field with massive neutrino matter, we obtain the standard thermal history of the Universe, with late-time cosmic acceleration as the last stage of evolution. In both cases, inflation and late-time acceleration are connected by a tracker solution.

We study Galileon scalar field model by considering the lowest order Galileon term in the lagrangian, $(\partial_{\mu} \phi)^2 \Box\phi$ by invoking a field potential. We use Statefinder hierarchy to distinguish the light mass galileon models with different potentials amongst themselves and from the $\Lambda$CDM behaviour. The $Om$ diagnostic is applied to cosmological dynamics and observational constraints on the model parameters are studied using SN+Hubble+BAO data.

We examine the Kaluza-Klein (KK) dimensional reduction from higher-dimensional space-time and the properties of the resultant Bergmann-Wagoner general action of scalar-tensor theories. With the analysis of the perturbations, we also investigate the stability of the anti-de Sitter (AdS) space-time in the $D\in\mathcal{N}$-dimensional Einstein gravity with the negative cosmological constant. Furthermore, we derive the conditions for the dimensional reduction to successfully be executed and present the KK compactification mechanism.

This paper is devoted to study the possibility of forming anisotropic compact stars in modified Gauss-Bonnet, namely called as $f(G)$ theory of gravity which is one of the strong candidates,responsible for the accelerated expansion of the universe. For this purpose, we have used analytical solution of Krori and Barua metric to the Einstein field equations with anisotropic form of matter and power law model of $f(G)$ gravity. To determine the unknown constants in Krori and Barua metric, we have used the sample of compact stars, 4$U$1820-30, Her X-1, SAX J 1808-3658. The physical behavior of these stars have been analyzed with the observational data. In this setting, we have checked all the regularity conditions and stability of the compact stars 4$U$1820-30, Her X-1, SAX J1808-3658.

We consider higher-derivative quantum gravity where renormalization group improved effective action beyond one-loop approximation is derived. Using this effective action, the quantum-corrected FRW equations are analyzed. De Sitter universe solution is found. It is demonstrated that such de Sitter inflationary universe is instable. The slow-roll inflationary parameters are calculated. The contribution of renormalization group improved Gauss-Bonnet term to quantum-corrected FRW equations as well as to instability of de Sitter universe is estimated. It is demonstrated that in this case the spectral index and tensor-to-scalar ratio are consistent with Planck data.

We investigate the turbulence effect in dark fluid universe with linear inhomogeneous equation of state. Attention is attached to two physical situations. First, we perform the perturbative analysis of turbulence and check its effects around the Big Rip. Later, treating the turbulence energy density as a part of total dark fluid, we study the stability of the system. The result shows that the stability is achieving as the energy density of turbulence decreases, changing into heat (the radiation), in perfect agreement with the avoidance of the Big Rip.

In this paper, we examine observational constraints on the power law cosmology; essentially dependent on two parameters $H_0$ (Hubble constant) and $q$ (deceleration parameter). We investigate the constraints on these parameters using the latest 28 points of H(z) data and 580 points of Union2.1 compilation data and, compare the results with the results of $\Lambda$CDM. We also forecast constraints using a simulated data set for the future JDEM, supernovae survey. Our studies give better insight into power law cosmology than the earlier done analysis by Kumar [arXiv:1109.6924] indicating it tuning well with Union2.1 compilation data but not with H(z) data. However, the constraints obtained on $<H_0>$ and $<q>$ i.e. $H_0$ average and $q$ average using the simulated data set for the future JDEM, supernovae survey are found to be inconsistent with the values obtained from the H(z) and Union2.1 compilation data. We also perform the statefinder analysis and find that the power-law cosmological models approach the standard $\Lambda$CDM model as $q\rightarrow -1$. Finally, we observe that although the power law cosmology explains several prominent features of evolution of the Universe, it fails in details.

This pedagogical review is devoted to quintessential inflation, which refers to unification of inflation and dark energy using a single scalar field. We present a brief but concise description of the concepts needed to join the two ends, which include discussion on scalar field dynamic, conformal coupling, instant preheating and relic gravitational waves. Models of quintessential inflation broadly fall into two classes, depending upon the early and late time behavior of the field potential. In the first type we include models in which the field potential is steep for most of the history of the Universe but turn shallow at late times, whereas in the second type the potential is shallow at early times followed by a steep behavior thereafter. In models of the first category inflation can be realized by invoking high-energy brane-induced damping, which is needed to facilitate slow roll along a steep potential. In models of second type one may invoke a non-minimal coupling of the scalar field with massive neutrino matter, which might induce a minimum in the potential at late times as neutrinos turn non-relativistic. In this category we review a class of models with non-canonical kinetic term in the Lagrangian, which can comply with recent B mode polarization measurements. The scenario under consideration is distinguished by the presence of a kinetic phase, which precedes the radiative regime, giving rise to blue spectrum of gravity waves generated during inflation. We highlight the generic features of quintessential inflation and also discuss on issues related to Lyth bound.

We consider the Noether Symmetry Approach for a cosmological model derived from a tachyon scalar field $T$ with a Dirac-Born-Infeld Lagrangian and a potential $V(T)$. Furthermore, we assume a coupled canonical scalar field $\phi$ with an arbitrary interaction potential $B(T,\phi)$. Exact solutions are derived consistent with the accelerated behavior of cosmic fluid.

Based on a stringy inspired Gauss-Bonnet (GB) modification of classical gravity, we constructed a model for neutron stars. We derived the modified forms of Tolman-Oppenheimer-Volkoff (TOV) equations for a generic function of $f(G)$ gravity. The hydrostatic equations remained unchanged but the dynamical equations for metric functions are modified due to the effects of GB term.

We study inflation induced by (power-low) scalar curvature corrections to General Relativity. The class of inflationary scalar potentials $V(\sigma)\sim\exp[n\,\sigma]$, $n$ general parameter, is investigated in the Einsein frame and the corresponding actions in the Jordan frame are derived. We found the conditions for which these potentials are able to reproduce viable inflation according with the last cosmological data and lead to large scalar curvature corrections which emerge only at a mass scale larger than the Planck mass. Cosmological constant may appear or be set equal to zero in the Jordan frame action without changing the behaviour of the model during inflation. Moreover, polynomial corrections to General Relativity are analyzed in detail. When de Sitter space-time emerges as an exact solution of the models, it is necessary to use perturbative equations in the Jordan framework to study their dynamics during the inflation. In this case, we demonstrate that the Ricci scalar decreases after a correct amount of inflation, making the models consistent with the observable evolution of the universe.

Quintessential inflation refers to an attempt to unify inflation and late-time cosmic acceleration using a single scalar field. In this letter we consider two different classes of quintessential inflation, one of which is based upon a Lagrangian with non-canonical kinetic term $k^2(\phi)\partial^\mu \phi \partial_\mu \phi$ and a steep exponential potential while the second class uses the concept of steep brane world inflation. We show that in both cases the Lyth bound can be evaded, despite the large tensor-to-scalar ratio of perturbations. The post-inflationary dynamics is consistent with nucleosynthesis constraint in these cases.

A modified Mimetic gravity (MMG) is proposed as a generalization of general relativity. The model contain a physical metric which is function of an auxiliary (unphysical) metric and a Lyra's metric. We construct different kinds of conformally invariant models in different levels of the expansion parameter $\lambda$. This model phenomenologically has been extended to higher order forms. Cosmology of a certain class of such models has been investigated in details. A cosmological solution has been proposed in inhomogeneous form of scalar field. For homogenous case, energy conditions are widely investigated. We have shown that the system evaluated at intervals shorter than a certain time $T_c$ meets all the energy conditions.

As an extension of Dvali-Gabadadze-Porrati (DGP) model, the Galileon theory has been proposed to explain the "self-accelerating problem" and "ghost instability problem". In this Paper, we extend the Galileon theory by considering a non-minimally coupled Galileon scalar with gravity. The statefinder analysis, $Om(z)$ diagnostic and constraint the model parameters have been investigate , we find $\Omega_\text{m0}=0.263_{-0.031}^{+0.031}$, $n=1.53_{-0.37}^{+0.21}$ (at the $95\%$ confidence level) with $\chi^2_{\text{min}}=473.376$. Further we show that due to the SNe Ia + BAO data ,our model behaves like a phantom-like dark energy.

We explore conformal-anomaly driven inflation in $F(R)$ gravity without invoking the scalar-tensor representation. We derive the stress-energy tensor of the quantum anomaly in the flat homogeneous and isotropic universe. We investigate a suitable toy model of exponential gravity plus the quantum contribution due to the conformal anomaly, which leads to the de Sitter solution. It is shown that in $F(R)$ gravity model, the curvature perturbations with its enough amplitude consistent with the observations are generated during inflation. We also evaluate the number of $e$-folds at the inflationary stage and the spectral index $n_\mathrm{s}$ of scalar modes of the curvature perturbations by analogy with scalar tensor theories, and compare them with the observational data. As a result, it is found that the Ricci scalar decreases during inflation and the standard evolution history of the universe is recovered at the small curvature regime. Furthermore, it is demonstrated that in our model, the tensor-to-scalar ratio of the curvature perturbations can be a finite value within the $68\%\,\mathrm{CL}$ error of the very recent result found by the BICEP2 experiment.

In this paper, we investigate a scenario of variable gravity and apply it to the unified description of inflation and late time cosmic acceleration dubbed quintessential inflation. The scalar field called "cosmon" which in this model unifies both the concepts reduces to inflaton at early epochs. We calculate the slow-roll parameters, the Hubble parameter at the end of inflation, the reheating temperature,the tensor-to-scalar ration and demonstrate the agreement of the model with observations and the Planck data.As for the post inflationary dynamics, cosmon tracks the background before it exits the scaling regime at late times. The scenario gives rise to correct epoch sequence of standard cosmology, namely, radiative regime, matter phase and dark-energy. We show that the long kinetic regime after inflation gives rise to enhancement of relic gravity wave amplitude resulting into violation of nucleosynthesis constraint at the commencement of radiative regime in case of an inefficient reheating mechanism such as gravitational particle production. Instant preheating is implemented to successfully circumvent the problem. As a generic feature, the scenario gives rise to a blue spectrum for gravity waves on scales smaller than the comoving horizon scale at the commencement of the radiative regime.

In this paper, we focus on general features of quintessential inflation which is an effort to unify inflation and dark energy using a single scalar field. We describe a class of models of quintessential inflation which can give rise to the tensor to scalar ratio of perturbations consistent with recent BICEP2 measurements. The scale of inflation in the model is around the GUT scale and there is large parameter space consistent with the recent findings.

In this paper we investigate the case of non-minimal coupling in the (extended) nonlinear massive gravity theories. We first consider massive gravity in the Brans-Dicke background such that the graviton mass is replaced by $A^2(\sigma)m$ where $\sigma$ is the Brans-Dicke field and $A(\sigma)$ is conformal coupling and show that there is no viable thermal history of the universe in this case. We then invoke a cubic galileon term as nonlinear completion of the $\sigma$ Lagrangian and show that there is a stable de Sitter solution in this case. However, the de Sitter is blocked by the matter phase which is also a simultaneous attractor of the dynamics. The de Sitter phase can, however, be realized by invoking unnatural fine tunings. We next investigate cosmology of quasi-dilaton nonlinear massive gravity with non-minimal coupling. As a generic feature of the non-minimal coupling, we show that the model exhibits a transient phantom phase which is otherwise impossible. While performing the observational data analysis on the models, we find that a small value of coupling constant is allowed for quasi-dilaton nonlinear massive gravity. For both the cases under consideration, it is observed that we have an effective pressure of matter which comes from the constraint equation. For mass-varying nonlinear massive gravity in the Brans-Dicke background, the effective pressure of matter is non zero which affects the evolution of the Hubble parameter thereby spoiling consistency of the model with data. As for, quasi-dilaton nonlinear massive gravity, the effective pressure of matter can be kept around zero by controlling the coupling constant, the model is shown to be fit well with observations.

In this paper, we consider a simple form of expansion history of Universe referred to as the hybrid expansion law - a product of power-law and exponential type of functions. The ansatz by construction mimics the power-law and de Sitter cosmologies as special cases but also provides an elegant description of the transition from deceleration to cosmic acceleration. We point out the Brans-Dicke realization of the cosmic history under consideration. We construct potentials for quintessence, phantom and tachyon fields, which can give rise to the hybrid expansion law in general relativity. We investigate observational constraints on the model with hybrid expansion law applied to late time acceleration as well as to early universe a la nucleosynthesis.

The presence of spinor fields is considered in the framework of some extensions of teleparallel gravity, where the Weitzenb\"ock connection is assumed. Some well known models as the Chaplygin gas and its generalizations are reconstructed in terms of a spinor field in the framework of teleparallel gravity. In addition, the {\Lambda}CDM model is also realized with the presence of a spinor field where a simple self-intercating term is considered and the corresponding action is reconstructed. Other cosmological solutions and the reconstruction of the gravitational action in terms of the scalar torsion is studied.

We investigate the future evolution of Nariai black hole which is extremal limit of Schwarzschild-de Sitter one in modified gravity. The perturbations equations around Nariai black hole are derived in static and cosmological patches for general $F(R)$-gravity. The analytical and numerical study of several realistic $F(R)$-models shows the occurence of rich variety of scenarios: instabilities, celebrated Hawking evaporation and anti-evaporation of black hole. The realization of specific scenario depends on the model under consideration. It is remarkable that the presence of such primordial black holes at current universe may indicate towards the modified gravity which supports the anti-evaporation as preferrable model. As some generalization we extend the study of Nariai black hole evolution to modified Gauss-Bonnet gravity. The corresponding perturbations equations turn out to be much more complicated than in the case of $F(R)$ gravity. For specific example of modified Gauss-Bonnet gravity we demonstrate that Nariai solution maybe stable.

In this paper we investigate charged static black holes in 4D for generalized teleparallel models of gravity, based on torsion as the geometric object for describing gravity according to the equivalence principle. As a motivated idea, we introduce a set of non-diagonal tetrads and derive the full system of non linear differential equations. We prove that the common Schwarzschild gauge is applicable only when we study linear f(T) case. We reobtain the Reissner-Nordstrom-de Sitter (or RN-AdS) solution for the linear case of f(T) and perform a parametric cosmological reconstruction for two nonlinear models. We also study in detail a type of the no-go theorem in the framework of this modified teleparallel gravity.

In literature usual point like symmetries of the Lagrangian have been introduced to study the symmetries and the structure of the fields. This kind of Noether symmetry is a subclass of a more general family of symmetries, called Noether Gauge Symmetries (NGS). Motivated by this mathematical tool, in this article, we discuss the generalized Noether symmetry of Quintom model of dark energy, which is a two component fluid model of quintessence and phantom fields. Our model is a generalization of the Noether symmetries of a single and multiple components which have been investigated in detail before. We found the general form of the quintom potential in which the whole dynamical system has a point like symmetry. We investigated different possible solutions of the system for diverse family of gauge function. Specially, we discovered two family of potentials, one corresponds to a free quintessence (phantom) and the second is in the form of quadratic interaction between two components. These two families of potential functions are proposed from the symmetry point of view, but in the quintom models they are used as phenomenological models without clear mathematical justification. From integrability point of view, we found two forms of the scale factor: one is power law and second is de-Sitter. Some cosmological implications of the solutions have been investigated

Motivated by earlier works on reconstruction of modified gravity models with dark energy components, we extend them by considering a newly proposed model of (m, n)- type of holographic dark energy for two models of modified gravity, f(R) and f(T) theories, where R and T represent Ricci scalar and torsion scalar respectively. Specifically we reconstruct the two later gravity models and discuss their viability and cosmography. The obtained gravity models are ghost free, compatible with local solar system tests and describe effective positive gravitational constant.

A generic feature of viable exponential $F(R)$-gravity is investigated. An additional modification to stabilize the effective dark energy oscillations during matter era is proposed and applied to two viable models. An analysis on the future evolution of the universe is performed. Furthermore, a unified model for early and late-time acceleration is proposed and studied.

This paper is devoted to the study of Noether gauge symmetries of $f(T)$ gravity minimally coupled with a canonical scalar field. We explicitly determine the unknown functions of the theory $f(T),V(\phi), W(\phi)$. We have shown that there are two invariants for this model, one of which defines the Hamiltonian $H$ under time invariance (energy conservation) and the other is related to scaling invariance. We show that the equation of state parameter in the present model can cross the cosmological constant boundary. The behavior of Hubble parameter in our model closely matches to that of $\Lambda$CDM model, thus our model is an alternative to the later.

Cosmological reconstruction of Little Rip model in $f(R, T)$ gravity is investigated, where $R$ is the curvature scalar and $T$ the trace of the energy momentum tensor. The model perfectly reproduces the present stage of the universe, characterized by the $\Lambda CDM$ model, without singularity at future finite-time (without the Big Rip). The input parameters are determined according to Supernovae Cosmology data and perfectly fit with the WMAP around the Little Rip. Moreover, the thermodynamics is considered in this Little Rip cosmology and it is illustrated that the second law of thermodynamics is always satisfied around the Little Rip universe for the temperature inside the horizon being the same as that of the apparent horizon. Moreover, we show the existence of a stable fixed point in the Little Rip universe which confirms that this is actually a late-time attractor in the phantom-dominated universe. The linear perturbation analysis is performed around the critical points, showing that the Little Rip model obtained is stable.

We study the generalized second law (GSL) of thermodynamics in $f(T)$ cosmology. We consider the universe as a closed bounded system filled with $n$ component fluids in the thermal equilibrium with the cosmological boundary. We use two different cosmic horizons: the future event horizon and the apparent horizon. We show the conditions under which the GSL will be valid in specific scenarios of the quintessence and the phantom energy dominated eras. Further we associate two different entropies with the cosmological horizons: with a logarithmic correction term and a power-law correction term. We also find the conditions for the GSL to be satisfied or violated by imposing constraints on model parameters.

In this paper, we consider a theory of gravity with a metric-dependent torsion namely the $F(R,T)$ gravity, where $R$ is the curvature scalar and $T$ is the torsion scalar. We study a geometric root of such theory. In particular we give the derivation of the model from the geometrical point of view. Then we present the more general form of $F(R,T)$ gravity with two arbitrary functions and give some of its particular cases. In particular, the usual $F(R)$ and $F(T)$ gravity theories are the particular cases of the $F(R,T)$ gravity. In the cosmological context, we find that our new gravitational theory can describes the accelerated expansion of the universe.

A generic feature of viable $F(R)$ gravity is investigated: It is demonstrated that during the matter dominated era the large frequency oscillations of the effective dark energy may influence the behavior of higher derivatives of the Hubble parameter with the risk to produce some singular unphysical solutions at high redshift. This behavior is explicitly analyzed for realistic $F(R)$ models, in particular, exponential gravity and a power form model. To stabilize such oscillations, we consider the additional modification of the models via a correction term which does not destroy the viability properties. A detailed analysis on the future evolution of the universe and the evolution history of the growth index of the matter density perturbations are performed. Furthermore, we explore two applications of exponential gravity to the inflationary scenario. We show how it is possible to obtain different numbers of $e$-folds during the early-time acceleration by making different choices of the model parameters in the presence of ultrarelativistic matter, which destabilizes inflation and eventually leads to the exit from the inflationary stage. We execute the numerical analysis of inflation in two viable exponential gravity models. It is proved that at the end of the inflation, the effective energy density and curvature of the universe decrease and thus a unified description between inflation and the $\Lambda$CDM-like dark energy dominated era can be realized.

In this paper, we investigate the Noether symmetries of $F(T)$ cosmology involving matter and dark energy. In this model, the dark energy is represented by a canonical scalar field with a potential. Two special cases for dark energy are considered including phantom energy and quintessence. We obtain $F(T)\sim T^{3/4},$ and the scalar potential $V(\phi)\sim\phi^2$ for both models of dark energy and discuss quantum picture of this model. Some astrophysical implications are also discussed.

The Ekpyrotic scenario is studied in the context of some extensions of Ho\v{r}ava-Lifshitz gravity. Some particular solutions that lead to cyclic Hubble parameters are analyzed, where the corresponding gravitational actions are reconstructed by using several techniques and auxiliary fields. Its comparison with standard $F(R)$ gravity is performed. In addition, the so-called Little Rip, a stage of the universe evolution when some bounded systems may be dissolute, is also studied in this frame of theories.

In this paper, we have considered the g-essence and its particular cases, k-essence and f-essence, within the framework of the Einstein-Cartan theory. We have shown that a single fermionic field can give rise to the accelerated expansion within the Einstein-Cartan theory. The exact analytical solution of the Einstein-Cartan-Dirac equations is found. This solution describes the accelerated expansion of the Universe with the equation of state parameter $w=-1$ as in the case of $\Lambda$CDM model.

It is shown that the acceleration of the universe can be understood by considering a F(T) gravity models. For these F(T) gravity models, a variant of the accelerating cosmology reconstruction program is developed. Some explicit examples of F(T) are reconstructed from the background FRW expansion history.

We study the two particular models of g-essence with Yukawa type interactions between a scalar field $\phi$ and a classical Dirac field $\psi$. For the homogeneous, isotropic and flat Friedmann-Robertson-Walker universe filled with the such g-essence, some exact solutions of these models are found. Moreover, we reconstruct the corresponding scalar and fermionic potentials.

We study the g-essence model with Yukawa interactions between a scalar field $\phi$ and a Dirac field $\psi$. For the homogeneous, isotropic and flat Friedmann-Robertson-Walker universe filled with the such g-essence, the exact solution of the model is found. Moreover, we reconstruct the corresponding scalar and fermionic potentials which describe the coupled dynamics of the scalar and fermionic fields. It is shown that some particular g-essence models with Yukawa interactions correspond to the usual and generalized Chaplygin gas unified models of dark energy and dark matter. Also we present some scalar-fermionic Dirac-Born-Infeld models corresponding g-essence models with Yukawa interactions which again describe the unified dark energy-dark matter system.

In this paper, we consider three types of k-essence. These k-essence models were presented in the parametric forms. The exact analytical solutions of the corresponding equations of motion are found. It is shown that these k-essence models for the presented solutions can give rise to cosmic acceleration.

In the context of the so-called Gauss-Bonnet gravity, where the gravitational action includes function of the Gauss-Bonnet invariant, we study cosmological solutions, especially the well-known $\Lambda$CDM model. It is shown that the dark energy contribution and even the inflationary epoch can be explained in the frame of this kind of theories with no need of any other kind of component. Other cosmological solutions are constructed and the rich properties that this kind of theories provide are explored.

In the present work we analyze the g-essence model for the particular Lagrangian: $L=R+2[\alpha X^n+\epsilon Y-V(\psi,\bar{\psi})]$. The g-essence models were proposed recently as an alternative and as a generalization to the scalar k-essence. We have presented the 3 types solutions of the g-essence model. We reconstructed the corresponding potentials and the dynamics of the scalar and fermionic fields according the evolution of the scale factor. The obtained results shows that the g-essence model can describes the decelerated and accelerated expansion phases of the universe.

Modified teleparallel gravity theory with the torsion scalar have recently gained a lot of attention as a possible explanation of dark energy. We perform a thorough reconstruction analysis on the so-called $F(T)$ models, where $F(T)$ is some general function of the torsion term, and derive conditions for the equivalence between of $F(T)$ models with purely kinetic k-essence. We present a new class models of $F(T)$ - gravity and k-essence.

A self-similar motion of the generalized Chaplygin gas in its own gravitational field is considered. The problem is being studied numerically and analytically (in the limiting cases). It is shown that the model under consideration admits only expanding solutions. As the astrophysical application of the model, a description of rotating curves of spiral galaxies is suggested by using the analytical solutions obtained in the paper.

This work deals with $F(T)$ gravity models driven by real scalar fields with usual and phantom dynamics. We illustrate the results with examples of current interest, and we find some analytical solutions for scale factors and scalar fields. The results indicate that torsion-scalar models also admit the accelerated expansion of the universe.

The model of dark matter is presented where the dark matter is a \emph{classical} gauge field. A spherical symmetric solution of Yang-Mills equation is obtained. The asymptotic behavior of the gauge fields and matter density is investigated. It is shown that the distribution of the matter density allows us interpret it as the dark matter. The fitting of a typical rotational curve with the rotational curve created by the spherical solution of SU(3) Yang-Mills equation is made.