We use observational data from Supernovae (SNIa) Pantheon sample, as well as from direct measurements of the Hubble parameter from the cosmic chronometers (CC) sample, in order to extract constraints on the scenario of Barrow holographic dark energy. The latter is a holographic dark energy model based on the recently proposed Barrow entropy, which arises from the modification of the black-hole surface due to quantum-gravitational effects. We first consider the case where the new deformation exponent $\Delta$ is the sole model parameter, and we show that although the standard value $\Delta=0$, which corresponds to zero deformation, lies within the 1$\sigma$ region, a deviation is favored. In the case where we let both $\Delta$ and the second model parameter to be free we find that a deviation from standard holographic dark energy is preferred. Additionally, applying the Akaike, Bayesian and Deviance Information Criteria, we conclude that the one-parameter model is statistically compatible with $\Lambda$CDM paradigm, and preferred comparing to the two-parameter one. Finally, concerning the present value of the Hubble parameter we find that it is close to the Planck value.

We present modified cosmological scenarios that arise from the application of the "gravity-thermodynamics" conjecture, using the Barrow entropy instead of the usual Bekenstein-Hawking one. The former is a modification of the black hole entropy due to quantum-gravitational effects that deform the black-hole horizon by giving it an intricate, fractal structure. We extract modified cosmological equations which contain new extra terms that constitute an effective dark-energy sector, and which coincide with the usual Friedmann equations in the case where the new Barrow exponent acquires its Bekenstein-Hawking value. We present analytical expressions for the evolution of the effective dark energy density parameter, and we show that the universe undergoes through the usual matter and dark-energy epochs. Additionally, the dark-energy equation-of-state parameter is affected by the value of the Barrow deformation exponent and it can lie in the quintessence or phantom regime, or experience the phantom-divide crossing. Finally, at asymptotically large times the universe always results in the de-Sitter solution.

We formulate Barrow holographic dark energy, by applying the usual holographic principle at a cosmological framework, but using the Barrow entropy instead of the standard Bekenstein-Hawking one. The former is an extended black-hole entropy that arises due to quantum-gravitational effects which deform the black-hole surface by giving it an intricate, fractal form. We extract a simple differential equation for the evolution of the dark energy density parameter, which possesses standard holographic dark energy as a limiting sub-case, and we show that the scenario can describe the universe thermal history, with the sequence of matter and dark energy eras. Additionally, the new Barrow exponent $\Delta$ significantly affects the dark energy equation of state, and according to its value it can lead it to lie in the quintessence regime, in the phantom regime, or experience the phantom-divide crossing during the evolution.

We investigate the validity of the generalized second law of thermodynamics, applying Barrow entropy for the horizon entropy. The former arises from the fact that the black-hole surface may be deformed due to quantum-gravitational effects, quantified by a new exponent $\Delta$. We calculate the entropy time-variation in a universe filled with the matter and dark energy fluids, as well as the corresponding quantity for the apparent horizon. We show that although in the case $\Delta=0$, which corresponds to usual entropy, the sum of the entropy enclosed by the apparent horizon plus the entropy of the horizon itself is always a non-decreasing function of time and thus the generalized second law of thermodynamics is valid, in the case where quantum-gravitational corrections switch on this is not true anymore, and the generalized second law of thermodynamics may be violated, depending on the universe evolution.

We revisit inflation with non-canonical scalar fields by applying deformed-steepness exponential potentials. We show that the resulting scenario can lead to inflationary observables, and in particular to scalar spectral index and tensor-to-scalar ratio, in remarkable agreement with observations. Additionally, a significant advantage of the scenario is that the required parameter values, such as the non-canonicality exponent and scale, as well as the potential exponent and scale, do not need to acquire unnatural values and hence can accept a theoretical justification. Hence, we obtain a significant improvement with respect to alternative schemes, and we present distinct correlations between the model parameters that better fit the data, which can be tested in future probes. This combination of observational efficiency and theoretical justification makes the scenario at hand a good candidate for the description of inflation.

We propose a new approach to investigate inflation in a model-independent way, and in particular to elaborate the involved observables, through the introduction of the "scale factor potential". Through its use one can immediately determine the inflation end, which corresponds to its first (and global) minimum. Additionally, we express the inflationary observables in terms of its logarithm, using as independent variable the e-folding number. As an example, we construct a new class of scalar potentials that can lead to the desired spectral index and tensor-to-scalar ratio, namely $n_s \approx 0.965$ and $r \sim 10^{-4}$ for 60 $e$-folds, in agreement with observations.

A violation of the distance-duality relation is directly linked with a temporal variation of the electromagnetic fine-structure constant. We consider a number of well-studied $f(T)$ gravity models and we revise the theoretical prediction of their corresponding induced violation of the distance-duality relationship. We further extract constraints on the involved model parameters through fine-structure constant variation data, alongside with supernovae data, and Hubble parameter measurements. Moreover, we constrain the evolution of the effective $f(T)$ gravitational constant. Finally, we compare with revised constraints on the phenomenological parametrisations of the violation of the equivalence principle in the electromagnetic sector.

We present a new class of non-singular bounce cosmology free from instabilities, using a single scalar field coupled to gravity within the framework of the Degenerate Higher-Order Scalar-Tensor (DHOST) theories. In this type of scenarios, the gradient instability that widely exists in non-singular bounce cosmologies in the framework of scalar-tensor and Horndeski/Galileon theories is removed by the effects of new operators introduced by the DHOST, due to the modification that they later bring about to the dispersion relation of perturbations. Hence, our results demonstrate that there is indeed a loophole for this type of bounce scenarios to be free from pathologies when primordial perturbations evolve through the bounce phase, and thus the theoretical {\it no-go} theorem for non-singular bounce cosmology of Horndeski/Galileon theories can be delicately evaded in DHOST extensions.

We use data from galaxy-galaxy weak lensing to perform a novel test on General Relativity (GR). In particular, we impose strong constraints using the torsional (teleparallel) formulation of gravity in which the deviation from GR is quantified by a single parameter $\alpha$, an approximation which is always valid at low-redshift Universe and weak gravitational fields. We calculate the difference in the deflection angle and eventually derive the modified Excess Surface Density profile, which is mainly affected at small scales. Hence, confronting the predictions with weak lensing data from the Sloan Digital Sky Survey Data Release 7 we obtain the upper bound on the deviation parameter, which, expressed via the dimensionless percentage in the universe energy content, reads as ${\rm{log_{10}}}\Omega_\alpha\le -18.52_{-0.42}^{+0.80}$({\it stat.})$_{-0.39}^{+1.50}$({\it sys.}) with systematics mainly arising from the modelling. To our knowledge, this is the first time that GR is verified at such an accuracy at the corresponding scales.

We apply Gaussian processes and Hubble function data in $f(T)$ cosmology, to reconstruct for the first time the $f(T)$ form in a model-independent way. In particular, using $H(z)$ datasets coming from cosmic chronometers as well as from the radial BAO method, alongside the latest released local value $H_{0} = 73.52 \pm 1.62$ km/s/Mpc, we reconstruct $H(z)$ and its derivatives, resulting eventually in a reconstructed region for $f(T)$, without any assumption. Although the cosmological constant lies in the central part of the reconstructed region, the obtained mean curve follows a quadratic function. Inspired by this we propose a new $f(T)$ parametrization, i.e. $f(T) = -2\Lambda +\xi T^2$, with $\xi$ the sole free parameter that quantifies the deviation from $\Lambda$CDM cosmology. Additionally, we confront three viable one-parameter $f(T)$ models of the literature, which respectively are the power-law, the square-root exponential, and the exponential one, with the reconstructed $f(T)$ region, and then we extract significantly improved constraints for their model parameters, comparing to the constraints that arise from usual observational analysis. Finally, we argue that since we are using the direct Hubble measurements and the local value for $H_0$ in our analysis, with the above reconstruction of $f(T)$, the $H_0$ tension can be efficiently alleviated.

We investigate the cosmological applications of scalar-tensor theories that arise effectively from the Lorentz fiber bundle of a Finsler-like geometry. We first show that the involved nonlinear connection induces a new scalar degree of freedom and eventually a scalar-tensor theory. Using both a holonomic and a nonholonomic basis, we show the appearance of an effective dark-energy sector, which additionally acquires an explicit interaction with the matter sector, arising purely from the internal structure of the theory. Applying the theory at late times we find that we can obtain the thermal history of the Universe, namely the sequence of matter and dark-energy epochs, and moreover the effective dark-energy equation-of-state parameter can be quintessencelike, phantomlike, or experience the phantom-divide crossing during cosmological evolution. Furthermore, applying the scenario at early times, we see that one can acquire an exponential de Sitter solution as well as obtain an inflationary realization with the desired scale-factor evolution. These features arise purely from the intrinsic geometrical structure of Finsler-like geometry and reveal the capabilities of the construction.

We investigate the cosmological applications of Myrzakulov $F(R,T)$ gravity. In this theory ones uses a specific but non-special connection, and thus both curvature and torsion are dynamical fields related to gravity. We introduce a parametrization that quantifies the deviation of curvature and torsion scalars form their corresponding values obtained using the special Levi-Civita and Weitzenb\"{o}ck connections, and we extract the cosmological field equations following the mini-super-space procedure. Even for the simple case where the action of the theory is linear in $R$ and $T$, we find that the Friedmann equations contain new terms of geometrical origin, reflecting the non-special connection. Applying the theory at late times we find that we can acquire the thermal history of the universe, where dark energy can be quintessence-like or phantom-like, or behave exactly as a cosmological constant and thus reproducing $\Lambda$CDM cosmology. Furthermore, we show that these features are obtained for other Lagrangian choices, too. Finally, early-time application leads to the de Sitter solution, as well as to an inflationary realization with the desired scale-factor evolution.

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.

We investigate interacting scenarios which belong to a wider class, since they include a dynamical dark energy component whose equation of state follows various one-parameter parametrizations. We confront them with the latest observational data from Cosmic Microwave Background (CMB), Joint light-curve (JLA) sample from Supernovae Type Ia, Baryon Acoustic Oscillations (BAO), Hubble parameter measurements from Cosmic Chronometers (CC) and a gaussian prior on the Hubble parameter $H_0$. In all examined scenarios we find a non-zero interaction, nevertheless the non-interacting case is allowed within 2$\sigma$. Concerning the current value of the dark energy equation of state for all combination of datasets it always lies in the phantom regime at more than two/three standard deviations. Finally, for all interacting models, independently of the combination of datasets considered, the estimated values of the present Hubble parameter $H_0$ are greater compared to the $\Lambda$CDM-based Planck's estimation and close to the local measurements, thus alleviating the $H_0$ tension.

We investigate various scenarios which include interaction forms between dark matter and dark energy that exhibit sign reverse, namely where the transfer of energy between the dark fluids changes sign during evolution. We study the large-scale inhomogeneities in such interacting scenarios and we confront them with the latest astronomical data. Our analysis shows that the sign-changeable interaction models are able to produce stable perturbations. Additionally, the data seem to slightly favor a non-zero interaction, however, within $1\sigma$ confidence level (CL) the scenarios cannot be distinguished from non-interacting cosmologies. We find that the best-fit value of the dark-energy equation-of-state parameter lies in the phantom regime, while the quintessence region is also allowed nevertheless at more than 2$\sigma$ CL. Examining the effect of the interaction on the CMB TT and matter power spectra we show that while from the simple spectra it is hard to distinguish the interacting case from $\Lambda$CDM scenario, in the residual graphs the interaction is indeed traceable. Moreover, we find that sign-changeable interaction models can reconcile the $H_0$ tension, however the $\sigma_8$ tension is still persisting. Finally, we examine the validity of the laws of thermodynamics and we show that the generalized second law is always satisfied, while the second derivative of the total entropy becomes negative at late times which implies that the universe tends towards thermodynamic equilibrium.

We use observational data from Supernovae (SNIa) Pantheon sample, from direct Hubble constant measurements with cosmic chronometers (CC), from the Cosmic Microwave Background shift parameter $\text{CMB}_{\text{shift}}$, and from redshift space distortion ($f\sigma_8$) measurements, in order to constrain $f(T)$ gravity. We do not follow the common $\gamma$ parameterization within the semi-analytical approximation of the growth rate, in order to avoid model-dependent uncertainties. Up to our knowledge this is the first time that $f(T)$ gravity is analyzed within a Bayesian framework, and with background and perturbation behaviour considered jointly. We show that all three examined $f(T)$ models are able to describe adequately the $f\sigma_8$ data. Furthermore, applying the Akaike, Bayesian and Deviance Information Criteria, we conclude that all considered models are statistically equivalent, however the most efficient candidate is the exponential model, which additionally presents a small deviation from $\Lambda$CDM paradigm.

We investigate the bounce realization arising from the application of the holographic principle in the early universe, inspired by its well-studied late-time application. We first consider as Infrared cutoffs the particle and future event horizons, and we show that the decrease of the horizons at early times naturally increases holographic energy density at bouncing scales, while we additionally obtain the necessary null energy condition violation. Furthermore, adding a simple correction to the horizons due to the Ultraviolet cutoff we analytically obtain improved nonsingular bouncing solutions, in which the value of the minimum scale factor is controlled by the UV correction. Finally, we construct generalized scenarios, arisen from the use of extended Infrared cutoffs, and as specific examples we consider cutoffs that can reproduce $F(R)$ gravity, and the bounce realization within it.

Cosmological tensions can arise within $\Lambda$CDM scenario amongst different observational windows, such as measurements of the $H_0$ and $\sigma_8$ parameters. These tensions, if finally confirmed by measurements, may indicate new physics beyond the standard paradigm. In this Letter, we report how to alleviate both the $H_0$ and $\sigma_8$ tensions simultaneously within torsional gravity from the perspective of effective field theory (EFT). We apply the EFT approach, which allows to investigate the evolution equations at the background and perturbation levels in a systematic way, and examine the conditions followed by the coefficients of various possibly involved operators such that cosmological tensions can be relaxed. Following these observations we construct concrete models of Lagrangians of torsional gravity. Specifically, we consider the parametrization $f(T)=-T -2\Lambda/M_P^2 +\alpha T^\beta$, where two out of the three parameters are independent (in which an additional term of the form $cT^{1/2}$ can be added). This model can efficiently fit observations solving all statistical tensions. To our knowledge, this is the first time where a modified gravity theory can alleviate both $H_0$ and $\sigma_8$ tensions simultaneously, hence offering an additional argument in favor of gravitational modification.

We perform for the first time a dynamical system analysis of both the background and perturbation equations, of $\Lambda$CDM cosmology and quintessence scenario with an exponential potential. In the former case the perturbations do not change the stability of the late-time attractor of the background equations, and the system still results in the dark-energy dominated, de Sitter solution, having passed from the correct dark-matter era with $\gamma\approx6/11$. However, in the case of quintessence the incorporation of perturbations changes the stability and properties of the background evolution, and the only conditionally stable points present either an exponentially increasing matter clustering not favored by observations, or Laplacian instabilities, and thus not physically interesting. This result is a severe disadvantage of quintessence cosmology comparing to $\Lambda$CDM paradigm.

Conformal symmetries appear in many parts of physics and play a unique role in exploring the Universe. In this work, we consider the possibility of constructing conformal theories of gravity in the Symmetric Teleparallel Gravity framework, where gravitation is expressed through non-metricity rather than curvature or torsion. We demonstrate that it is possible to formulate a family of conformal theories with second-order field equations and have the metric tensor as the fundamental variable. Moreover, we consider the addition of a scalar field to the theory, and compare the results to the Teleparallel Gravity setting. Finally, we present the scalar-tensor representation of modified symmetric teleparallel theories, and amongst others, we show that there is an interesting sub-case that has the possibility of playing the role of dark energy.

Symmetric teleparallel gravity (STG) offers an interesting avenue to formulate a theory of gravitation that relies neither on curvature nor torsion but only on non-metricity Q. Given the growing number of observations of gravitational waves (GWs) and their use to explore gravitational theories, in this work we investigate the GWs in various extensions of STG, focusing on their speed and polarization. For the simple STG, for theories that arise from the generalized irreducible decomposition of STG, and for f(Q) gravity, we obtain the same speed and polarizations with general relativity. For scalar - non-metricity theories, where a scalar field is nonminimally coupled to f(Q), we find that GWs propagate in general with a speed different than the one of light. Finally, for the case of f(Q,B) gravity we find that new polarizations do appear.

A model of inflation realization driven by fermions with curvature-dependent mass is studied. Such a term is derived from the Covariant Canonical Gauge Theory of gravity (CCGG) incorporating Dirac fermions. We obtain an initial de Sitter phase followed by a successful exit, and moreover, we acquire the subsequent thermal history, with an effective matter era, followed finally by a dark-energy epoch. This behavior is a result of the effective `weakening' of gravity at early times, due to the decreased curvature-dependent fermion mass. Through the perturbation level, we obtain the scalar spectral index and the tensor-to-scalar ratio, which are in agreement with the Planck observations. The efficiency of inflation from fermions with curvature-dependent mass, at both the background and perturbation level, reveals the capabilities of the scenario and makes it a good candidate for the description of nature.

We apply the holographic principle at the early universe, obtaining an inflation realization of holographic origin. Such a consideration has equal footing with its well-studied late-time application, and moreover the decrease of the horizons at early times naturally increases holographic energy density at inflationary scales. Taking as Infrared cutoff the particle or future event horizons, and adding a simple correction due to the Ultraviolet cutoff, whose role is non-negligible at the high energy scales of inflation, we result in a holographic inflation scenario that is very efficient in incorporating inflationary requirements and predictions. We first extract analytically the solution of the Hubble function in an implicit form, which gives a scale factor evolution of the desired e-foldings. Furthermore, we analytically calculate the Hubble slow-roll parameters and then the inflation-related observables, such as the scalar spectral index and its running, the tensor-to-scalar ratio, and the tensor spectral index. Confronting the predictions with Planck 2018 observations we show that the agreement is perfect and in particular deep inside the 1$\sigma$ region.

We construct cubic gravity and its $f(P)$ extension and we investigate their early- and late-time cosmological applications. Cubic gravity is based on a particular invariant $P$, constructed from cubic contractions of the Riemann tensor, under three requirements: (i) the resulting theory possesses a spectrum identical to that of general relativity, (ii) it is neither topological nor trivial in four dimensions, and (iii) it is defined such that it is independent of the dimensions. Relaxing the last condition and restricting the parameters of cubic gravity we can obtain second-order field equations in a cosmological background. We show that at early times one can obtain inflationary, de Sitter solutions, which are driven by an effective cosmological constant constructed purely from the cubic terms of the simple cubic or $f(P)$ gravity. Concerning late-time evolution, the new terms constitute an effective dark-energy sector and we show that the Universe experiences the usual thermal history and the onset of late-time acceleration. In the case of $f(P)$ gravity, depending on the choice of parameters, we find that the dark-energy equation-of-state parameter can be quintessencelike, phantomlike or it can experience the phantom-divide crossing during the evolution, even if an explicit cosmological constant is absent.

We show that a nonsingular bounce, free of ghosts and gradient instabilities, can be realized in the framework of Horndeski or generalized Galileon cosmology. In particular, we first review that the theoretical no-go theorem, which states that the above is impossible, is based on two very strong assumptions, namely that a particular quantity cannot be discontinuous during the bounce, and that there is only one bounce. However, as we show in the present work, the first assumption not only can be violated in a general Horndeski/Galileon scenario, but also it is necessarily violated at the bounce point within the subclass of Horndeski/Galileon gravity in which $K(\phi,X)$ becomes zero at $X=0$. Additionally, concerning the second assumption, which is crucial in improved versions of the theorem which claim that even if a nonlinear free of pathologies can be realized it will lead to pathologies in the infinite past or infinite future, we show that if needed it can be evaded by considering cyclic cosmology, with an infinite sequence of nonsingular bounces free of pathologies, which forbids the universe to reach the "problematic" regime at infinite past or infinite future. Finally, in order to make the analysis more transparent we provide explicit examples where nonsingular bounces without theoretical pathologies can be achieved.

We present a modified cosmological scenario that arises from the application of non-extensive thermodynamics with varying exponent. We extract the modified Friedmann equations, which contain new terms quantified by the non-extensive exponent, possessing standard $\Lambda$CDM cosmology as a subcase. Concerning the universe evolution at late times we obtain an effective dark energy sector, and we show that we can acquire the usual thermal history, with the successive sequence of matter and dark-energy epochs, with the effective dark-energy equation-of-state parameter being in the quintessence or in the phantom regime. The interesting feature of the scenario is that the above behaviors can be obtained even if the explicit cosmological constant is set to zero, namely they arise purely from the extra terms. Additionally, we confront the model with Supernovae type Ia and Hubble parameter observational data, and we show that the agreement is very good. Concerning the early-time universe we obtain inflationary de Sitter solutions, which are driven by an effective cosmological constant that includes the new terms of non-extensive thermodynamics. This effective screening can provide a description of both inflation and late-time acceleration with the same parameter choices, which is a significant advantage.

We investigate the bounce realization in the framework of generalized modified gravities arising from Finsler and Finsler-like geometries. In particular, the richer intrinsic geometrical structure is reflected in the appearance of extra degrees of freedom in the Friedmann equations that can drive the bounce. We examine various Finsler and Finsler-like constructions. In the cases of general very special relativity as well as of Finsler-like gravity on the tangent bundle we show that a bounce cannot be easily obtained. However, in the Finsler-Randers space the induced scalar anisotropy can fulfill the bounce conditions and bouncing solutions are easily obtained. Finally, for the general class of theories that include a nonlinear connection a new scalar field is induced, leading to a scalar-tensor structure that can easily drive a bounce. These features reveal the capabilities of Finsler and Finsler-like geometries.

The phenomenological parametrizations of dark-energy (DE) equation of state can be very helpful, since they allow for the investigation of its cosmological behavior despite the fact that its underlying theory is unknown. However, although there has been a large amount of research on DE parametrizations which involve two or more free parameters, the one-parameter parametrizations seem to be underestimated. We perform a detailed observational confrontation of five one-parameter DE models, with observational data from cosmic microwave background (CMB), Joint light-curve analysis sample from Supernovae Type Ia observations (JLA), baryon acoustic oscillations (BAO) distance measurements, and cosmic chronometers (CC). We find that all models favor a phantom DE equation of state at present time, while they lead to $H_0$ values in perfect agreement with its direct measurements and therefore they offer an alleviation to the $H_0$-tension. Finally, performing a Bayesian analysis we show that although $\Lambda$CDM cosmology is still favored, one-parameter DE models have similar or better efficiency in fitting the data comparing to two-parameter DE parametrizations, and thus they deserve a thorough investigation.

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 construct modified cosmological scenarios through the application of the first law of thermodynamics in the universe horizon, but using the generalized, nonextensive Tsallis entropy instead of the usual Bekenstein-Hawking one. We result to modified cosmological equations that possess the usual ones as a particular limit, but which in the general case contain extra terms that appear for the first time, that constitute an effective dark energy sector quantified by the nonextensive parameter $\delta$. When the matter sector is dust, we extract analytical expressions for the dark energy density and equation-of-state parameters, and we extend these solutions to the case where radiation is present too. We show that the universe exhibits the usual thermal history, with the sequence of matter and dark-energy eras, and according to the value of $\delta$ the dark-energy equation-of-state parameter can be quintessence-like, phantom-like, or experience the phantom-divide crossing during the evolution. Even in the case where the explicit cosmological constant is absent, the scenario at hand can very efficiently mimic $\Lambda$CDM cosmology, and is in excellent agreement with Supernovae type Ia observational data.

We investigate the new observational constraints on $f(T)$ gravity that arise from the effects of primordial gravitational waves (GWs) on the cosmic microwave background (CMB) anisotropies and the BB spectrum. We first show that on the GWs propagation in $f(T)$ gravity we obtain only an amplitude modification and not a phase one, comparing to the case of general relativity in the background of $\Lambda$CDM cosmology. Concerning primordial GWs we find that the more the model departs from general relativity the larger is the GW amplitude decay, and thus a possible future detection would bring the viable $f(T)$ gravity models five orders of magnitude closer to $\Lambda$CDM cosmology comparing to standard cosmological constraints. Additionally, we use the CLASS code and both data from the Planck probe, as well as forecasts from the near-future CORE collaboration, and we show that possible non-trivial constraints on the tensor-to-scalar ratio would offer a clear signature of $f(T)$ gravity. Finally, we discuss on the possibility to use the properties of the GWs that arise from neutron stars mergers in order to extract additional constrains on the theory.

We investigate the generalized Chevallier-Polarski-Linder (CPL) parametrization, which contains the pivoting redshift $z_p$ as an extra free parameter. We use various data combinations from cosmic microwave background (CMB), baryon acoustic oscillations (BAO), redshift space distortion (RSD), weak lensing (WL), joint light curve analysis (JLA), cosmic chronometers (CC), and we include a Gaussian prior on the Hubble constant value, in order to extract the observational constraints on various quantities. For the case of free $z_p$ we find that for all data combinations it always remains unconstrained, and there is a degeneracy with the current value of the dark energy equation of state $w_0$. For the case where $z_p$ is fixed to specific values, and for the full data combination, we find that with increasing $z_p$ the mean value of $w_0$ slowly moves into the phantom regime, however the cosmological constant is always allowed within 1$\sigma$ confidence-level. However, the significant effect is that with increasing $z_p$ the correlations between $w_0$ and $w_a$ change from negative to positive, with the case $z_p =0.35$ corresponding to no correlation. This feature indeed justifies why a non-zero pivoting redshift should be taken into account.

We develop the effective field theory approach to torsional modified gravities, a formalism that allows for the systematic investigation of the background and perturbation levels separately. Starting from the usual effective field theory approach to curvature-based gravity, we suitably generalize it at the background level by including terms of the contracted torsion tensor, and at the perturbation level by including pure torsion perturbative terms and mixed perturbative terms of torsion and curvature. Having constructed the effective field theory action of general torsional modified gravity, amongst others we focus on $f(T)$ gravity and we perform a cosmological application. We investigate the scalar perturbations up to second order, and we derive the expressions of the Newtonian constant and the post Newtonian parameter $\gamma$. Finally, we apply this procedure to two specific and viable $f(T)$ models, namely the power-law and the exponential ones, introducing a new parameter that quantifies the deviation from general relativity and depends on the model parameters. Since this parameter can be expressed in terms of the scalar perturbation mode, a precise measurement of its evolution could be used as an alternative way to impose constraints on $f(T)$ gravity and break possible degeneracies between different $f(T)$ models.

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 perform a detailed confrontation of various oscillating dark-energy parame-trizations with the latest sets of observational data. In particular, we use data from Joint Light Curve analysis (JLA) sample from Supernoave Type Ia, Baryon Acoustic Oscillations (BAO) distance measurements, Cosmic Microwave Background (CMB) observations, redshift space distortion, weak gravitational lensing, Hubble parameter measurements from cosmic chronometers, and we impose constraints on four oscillating models. From the analyses we find that the best-fit characters of almost all models are bent towards the phantom region, nevertheless in all of them the quintessential regime is also allowed within 1$\sigma$ confidence-level. Furthermore, the deviations from $\Lambda$CDM cosmology are not significant, however for two of the models they could be visible at large scales, through the impact on the temperature anisotropy of the CMB spectra and on the matter power spectra. Finally, we peform the Bayesian analysis, which shows that the current observational data support the $\Lambda$CDM paradigm over this set of oscillating dark-energy parametrizations.

We extract observational constraints on $f(T)$ gravity, using the recently proposed statistical method which is not affected by the value of $H_0$ and thus it bypasses the problem of the disagreement in its exact numerical value between Planck and direct measurements. We use direct measurements of the Hubble parameter with the corresponding covariance matrix, and for completeness we perform a joint analysis using the latest data from Supernovae type Ia based on JLA sample, quasi-stellar objects, and Cosmic Microwave Background shift parameter from Planck. We analyze a large family of $f(T)$ models, and we compare the fitting results with $\Lambda$CDM cosmology using the AIC statistical test. Utilizing only the Hubble parameter data we find that in the case of the power-law $f(T)$ model a small but non-zero deviation from $\Lambda$CDM cosmology is slightly favored at 1-$\sigma$, nevertheless the corresponding AIC value shows a statistical equivalence with it. Finally, the join analysis reveals that all $f(T)$ models are very efficient and in very good agreement with observations.

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 investigate the gravitational waves and their properties in various modified teleparallel theories, such as $f(T), f(T,B)$ and $f(T,T_G)$ gravities. We perform the perturbation analysis both around a Minkowski background, as well as in the case where a cosmological constant is present, and for clarity we use both the metric and tetrad language. For $f(T)$ gravity we verify the result that no further polarization modes comparing to general relativity are present at first-order perturbation level, and we show that in order to see extra modes one should look at third-order perturbations. For non-trivial $f(T,B)$ gravity, by examining the geodesic deviation equations, we show that extra polarization models, namely the longitudinal and breathing modes, do appear at first-order perturbation level, and the reason for this behavior is the fact that although the first-order perturbation does not have any effect on $T$, it does affect the boundary term $B$. Finally, for $f(T,T_G)$ gravity we show that at first-order perturbations the gravitational waves exhibit the same behavior to those of $f(T)$ gravity. Since different modified teleparallel theories exhibit different gravitational wave properties, the dvancing gravitational-wave astronomy would help to alleviate the degeneracy not only between curvature and torsional modified gravity but also between different subclasses of modified teleparallel gravities.

We present a model of holographic dark energy in which the Infrared cutoff is determined by both the Ricci and the Gauss-Bonnet invariants. Such a construction has the significant advantage that the Infrared cutoff, and consequently the holographic dark energy density, does not depend on the future or the past evolution of the universe, but only on its current features, and moreover it is determined by invariants, whose role is fundamental in gravitational theories. We extract analytical solutions for the behavior of the dark energy density and equation-of-state parameters as functions of the redshift. These reveal the usual thermal history of the universe, with the sequence of radiation, matter and dark energy epochs, resulting in the future to a complete dark energy domination. The corresponding dark energy equation-of-state parameter can lie in the quintessence or phantom regime, or experience the phantom-divide crossing during the cosmological evolution, and its asymptotic value can be quintessence-like, phantom-like, or be exactly equal to the cosmological-constant value. Finally, we extract the constraints on the model parameters that arise from Big Bang Nucleosynthesis.

We use dynamical system methods to explore the general behaviour of $f(T)$ cosmology. In contrast to the standard applications of dynamical analysis, we present a way to transform the equations into a one-dimensional autonomous system, taking advantage of the crucial property that the torsion scalar in flat FRW geometry is just a function of the Hubble function, thus the field equations include only up to first derivatives of it, and therefore in a general $f(T)$ cosmological scenario every quantity is expressed only in terms of the Hubble function. The great advantage is that for one-dimensional systems it is easy to construct the phase space portraits, and thus extract information and explore in detail the features and possible behaviours of $f(T)$ cosmology. We utilize the phase space portraits and we show that $f(T)$ cosmology can describe the universe evolution in agreement with observations, namely starting from a Big Bang singularity, evolving into the subsequent thermal history and the matter domination, entering into a late-time accelerated expansion, and resulting to the de Sitter phase in the far future. Nevertheless, $f(T)$ cosmology can present a rich class of more exotic behaviours, such as the cosmological bounce and turnaround, the phantom-divide crossing, the Big Brake and the Big Crunch, and it may exhibit various singularities, including the non-harmful ones of type II and type IV. We study the phase space of three specific viable $f(T)$ models offering a complete picture. Moreover, we present a new model of $f(T)$ gravity that can lead to a universe in agreement with observations, free of perturbative instabilities, and applying the Om(z) diagnostic test we confirm that it is in agreement with the combination of SNIa, BAO and CMB data at 1$\sigma$ confidence level.

We present a unified description of the dark matter and the dark energy sectors, in the framework of shift-symmetric generalized Galileon theories. Considering a particular combination of terms in the Horndeski Lagrangian in which we have not introduced a cosmological constant or a matter sector, we obtain an effective unified cosmic fluid whose equation of state $w_U$ is zero during the whole matter era, namely from redshifts $z\sim3000$ up to $z\sim2-3$. Then at smaller redshifts it starts decreasing, passing the bound $w_U=-1/3$, which marks the onset of acceleration, at around $z\sim0.5$. At present times it acquires the value $w_U=-0.7$. Finally, it tends toward a de-Sitter phase in the far future. This behaviour is in excellent agreement with observations. Additionally, confrontation with Supernovae type Ia data leads to a very efficient fit. Examining the model at the perturbative level, we show that it is free from pathologies such as ghosts and Laplacian instabilities, at both scalar and tensor sectors, at all times.

The combined observation of GW170817 and its electromagnetic counterpart GRB170817A reveals that gravitational waves propagate at the speed of light in high precision. We apply the effective field theory approach to investigate the experimental consequences for the theory of $f(T)$ gravity. We find that the speed of gravitational waves within $f(T)$ gravity is exactly equal to the light speed, and hence the constraints from GW170817 and GRB170817A are trivially satisfied. The results are verified through the standard analysis of cosmological perturbations. Nevertheless, by examining the dispersion relation and the frequency of cosmological gravitational waves, we observe a deviation from the results of General Relativity, quantified by a new parameter. Although its value is relatively small in viable $f(T)$ models, its possible future measurement in advancing gravitational-wave astronomy would be the smoking gun of testing this type of modified gravity.

We investigate the realization of the emergent universe scenario in theories with natural UV cutoffs, namely a minimum length and a maximum momentum, quantified by a new deformation parameter in the generalized uncertainty principle. We extract the Einstein static universe solutions and we examine their stability through a phase-space analysis. As we show, the role of the new deformation parameter is crucial in a twofold way: Firstly, it leads to the appearance of new Einstein static universe critical points, that are absent in standard cosmology. Secondly, it provides a way for a graceful exit from the Einstein static universe into the expanding thermal history, that is needed for a complete and successful realization of the emergent universe scenario.

We study quintessential inflation using a generalized exponential potential $V(\phi)\propto exp(-\lambda \phi^n/Mpl^n), n>1$, the model admits slow-roll inflation at early times and leads to close-to-scaling behaviour in the post inflationary era with an exit to dark energy at late times. We present detailed investigations of the inflationary stage in the light of the Planck 2015 results, study post-inflationary dynamics and analytically confirm the existence of an approximately scaling solution. Additionally, assuming that standard massive neutrinos are non-minimally coupled, makes the field $\phi$ dominant once again at late times giving rise to present accelerated expansion of the Universe. We derive observational constraints on the field and time-dependent neutrino masses. In particular, for $n=6 (8)$, the parameter $\lambda$ is constrained to be,$\log \lambda > -7.29 (-11.7)$; the model produces the spectral index of the power spectrum of primordial scalar (matter density) perturbations as $ n_s = 0.959 \pm 0.001 (0.961 \pm 0.001)$ and tiny tensor-to-scalar ratio, $r<1.72 \times 10^{-2} (2.32 \times 10^{-2})$ respectively. Consequently, the upper bound on possible values of the sum of neutrino masses $\Sigma m_{\nu} \lesssim 2.5$ eV significantly enhances compared to that in the standard $\Lambda$CDM model.

We investigate the inflationary realization in the context of unimodular $F(T)$ gravity, which is based on the $F(T)$ modification of teleparallel gravity, in which one imposes the unimodular condition through the use of Lagrange multipliers. We develop the general reconstruction procedure of the $F(T)$ form that can give rise to a given scale-factor evolution, and then we apply it in the inflationary regime. We extract the Hubble slow-roll parameters that allow us to calculate various inflation-related observables, such as the scalar spectral index and its running, the tensor-to-scalar ratio, and the tensor spectral index. Then, we examine the particular cases of de Sitter and power-law inflation, of Starobinsky inflation, as well as inflation in a specific model of unimodular $F(T)$ gravity. As we show, in all cases the predictions of our scenarios are in a very good agreement with Planck observational data. Finally, inflation in unimodular $F(T)$ gravity has the additional advantage that it always allows for a graceful exit for specific regions of the model parameters.

We use observations related to the variation of fundamental constants, in order to impose constraints on the viable and most used $f(T)$ gravity models. In particular, for the fine-structure constant we use direct measurements obtained by different spectrographic methods, while for the effective Newton's constant we use a model-dependent reconstruction, using direct observational Hubble parameter data, in order to investigate its temporal evolution. We consider two $f(T)$ models and we quantify their deviation from $\Lambda$CDM cosmology through a sole parameter. Our analysis reveals that this parameter can be slightly different from its $\Lambda$CDM value, however the best-fit value is very close to the $\Lambda$CDM one. Hence, $f(T)$ gravity is consistent with observations, nevertheless, as every modified gravity, it may exhibit only small deviations from $\Lambda$CDM cosmology, a feature that must be taken into account in any $f(T)$ model-building.

We investigate the bounce and cyclicity realization in the framework of weakly broken galileon theories. We study bouncing and cyclic solutions at the background level, reconstructing the potential and the galileon functions that can give rise to a given scale factor, and presenting analytical expressions for the bounce requirements. We proceed to a detailed investigation of the perturbations, which after crossing the bouncing point give rise to various observables, such as the scalar and tensor spectral indices and the tensor-to-scalar ratio. Although the scenario at hand shares the disadvantage of all bouncing models, namely that it provides a large tensor-to-scalar ratio, introducing an additional light scalar significantly reduces it through the kinetic amplification of the isocurvature fluctuations.

We use BBN observational data on primordial abundance of ${}^4He$ to constrain f(T) gravity. The three most studied viable $f(T)$ models, namely the power law, the exponential and the square-root exponential are considered, and the BBN bounds are adopted in order to extract constraints on their free parameters. For the power-law model, we find that the constraints are in agreement with those acquired using late-time cosmological data. For the exponential and the square-root exponential models, we show that for realiable regions of parameters space they always satisfy the BBN bounds. We conclude that viable f(T) models can successfully satisfy the BBN constraints.

We construct gravitational modifications that go beyond Horndeski, namely theories with extended nonminimal derivative couplings, in which the coefficient functions depend not only on the scalar field but also on its kinetic energy. Such theories prove to be ghost-free in a cosmological background. We investigate the early-time cosmology and show that a de Sitter inflationary phase can be realized as a pure result of the novel gravitational couplings. Additionally, we study the late-time evolution, where we obtain an effective dark energy sector which arises from the scalar field and its extended couplings to gravity. We extract various cosmological observables and analyse their behavior at small redshifts for three choices of potentials, namely, for the exponential, the power-law, and the Higgs potential. We show that the Universe passes from deceleration to acceleration in the recent cosmological past, while the effective dark-energy equation-of-state parameter tends to the cosmological-constant value at present. Finally, the effective dark energy can be phantom-like, although the scalar field is canonical, which is an advantage of the model.

We use the recently released cosmic chronometer data and the latest measured value of the local Hubble parameter, combined with the latest joint light curves of Supernovae Type Ia, and Baryon Acoustic Oscillation distance measurements, in order to impose constraints on the viable and most used $f(R)$ gravity models. We consider four $f(R)$ models, namely the Hu-Sawicki, the Starobinsky, the Tsujikawa, and the exponential one, and we parametrize them introducing a distortion parameter $b$ that quantifies the deviation from $\Lambda$CDM cosmology. Our analysis reveals that a small but non-zero deviation from $\Lambda$CDM cosmology is slightly favored, with the corresponding fittings exhibiting very efficient $AIC$ and $BIC$ Information Criteria values. Clearly, $f(R)$ gravity is consistent with observations, and it can serve as a candidate for modified gravity.

We investigate how baryogenesis can occur by the presence of an $f(T)$-related gravitational term. We study various cases of $f(T)$ gravity and we discuss in detail the effect of the novel terms on the baryon-to-entropy ratio. Additionally, we study the constraints imposed by the observational values of the baryon-to-entropy ratio and we discuss how more generalized cosmologies can contribute successfully, in a viable and consistent way, in the gravitational baryogenesis mechanism.

We discuss the cosmological implications of an extended Brans-Dicke theory presented recently, in which there is an energy exchange between the scalar field and ordinary matter, determined by the theory. A new mass scale is generated in the theory which modifies the Friedmann equations with field-dependent corrected kinetic terms. In a radiation universe the general solutions are found and there are branches with complete removal of the initial singularity, while at the same time a transient accelerating period can occur within deceleration. Entropy production is also possible in the early universe. In the dust era, late-times acceleration has been found numerically in agreement with the correct behaviour of the density parameters and the dark energy equation of state, while the gravitational constant has only a slight variation over a large redshift interval in agreement with observational bounds.

We construct $F(T,\left(\nabla{T}\right)^2,\Box {T})$ gravitational modifications, which are novel classes of modified theories arising from higher-derivative torsional terms in the action, and are different than their curvature analogue. Applying them in a cosmological framework we obtain an effective dark energy sector that comprises of the novel torsional contributions. We perform a detailed dynamical analysis for two specific examples, extracting the stable late-time solutions and calculating the corresponding observables. We show that the thermal history of the universe can be reproduced, and it can result in a dark-energy dominated, accelerating universe, where the dark-energy equation-of-state parameter lies in the quintessence regime, or may exhibit the phantom-divide crossing during the cosmological evolution. Finally, the scale factor behaves asymptotically either as a power-law or as an exponential, in agreement with observations.

Over the past decades, the role of torsion in gravity has been extensively investigated along the main direction of bringing gravity closer to its gauge formulation and incorporating spin in a geometric description. Here we review various torsional constructions, from teleparallel, to Einstein-Cartan, and metric-affine gauge theories, resulting in extending torsional gravity in the paradigm of f(T) gravity, where f(T) is an arbitrary function of the torsion scalar. Based on this theory, we further review the corresponding cosmological and astrophysical applications. In particular, we study cosmological solutions arising from f(T) gravity, both at the background and perturbation levels, in different eras along the cosmic expansion. The f(T) gravity construction can provide a theoretical interpretation of the late-time universe acceleration, and it can easily accommodate with the regular thermal expanding history including the radiation and cold dark matter dominated phases. Furthermore, if one traces back to very early times, a sufficiently long period of inflation can be achieved and hence can be investigated by cosmic microwave background observations, or alternatively, the Big Bang singularity can be avoided due to the appearance of non-singular bounces. Various observational constraints, especially the bounds coming from the large-scale structure data in the case of f(T) cosmology, as well as the behavior of gravitational waves, are described in detail. Moreover, the spherically symmetric and black hole solutions of the theory are reviewed. Additionally, we discuss various extensions of the f(T) paradigm. Finally, we consider the relation with other modified gravitational theories, such as those based on curvature, like f(R) gravity, trying to enlighten the subject of which formulation might be more suitable for quantization ventures and cosmological applications.

We show that the Pauli exclusion principle in a system of $M0$-branes can give rise to the expansion and contraction of the universe which is located on an $M3$-brane. We start with a system of $M0$-branes with high symmetry, which join mutually and form pairs of $M1$-anti-$M1$-branes. The resulting symmetry breaking creates gauge fields that live on the $M1$-branes and play the role of graviton tensor modes, which induce an attractive force between the $M1$ and anti-$M1$ branes. Consequently, the gauge fields that live on the $M1$-branes, and the scalar fields which are attached symmetrically to all parts of these branes, decay to fermions that attach anti-symmetrically to the upper and lower parts of the branes, and hence the Pauli exclusion principle emerges. By closing $M1$-branes mutually, the curvatures produced by parallel spins will be different from the curvatures produced by anti-parallel spins, and this leads to an inequality between the number of degrees of freedom on the boundary surface and the number of degrees of freedom in the bulk region. This behavior is inherited in the $M3$-brane on which the universe is located, and hence this leads to the emergence of the universe expansion and contraction. In this sense, the Pauli exclusion principle rules the cosmic dynamics.

We use the local value of the Hubble constant recently measured with 2.4% precision, as well as the latest compilation of cosmic chronometers data, together with standard probes such as Supernovae Type Ia and Baryon Acoustic Oscillation distance measurements, in order to impose constraints on the viable and most used f(T) gravity models, where T is the torsion scalar in teleparallel gravity. In particular, we consider three f(T) models with two parameters, out of which one is independent, and we quantify their deviation from $\Lambda$CDM cosmology through a sole parameter. Our analysis reveals that for one of the models a small but non-zero deviation from $\Lambda$CDM cosmology is slightly favored, while for the other models the best fit is very close to $\Lambda$CDM scenario. Clearly, f(T) gravity is consistent with observations, and it can serve as a candidate for modified gravity.

We use the latest compilation of observational Hubble parameter measurements estimated with the differential evolution of cosmic chronometers, in combination with the local value of the Hubble constant recently measured with 2.4% precision, to constrain the cosmological scenario where dark energy interacts directly with the dark matter sector. To diminish the degeneracy between the parameters we additionally consider standard probes, such as Supernovae Type Ia from joint light curves (JLA) sample, Baryon Acoustic Oscillation distance measurements (BAO), and cosmic microwave background data from Planck 2015 estimations. Our analysis shows that the direct interaction between dark energy and dark matter is mildly flavored, while the dark energy equation-of-state parameter is w < - 1 at 3$\sigma$ confidence level.

We investigate the cosmological implications of a new class of modified gravity, where the field equations generically include higher-order derivatives of the matter fields, arising from the introduction of non-dynamical auxiliary fields in the action. Imposing a flat, homogeneous and isotropic geometry we extract the Friedmann equations, obtaining an effective dark-energy sector containing higher derivatives of the matter energy density and pressure. For the cases of dust, radiation, and stiff matter we analyze the cosmological behavior, finding accelerating, de Sitter, and non-accelerating phases, dominated by matter or dark energy. Additionally, the effective dark-energy equation-of-state parameter can be quintessence-like, cosmological-constant-like, or even phantom-like. The detailed study of these scenarios may provide signatures that could distinguish them from other candidates of modified gravity.

We investigate the cosmological applications of new gravitational scalar-tensor theories, which are novel modifications of gravity possessing 2+2 propagating degrees of freedom, arising from a Lagrangian that includes the Ricci scalar and its first and second derivatives. Extracting the field equations we obtain an effective dark energy sector that consists of both extra scalar degrees of freedom, and we determine various observables. We analyze two specific models and we obtain a cosmological behavior in agreement with observations, i.e. transition from matter to dark energy era, with the onset of cosmic acceleration. Additionally, for a particular range of the model parameters, the equation-of-state parameter of the effective dark energy sector can exhibit the phantom-divide crossing. These features reveal the capabilities of these theories, since they arise solely from the novel, higher-derivative terms.

We investigate the Einstein static universe (ESU) and the emergent universe scenario in the framework of Ho\v{r}ava-Lifshitz-like $F(R)$ gravity. We first perform a dynamical analysis in the phase space, and amongst others we show that a spatially open universe filled with matter satisfying the strong energy condition can exhibit a stable static phase. Additionally, we examine the behavior of the scenario under scalar perturbations and extract the conditions under which it is free of perturbative instabilities, showing that the obtained background ESU solutions are free of such instabilities. However, in order for the Einstein static universe to give rise to the emergent universe scenario we need to have an exotic matter sector that can lead the universe to depart from the stable static state and enter into its usual expanding thermal history.

The Schwarzschild-de Sitter (SdS) solution exists in the large majority of modified gravity theories, as expected, and in particular the effective cosmological constant is determined by the specific parameters of the given theory. We explore the possibility to use future extended radio-tracking data from the currently ongoing New Horizons mission in the outskirts peripheries of the Solar System, at about 40 au, in order to constrain this effective cosmological constant, and thus to impose constrain on each scenario's parameters. We investigate some of the recently most studied modified gravities, namely $f(R)$ and $f(T)$ theories, dRGT massive gravity, and Ho\v{r}ava-Lifshitz gravity, and we show that New Horizons mission may bring an improvement of one-two orders of magnitude with respect to the present bounds from planetary orbital dynamics.

We show that the well-known problem of frame dependence and violation of local Lorentz invariance in the usual formulation of $f(T)$ gravity is a consequence of neglecting the role of spin connection. We re-formulate $f(T)$ gravity starting, instead of the "pure-tetrad" teleparallel gravity, from the covariant teleparallel gravity, using both the tetrad and the spin connection as dynamical variables, resulting in the fully covariant, consistent, and frame-independent, version of $f(T)$ gravity, which does not suffer from the notorious problems of the usual, pure-tetrad, $f(T)$ theory. We present the method to extract solutions for the most physically important cases, such as the Minkowski, the FRW and the spherically-symmetric ones. We show that in the covariant $f(T)$ gravity we are allowed to use an arbitrary tetrad in an arbitrary coordinate system along with the corresponding spin connection, resulting always to the same physically relevant field equations.

We investigate the cosmological applications of a bi-scalar modified gravity that exhibits partial conformal invariance, which could become full conformal invariance in the absence of the usual Einstein-Hilbert term and introducing additionally either the Weyl derivative or properly rescaled fields. Such a theory is constructed by considering the action of a non-minimally conformally-coupled scalar field, and adding a second scalar allowing for a nonminimal derivative coupling with the Einstein tensor and the energy-momentum tensor of the first field. At a cosmological framework we obtain an effective dark-energy sector constituted from both scalars. In the absence of an explicit matter sector we extract analytical solutions, which for some parameter regions correspond to an effective matter era and/or to an effective radiation era, thus the two scalars give rise to "mimetic dark matter" or to "dark radiation" respectively. In the case where an explicit matter sector is included we obtain a cosmological evolution in agreement with observations, that is a transition from matter to dark energy era, with the onset of cosmic acceleration. Furthermore, for particular parameter regions, the effective dark-energy equation of state can transit to the phantom regime at late times. These behaviours reveal the capabilities of the theory, since they arise purely from the novel, bi-scalar construction and the involved couplings between the two fields.

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 this work, we perform a detailed dynamical analysis for the cosmological applications of a nonminimal torsion-matter coupled gravity. Two alternative formalisms are proposed, which enable one to choose between the easier approach for a given problem, and furthermore, we analyze six specific models. In general, we extract fixed points corresponding either to dark-matter dominated, scaling decelerated solutions, or to dark-energy dominated accelerated solutions. Additionally, we find that there is a small parameter region in which the model can experience the transition from the matter epoch to a dark-energy era. These features are in agreement with the observed universe evolution, and make the theory a successful candidate for the description of Nature.

We investigate the cosmological behavior in a universe governed by time asymmetric extensions of general relativity, which is a novel modified gravity based on the addition of new, time-asymmetric, terms on the Hamiltonian framework, in a way that the algebra of constraints and local physics remain unchanged. Nevertheless, at cosmological scales these new terms can have significant effects that can alter the universe evolution, both at early and late times, and the freedom in the choice of the involved modification function makes the scenario able to produce a huge class of cosmological behaviors. For basic ansatzes of modification, we perform a detailed dynamical analysis, extracting the stable late-time solutions. Amongst others, we find that the universe can result in dark-energy dominated, accelerating solutions, even in the absence of an explicit cosmological constant, in which the dark energy can be quintessence-like, phantom-like, or behave as an effective cosmological constant. Moreover, it can result to matter-domination, or to a Big Rip, or experience the sequence from matter to dark energy domination. Additionally, in the case of closed curvature, the universe may experience a cosmological bounce or turnaround, or even cyclic behavior. Finally, these scenarios can easily satisfy the observational and phenomenological requirements. Hence, time asymmetric cosmology can be a good candidate for the description of the universe.

We present a series of theoretical constraints on the potentially viable inflation models that might yield a blue spectrum for primordial tensor perturbations. By performing a detailed dynamical analysis we show that, while there exists such possibility, the corresponding phase space is strongly bounded. Our result implies that, in order to achieve a blue tilt for inflationary tensor perturbations, one may either construct a non-canonical inflation model delicately, or study the generation of primordial tensor modes beyond the standard scenario of single slow-roll field.

We investigate the realization of two bouncing paradigms, namely of the superbounce and the loop quantum cosmological ekpyrosis, in the framework of various modified gravities. In particular, we focus on the $F(R)$, $F(G)$ and $F(T)$ gravities, and we reconstruct their specific subclasses which lead to such universe evolutions. These subclasses constitute from power laws, polynomials, or hypergeometric ansatzes, which can be approximated by power laws. The qualitative similarity of different effective gravities which realize the above two bouncing cosmologies, indicates to some universality lying behind such a bounce. Finally, performing a linear perturbation analysis, we show that the obtained solutions are conditionally or fully stable.

We demonstrate that finite time singularities of Type IV can be consistently incorporated in the Universe's cosmological evolution, either appearing in the inflationary era, or in the late-time regime. While using only one scalar field instabilities can in principle occur at the time of the phantom-divide crossing, when two fields are involved we are able to avoid such instabilities. Additionally, the two-field scalar-tensor theories prove to be able to offer a plethora of possible viable cosmological scenarios, at which various types of cosmological singularities can be realized. Amongst others, it is possible to describe inflation with the appearance of a Type IV singularity, and phantom late-time acceleration which ends in a Big Rip. Finally, for completeness, we also present the Type IV realization in the context of suitably reconstructed $F(R)$ gravity.

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 extract constraints on the transition redshift $z_{tr}$, determining the onset of cosmic acceleration, predicted by an effective cosmographic construction, in the framework of $f(T)$ gravity. In particular, employing cosmography we obtain bounds on the viable $f(T)$ forms and their derivatives. Since this procedure is model independent, as long as the scalar curvature is fixed, we are able to determine intervals for $z_{tr}$. In this way we guarantee that the Solar-System constraints are preserved and moreover we extract bounds on the transition time and the free parameters of the scenario. We find that the transition redshifts predicted by $f(T)$ cosmology, although compatible with the standard $\Lambda$CDM predictions, are slightly smaller. Finally, in order to obtain observational constraints on $f(T)$ cosmology, we perform a Monte Carlo fitting using supernova data, involving the most recent union 2.1 data set.

Recently a new Lagrangian framework was introduced to describe interactions between scalar fields and relativistic perfect fluids. This allows two consistent generalizations of coupled quintessence models: non-vanishing pressures and a new type of derivative interaction. Here the implications of these to the formation of cosmological large-scale structure are uncovered at the linear order. The full perturbation equations in the two cases are derived in a unified formalism and their Newtonian, quasi-static limit is studied analytically. Requiring the absence of an effective sound speed for the coupled dark matter fluid restricts the Lagrangian to be a linear function of the matter number density. This still leaves new potentially viable classes of both algebraically and derivatively interacting models wherein the coupling may impact the background expansion dynamics and imprint signatures into the large-scale structure.

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 present an extension of $f(T)$ gravity, allowing for a general coupling of the torsion scalar $T$ with the trace of the matter energy-momentum tensor $\mathcal{T}$. The resulting $f(T,\mathcal{T})$ theory is a new modified gravity, since it is different from all the existing torsion or curvature based constructions. Applied to a cosmological framework, it leads to interesting phenomenology. In particular, one can obtain a unified description of the initial inflationary phase, the subsequent non-accelerating, matter-dominated expansion, and then the transition to a late-time accelerating phase. Additionally, the effective dark energy sector can be quintessence or phantom-like, or exhibit the phantom-divide crossing during the evolution. Moreover, in the far future the universe results either to a de Sitter exponential expansion, or to eternal power-law accelerated expansions. Finally, a detailed study of the scalar perturbations at the linear level reveals that $f(T,\mathcal{T})$ cosmology can be free of ghosts and instabilities for a wide class of ansatzes and model parameters.

We investigate the cosmological implications of the recently constructed 5-dimensional braneworld cosmology with gravitating Nambu-Goto matching conditions. Inserting both matter and radiation sectors, we first extract the analytical cosmological solutions. Additionally, we use observational data from Type Ia Supernovae (SNIa) and Baryon Acoustic Oscillations (BAO), along with requirements of Big Bang Nucleosynthesis (BBN), in order to impose constraints on the parameters of the model. We find that the scenario at hand is in very good agreement with observations, and thus a small departure from the standard Randall-Sundrum scenario is allowed.

We investigate the cosmological applications of $F(T,T_G)$ gravity, which is a novel modified gravitational theory based on the torsion invariant $T$ and the teleparallel equivalent of the Gauss-Bonnet term $T_{G}$. $F(T,T_{G})$ gravity differs from both $F(T)$ theories as well as from $F(R,G)$ class of curvature modified gravity, and thus its corresponding cosmology proves to be very interesting. In particular, it provides a unified description of the cosmological history from early-times inflation to late-times self-acceleration, without the inclusion of a cosmological constant. Moreover, the dark energy equation-of-state parameter can be quintessence or phantom-like, or experience the phantom-divide crossing, depending on the parameters of the model.

The theory of nonlinear massive gravity can be extended into the F(R) form as developed in Phys.Rev.D90, 064051 (2014). Being free of the Boulware-Deser ghost, such a construction has the additional advantage of exhibiting no linear instabilities around a cosmological background. We investigate various cosmological evolutions of a universe governed by this generalized massive gravitational theory. Specifically, under the Starobinsky ansantz, this model provides a unified description of the cosmological history, from early-time inflation to late-time self-acceleration. Moreover, under viable F(R) forms, the scenario leads to a very interesting dark-energy phenomenology, including the realization of the quintom scenario without any pathology. Finally, we provide a detailed analysis of the cosmological perturbations at linear order, as well as the Hamiltonian constraint analysis, in order to examine the physical degrees of freedom.

We propose a nonlinear massive gravitational theory which includes $F(R)$ modifications. This construction inherits the benefits of the de Rham-Gabadadze-Tolley model and is free of the Boulware-Deser ghost due to the existence of a Hamiltonian constraint accompanied by a nontrivial secondary one. The scalar perturbations in a cosmological background can be stabilized at the linear level for a wide class of the $F(R)$ models. The linear scalar mode arisen from the $F(R)$ sector can absorb the nonlinear longitudinal graviton, and hence, our scenario demonstrates the possibility of a gravitational Goldstone theorem. Finally, due to the combined contribution of the $F(R)$ and graviton-mass sectors, the proposed theory allows for a large class of cosmological evolutions, such as the simultaneous and unified description of inflation and late-time acceleration.

Inspired by the teleparallel formulation of General Relativity, whose Lagrangian is the torsion invariant T, we have constructed the teleparallel equivalent of Gauss-Bonnet gravity in arbitrary dimensions. Without imposing the Weitzenbock connection, we have extracted the torsion invariant T_G, equivalent (up to boundary terms) to the Gauss-Bonnet term G. T_G is constructed by the vielbein and the connection, it contains quartic powers of the torsion tensor, it is diffeomorphism and Lorentz invariant, and in four dimensions it reduces to a topological invariant as expected. Imposing the Weitzenbock connection, T_G depends only on the vielbein, and this allows us to consider a novel class of modified gravity theories based on F(T,T_G), which is not spanned by the class of F(T) theories, nor by the F(R,G) class of curvature modified gravity. Finally, varying the action we extract the equations of motion for F(T,T_G) gravity.

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.

We confront various nonsingular bouncing cosmologies with the recently released BICEP2 data and investigate the observational constraints on their parameter space. In particular, within the context of the effective field approach, we analyze the constraints on the matter bounce curvaton scenario with a light scalar field, and the new matter bounce cosmology model in which the universe successively experiences a period of matter contraction and an ekpyrotic phase. Additionally, we consider three nonsingular bouncing cosmologies obtained in the framework of modified gravity theories, namely the Ho\v{r}ava-Lifshitz bounce model, the $f(T)$ bounce model, and loop quantum cosmology.

The $f(T,T_G)$ class of gravitational modification, based on the quadratic torsion scalar $T$, as well as on the new quartic torsion scalar $T_G$ which is the teleparallel equivalent of the Gauss-Bonnet term, is a novel theory, different from both $f(T)$ and $f(R,G)$ ones. We perform a detailed dynamical analysis of a spatially flat universe governed by the simplest non-trivial model of $f(T,T_G)$ gravity which does not introduce a new mass scale. We find that the universe can result in dark-energy dominated, quintessence-like, cosmological-constant-like or phantom-like solutions, according to the parameter choices. Additionally, it may result to a dark energy - dark matter scaling solution, and thus it can alleviate the coincidence problem. Finally, the analysis "at infinity" reveals that the universe may exhibit future, past, or intermediate singularities depending on the parameters.

We construct an extension of f(T) gravity with the inclusion of a non-minimal torsion-matter coupling in the action. The resulting theory is a novel gravitational modification, since it is different from both f(T) gravity, as well as from the non-minimal curvature-matter-coupled theory. The cosmological application of this new theory proves to be very interesting. In particular, we obtain an effective dark energy sector whose equation-of-state parameter can be quintessence or phantom-like, or exhibit the phantom-divide crossing, while for a large range of the model parameters the Universe results in a de Sitter, dark-energy-dominated, accelerating phase. Additionally, we can obtain early-time inflationary solutions too, and thus provide a unified description of the cosmological history.

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.

Spherically symmetric solutions for f(T) gravity models are derived by the so called Noether Symmetry Approach. First, we present a full set of Noether symmetries for some minisuperspace models. Then, we compute analytical solutions and find that spherically symmetric solutions in f(T) gravity can be recast in terms of Schwarzschild-like solutions modified by a distortion function depending on a characteristic radius. The obtained solutions are more general than those obtained by the usual solution methods.

The self-gravitating gas in the Newtonian limit is studied in the presence of dark energy with a linear and constant equation of state. Entropy extremization associates to the isothermal Boltzmann distribution an effective density that includes `dark energy particles', which either strengthen or weaken mutual gravitational attraction, in case of quintessence or phantom dark energy, respectively, that satisfy a linear equation of state. Stability is studied for microcanonical (fixed energy) and canonical (fixed temperature) ensembles. Compared to the previously studied cosmological constant case, in the present work it is found that quintessence increases, while phantom dark energy decreases the instability domain under gravitational collapse. Thus, structures are more easily formed in a quintessence rather than in a phantom dominated Universe. Assuming that galaxy clusters are spherical, nearly isothermal and in hydrostatic equilibrium we find that dark energy with a linear and constant equation of state, for fixed radius, mass and temperature, steepens their total density profile. In case of a cosmological constant, this effect accounts for a 1.5% increase in the density contrast, that is the center to edge density ratio of the cluster. We also propose a method to constrain phantom dark energy.

We perform a detailed dynamical analysis of anisotropic scalar-field cosmologies, and in particular of the most significant Kantowski-Sachs, Locally Rotationally Symmetric (LRS) Bianchi I and LRS Bianchi III cases. We follow the new and powerful method of $f$-devisers, which allows us to perform the whole analysis for a wide range of potentials. Thus, one can just substitute the specific potential form in the final results and obtain the corresponding behavior, without the need of new calculations. We find a very rich behavior, and amongst others the universe can result in isotropized solutions with observables in agreement with observations, such as de Sitter, quintessence-like, or stiff-dark energy solutions. In particular, all expanding, accelerating, stable attractors are isotropic. Additionally, we prove that as long as matter obeys the null energy condition, bounce behavior is impossible. Finally, applying the general results to the well-studied exponential and power-law potentials, we find that some of the general stable solutions disappear. This feature may be an indication that such simple potentials might restrict the dynamics in scalar-field cosmology, opening the way to the introduction of more complicated ones.

We investigate the cosmology of interacting spin-2 particles, formulating the multi-gravitational theory in terms of vierbeins and without imposing any Deser-van Nieuwenhuizen-like constraint. The resulting multi-vierbein theory represents a wider class of gravitational theories if compared to the corresponding multi-metric models. Moreover, as opposed to its metric counterpart which in general seems to contain ghosts, it has already been proved to be ghost-free. We outline a discussion about the possible matter couplings and we focus on the study of cosmological scenarios in the case of three and four interacting vierbeins. We find rich behavior, including de Sitter solutions with an effective cosmological constant arising from the multi-vierbein interaction, dark-energy solutions and nonsingular bouncing solutions.

We investigate the cosmological predictions of several $f(T)$ models, with up to two parameters, at both the background and the perturbation levels. Using current cosmological observations (geometric supernovae type Ia, cosmic microwave background and baryonic acoustic oscillation and dynamical growth data) we impose constraints on the distortion parameter, which quantifies the deviation of these models from the concordance $\Lambda$ cosmology at the background level. In addition we constrain the growth index $\gamma$ predicted in the context of these models using the latest perturbation growth data in the context of three parametrizations for $\gamma$. The evolution of the best fit effective Newton constant, which incorporates the $f(T)$-gravity effects, is also obtained along with the corresponding $1\sigma$ error regions. We show that all the viable parameter sectors of the $f(T)$ gravity models considered practically reduce these models to $\Lambda$CDM. Thus, the degrees of freedom that open up to $\Lambda$CDM in the context of $f(T)$ gravity models are not utilized by the cosmological data leading to an overall disfavor of these models.

In this article we present the cosmological equivalence between the relativistic Finsler-Randers cosmology, with dark energy and modified gravity constructions, at the background level. Starting from a small deviation from the quadraticity of the Riemannian geometry, through which the local structure of General Relativity is modified and the curvature theory is extended, we extract the modified Friedmann equation. The corresponding extended Finsler-Randers cosmology is very interesting, and it can mimic dark-energy and modified gravity, describing a large class of scale-factor evolutions, from inflation to late-time acceleration, including the phantom regime. In this respect, the non-trivial universe evolution is not attributed to a new scalar field, or to gravitational modification, but it arises from the modification of the geometry itself.

We perform a combined perturbation and observational investigation of the scenario of non-minimal derivative coupling between a scalar field and curvature. First we extract the necessary condition that ensures the absence of instabilities, which is fulfilled more sufficiently for smaller coupling values. Then using Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB) observations, we show that, contrary to its significant effects on inflation, the non-minimal derivative coupling term has a negligible effect on the universe acceleration, since it is driven solely by the usual scalar-field potential. Therefore, the scenario can provide a unified picture of early and late time cosmology, with the non-minimal derivative coupling term responsible for inflation, and the usual potential responsible for late-time acceleration. Additionally, the fact that the necessary coupling term does not need to be large, improves the model behavior against instabilities.

The transient acceleration which the current cosmic acceleration is not eternal is possible by introducing the interaction between dark matter and dark energy. If the energy transfer is from dark energy to dark matter, then it is possible to realize the transient acceleration. We study the possibility of transient acceleration by considering two time-dependent phenomenological interaction forms so that the energy transfer increases as the universe evolves. Starting from a simple and extending to a more complicated ansatz, we obtain analytical expressions for the evolutions of the deceleration and the various energy density parameters. We find the ranges of the parameters in the models for a transient acceleration.

We investigate non-singular bounce realizations in the framework of ghost-free generalized Galileon cosmology, which furthermore can be free of the anisotropy problem. Considering an Ekpyrotic-like potential we can obtain a total Equation-of-State (EoS) larger than one in the contracting phase, which is necessary for the evolution to be stable against small anisotropic fluctuations. Since such a large EoS forbids the Galileon field to generate the desired form of perturbations, we additionally introduce the curvaton field which can in general produce the observed nearly scale-invariant spectrum. In particular, we provide approximate analytical and exact semi-analytical expressions under which the bouncing scenario is consistent with observations. Finally, the combined Galileon-curvaton system is free of the Big-Rip after the bounce.

The Nine-Year WMAP results combined with other cosmological data seem to indicate an enhanced favor for the phantom regime, comparing to previous analyses. This behavior, unless reversed by future observational data, suggests to consider the phantom regime more thoroughly. In this work we provide three modified gravitational scenarios in which we obtain the phantom realization without the appearance of ghosts degrees of freedom, which plague the naive approaches on the subject, namely the Brans-Dicke type gravity, the scalar-Einstein-Gauss-Bonnet gravity, and the $F(R)$ gravity, which are moreover free of perturbative instabilities. The phantom regime seems to favor the gravitational modification instead of the universe-content alteration.

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.

We perform a detailed dynamical analysis of various cosmological scenarios in extended (varying-mass) nonlinear massive gravity. Due to the enhanced freedom in choosing the involved free functions, this cosmological paradigm allows for a huge variety of solutions that can attract the universe at late times, comparing to scalar-field cosmology or usual nonlinear massive gravity. Amongst others, it accepts quintessence, phantom, or cosmological-constant-like late-time solutions, which moreover can alleviate the coincidence problem. These features seem to be general and non-sensitive to the imposed ansantzes and model parameters, and thus extended nonlinear massive gravity can be a good candidate for the description of nature.

We investigate a (super-)renormalizable and ghost-free theory of gravity, showing that under a natural (exponential) ansatz of the form factor and a suitable truncation it can give rise to the Starobinsky inflationary theory in cosmological frameworks, and thus offering a theoretical justification of its origin. We study the corresponding inflationary evolution and we examine the generation of curvature perturbations, adapting the $f(R)$-like equations in a symmetry-reduced FLRW metric. Furthermore, we analyze how the ultraviolet regime of a simply renormalizable and unitary theory of gravity is also compatible with the Starobinsky action, and hence we show that such a theory could account for an inflationary phase of the Universe in the ultraviolet regime.

We perform a detailed dynamical analysis of generalized Galileon cosmology, incorporating also the requirements of ghost and instabilities absence. We find that there are not any new stable late-time solutions apart from those of standard quintessence. Furthermore, depending on the model parameters the Galileons may survive at late times or they may completely disappear by the dynamics, however the corresponding observables are always independent of the Galileon terms, determined only by the usual action terms. Thus, although the Galileons can play an important role at inflationary or at recent times, in the future, when the universe will asymptotically reach its stable state, they will not have any effect on its evolution.

We investigate the cosmological evolution in a universe governed by the extended, varying-mass, nonlinear massive gravity, in which the graviton mass is promoted to a scalar-field. We find that the dynamics, both in flat and open universe, can lead the varying graviton mass to zero at late times, offering a natural explanation for its hugely-constrained observed value. Despite the limit of the scenario towards standard quintessence, at early and intermediate times it gives rise to an effective dark energy sector of a dynamical nature, which can also lie in the phantom regime, from which it always exits naturally, escaping a Big-Rip. Interestingly enough, although the motivation of massive gravity is to obtain an IR modification, its varying-mass extension in cosmological frameworks leads to early and intermediate times modification instead.

We find static charged black hole solutions in nonlinear massive gravity. In the parameter space of two gravitational potential parameters $(\alpha, \beta)$ we show that below the Compton wavelength the black hole solutions reduce to that of Reissner-Nordstr\"om via the Vainshtein mechanism in the weak field limit. In the simplest case with $\alpha=\beta=0$ the solution exhibits the vDVZ discontinuity but ordinary General Relativity is recovered deep inside the horizon due to the existence of electric charge. For $\alpha\neq0$ and $\beta=0$, the post-Newtonian parameter of the charged black hole evolves to that of General Relativity via the Vainshtein mechanism within a macroscopic distance; however, a logarithmic correction to the metric factor of the time coordinate is obtained. When $\alpha$ and $\beta$ are both nonzero, there exist two branches of solutions depending on the positivity of $\beta$. When $\beta<0$, the strong coupling of the scalar graviton weakens gravity at distances smaller than the Vainshtein radius. However, when $\beta>0$ the metric factors exhibit only small corrections compared to the solutions obtained in General Relativity, and under a particular choice of $\beta=\alpha^2/6$ the standard Reissner-Nordstr\"om-de Sitter solution is recovered.

We investigate non-singular bounce and cyclic cosmological evolutions in a universe governed by the extended nonlinear massive gravity, in which the graviton mass is promoted to a scalar-field potential. The extra freedom of the theory can lead to certain energy conditions violations and drive cyclicity with two different mechanisms: either with a suitably chosen scalar-field potential under a given Stuckelberg-scalar function, or with a suitably chosen Stuckelberg-scalar function under a given scalar-field potential. Our analysis shows that extended nonlinear massive gravity can alter significantly the evolution of the universe at both early and late times.

We use recent observations from solar system orbital motions in order to constrain f(T) gravity. In particular, imposing a quadratic f(T) correction to the linear-in-T form, which is a good approximation for every realistic case, we extract the spherical solutions of the theory. Using these spherical solutions to describe the Sun's gravitational field, we use recently determined supplementary advances of planetary perihelia, to infer upper bounds on the allowed f(T) corrections. We find that the maximal allowed divergence of the gravitational potential in f(T) gravity from that in the teleparallel equivalent of General Relativity is of the order of 6.2 \times 10^{-10}, in the applicability region of our analysis. This is much smaller than the corresponding (significantly small too) divergence that is predicted from cosmological observations, as expected. Such a tiny allowed divergence from the linear form should be taken into account in f(T) model building.

In the N=1 four-dimensional new-minimal supergravity framework, we supersymmetrise the coupling of the scalar kinetic term to the Einstein tensor. This coupling, although introduces a non-minimal derivative interaction of curvature to matter, it does not introduce harmful higher-derivatives. For this construction, we employ off-shell chiral and real linear multiplets. Physical scalars are accommodated in the chiral multiplet whereas curvature resides in a linear one.

We examine various entropic inflation scenarios, under the light of eternality. After describing the inflation realization and the normal condition for inflation to last at the background level, we investigate the conditions for eternal inflation with the effect of thermal fluctuations produced from standard radiation and from the holographic screen. Furthermore, we incorporate stochastic quantum fluctuations through a phenomenological, Langevin analysis, studying whether they can affect the inflation eternality. In single-holographic-screen scenarios eternality can be easily obtained, while in double-screen considerations inflation is eternal only in the high-energy regime. Thus, from the cosmological point of view, one should take these into account before he can consider entropic gravity as a candidate for the description of nature. However, form the string theory point of view, inflation eternality may form the background for the "Landscape" of string/M theory vacua, leading to new perspectives in entropy gravity.

We use data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB) observations to constrain the recently proposed teleparallel dark energy scenario based on the teleparallel equivalent of General Relativity, in which one adds a canonical scalar field, allowing also for a nonminimal coupling with gravity. Using the power-law, the exponential and the inverse hyperbolic cosine potential ansatzes, we show that the scenario is compatible with observations. In particular, the data favor a nonminimal coupling, and although the scalar field is canonical the model can describe both the quintessence and phantom regimes.

We investigate cosmological scenarios of generalized Chaplygin gas in a universe governed by Horava-Lifshitz gravity. We consider both the detailed and non-detailed balance versions of the gravitational background, and we include the baryonic matter and radiation sectors. We use observational data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis (BBN), to constrain the parameters of the model, and we provide the corresponding likelihood contours. We deduce that the present scenario is compatible with observations. Additionally, examining the evolution of the total equation-of-state parameter, we find in a unified way the succession of the radiation, matter, and dark energy epochs, consistently with the thermal history of the universe.

We investigate phantom cosmology in which the scale factor is a power law, and we use cosmological observations from Cosmic Microwave Background (CMB), Baryon Acoustic Oscillations (BAO) and observational Hubble data, in order to impose complete constraints on the model parameters. We find that the power-law exponent is $\beta\approx-6.51^{+0.24}_{-0.25}$, while the Big Rip is realized at $t_s\approx104.5^{+1.9}_{-2.0}$ Gyr, in 1$\sigma$ confidence level. Providing late-time asymptotic expressions, we find that the dark-energy equation-of-state parameter at the Big Rip remains finite and equal to $w_{DE}\approx -1.153$, with the dark-energy density and pressure diverging. Finally, we reconstruct the phantom potential.

Using the "teleparallel" equivalent of General Relativity as the gravitational sector, which is based on torsion instead of curvature, we add a canonical scalar field, allowing for a nonminimal coupling with gravity. Although the minimal case is completely equivalent to standard quintessence, the nonminimal scenario has a richer structure, exhibiting quintessence-like or phantom-like behavior, or experiencing the phantom-divide crossing. The richer structure is manifested in the absence of a conformal transformation to an equivalent minimally-coupled model.

We show that the f(T) gravitational paradigm, in which gravity is described by an arbitrary function of the torsion scalar, can provide a mechanism for realizing bouncing cosmologies, thereby avoiding the Big Bang singularity. After constructing the simplest version of an f(T) matter bounce, we investigate the scalar and tensor modes of cosmological perturbations. Our results show that metric perturbations in the scalar sector lead to a background-dependent sound speed, which is a distinguishable feature from Einstein gravity. Additionally, we obtain a scale-invariant primordial power spectrum, which is consistent with cosmological observations, but suffers from the problem of a large tensor-to-scalar ratio. However, this can be avoided by introducing extra fields, such as a matter bounce curvaton.

In this article we review recent developments of cyclic cosmology. A typical non-singular cyclic model within General Relativity requires a non-conventional fluid with negative effective energy density, in order to cancel the matter component and lead to a non-singular bounce. However, the existence of such a non-conventional fluid usually leads to quantum instabilities and makes the theory ill-defined. In the present work we follow the alternative way, obtaining two scenarios of non-singular cyclic cosmological evolutions in the context of gravitational theories beyond General Relativity. The degrees-of-freedom examination reveals that these two models are free of the Phantom Menace. Our analysis illustrates that, if cyclic cosmology describes overall universe, a theory of gravity beyond Einstein may be expected.

We investigate f(R)-gravity models performing the ADM-slicing of standard General Relativity. We extract the static, spherically-symmetric vacuum solutions in the general case, which correspond to either Schwarzschild de-Sitter or Schwarzschild anti-de-Sitter ones. Additionally, we study the cosmological evolution of a homogeneous and isotropic universe, which is governed by an algebraic and not a differential equation. We show that the universe admits solutions corresponding to acceleration at late cosmological epochs, without the need of fine-tuning the model-parameters or the initial conditions.

We construct general anisotropic cosmological scenarios governed by an $f(R)$ gravitational sector. Focusing then on Kantowski-Sachs geometries in the case of $R^n$-gravity, and modelling the matter content as a perfect fluid, we perform a detailed phase-space analysis. We find that at late times the universe can result to a state of accelerating expansion, and additionally, for a particular $n$-range ($2<n<3$) it exhibits phantom behavior. Furthermore, isotropization has been achieved independently of the initial anisotropy degree, showing in a natural way why the observable universe is so homogeneous and isotropic, without relying on a cosmic no-hair theorem. Moreover, contracting solutions have also a large probability to be the late-time states of the universe. Finally, we can also obtain the realization of the cosmological bounce and turnaround, as well as of cyclic cosmology. These features indicate that anisotropic geometries in modified gravitational frameworks present radically different cosmological behaviors comparing to the simple isotropic scenarios.

We investigate thermal inflation in double-screen entropic cosmology. We find that its realization is general, resulting from the system evolution from non-equilibrium to equilibrium. Furthermore, going beyond the background evolution, we study the primordial curvature perturbations arising from the universe interior, as well as from the thermal fluctuations generated on the holographic screens. We show that the power spectrum is nearly scale-invariant with a red tilt, while the tensor-to-scalar ratio is in agreement with observations. Finally, we examine the non-Gaussianities of primordial curvature perturbations, and we find that a sizable value of the non-linearity parameter is possible due to holographic statistics on the outer screen.

We investigate cyclic and singularity-free evolutions in a universe governed by Lagrange-multiplier modified gravity, either in scalar-field cosmology, as well as in $f(R)$ one. In the scalar case, cyclicity can be induced by a suitably reconstructed simple potential, and the matter content of the universe can be successfully incorporated. In the case of $f(R)$-gravity, cyclicity can be induced by a suitable reconstructed second function $f_2(R)$ of a very simple form, however the matter evolution cannot be analytically handled. Furthermore, we study the evolution of cosmological perturbations for the two scenarios. For the scalar case the system possesses no wavelike modes due to a dust-like sound speed, while for the $f(R)$ case there exist an oscillation mode of perturbations which indicates a dynamical degree of freedom. Both scenarios allow for stable parameter spaces of cosmological perturbations through the bouncing point.

We review some general aspects of Horava-Lifshitz cosmology. Formulating it in its basic version, we extract the cosmological equations and we use observational data in order to constrain the parameters of the theory. Through a phase-space analysis we extract the late-time stable solutions, and we show that eternal expansion, and bouncing and cyclic behavior can arise naturally. Concerning the effective dark energy sector we show that it can describe the phantom phase without the use of a phantom field. However, performing a detailed perturbation analysis, we see that Horava-Lifshitz gravity in its basic version suffers from instabilities. Therefore, suitable generalizations are required in order for this novel theory to be a candidate for the description of nature.

We investigate f(T) cosmology in both the background, as well as in the perturbation level, and we present the general formalism for reconstructing the equivalent one-parameter family of f(T) models for any given dynamical dark energy scenario. Despite the completely indistinguishable background behavior, the perturbations break this degeneracy and the growth histories of all these models differ from one another. As an application we reconstruct the f(T) equivalent for quintessence, and we show that the deviation of the matter overdensity evolution is strong for small scales and weak for large scales, while it is negligible for large redshifts.

We investigate the validity of the generalized second law of thermodynamics in a universe governed by Horava-Lifshitz gravity. Under the equilibrium assumption, that is in the late-time cosmological regime, we calculate separately the entropy time-variation for the matter fluid and, using the modified entropy relation, that of the apparent horizon itself. We find that under detailed balance the generalized second law is generally valid for flat and closed geometry and it is conditionally valid for an open universe, while beyond detailed balance it is only conditionally valid for all curvatures. Furthermore, we also follow the effective approach showing that it can lead to misleading results. The non-complete validity of the generalized second law could either provide a suggestion for its different application, or act as an additional problematic feature of Horava-Lifshitz gravity.

We investigate the cosmological perturbations in f(T) gravity. Examining the pure gravitational perturbations in the scalar sector using a diagonal vierbien, we extract the corresponding dispersion relation, which provides a constraint on the f(T) ansatzes that lead to a theory free of instabilities. Additionally, upon inclusion of the matter perturbations, we derive the fully perturbed equations of motion, and we study the growth of matter overdensities. We show that f(T) gravity with f(T) constant coincides with General Relativity, both at the background as well as at the first-order perturbation level. Applying our formalism to the power-law model we find that on large subhorizon scales (O(100 Mpc) or larger), the evolution of matter overdensity will differ from LCDM cosmology. Finally, examining the linear perturbations of the vector and tensor sectors, we find that (for the standard choice of vierbein) f(T) gravity is free of massive gravitons.

We investigate the new agegraphic dark energy scenario in a universe governed by Horava-Lifshitz gravity. We consider both the detailed and non-detailed balanced version of the theory, we impose an arbitrary curvature, and we allow for an interaction between the matter and dark energy sectors. Extracting the differential equation for the evolution of the dark energy density parameter and performing an expansion of the dark energy equation-of-state parameter, we calculate its present and its low-redshift value as functions of the dark energy and curvature density parameters at present, of the Horava-Lifshitz running parameter $\lambda$, of the new agegraphic dark energy parameter $n$, and of the interaction coupling $b$. We find that $w_0=-0.82^{+0.08}_{-0.08}$ and $w_1=0.08^{+0.09}_{-0.07}$. Although this analysis indicates that the scenario can be compatible with observations, it does not enlighten the discussion about the possible conceptual and theoretical problems of Horava-Lifshitz gravity.

We investigate the validity of the generalized second law of thermodynamics, in the cosmological scenario where dark energy interacts with both dark matter and radiation. Calculating separately the entropy variation for each fluid component and for the apparent horizon itself, we show that the generalized second law is always and generally valid, independently of the specific interaction form, of the fluids equation-of-state parameters and of the background geometry.

We have investigated the necessary conditions that prevent phantom inflation from being eternal. Allowing additionally for a nonminimal coupling between the phantom field and gravity, we present the slow-climb requirements, perform an analysis of the fluctuations, and finally we extract the overall conditions that are necessary in order to prevent eternality. Furthermore, we verify our results by solving explicitly the cosmological equations in a simple example of an exponential potential, formulating the classical motion plus the stochastic effect of the fluctuations through Langevin equations. Our analysis shows that phantom inflation can be finite without the need of additional exotic mechanisms.

We examine an extension of the quintom scenario of dark energy, in which a canonical scalar field and a phantom field are coupled through a kinetic interaction. We perform a phase space analysis and show that the kinetic coupling gives rise to novel cosmological behaviour. In particular, we obtain both quintessence-like and phantom-like late-time solutions, as well as solutions that cross the phantom divide during the evolution of the universe.

We use observational data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis (BBN), to constrain the running parameter $\lambda$ of Horava-Lifshitz gravity, which determines the flow between the Ultra-Violet and the Infra-Red. We consider both the detailed and non-detailed balance versions of the gravitational sector, and we include the matter and radiation sectors. Allowing for variation of all the parameters of the theory, we construct the likelihood contours and we conclude that in $1\sigma$ confidence $\lambda$ is restricted to $|\lambda-1|\lesssim0.02$, while its best fit value is $|\lambda_{b.f}-1|\approx0.002$. Although this observational analysis restricts the running parameter $\lambda$ very close to its IR value 1, it does not enlighten the discussion about the theory's possible conceptual and theoretical problems.

We review the paradigm of quintom cosmology. This scenario is motivated by the observational indications that the equation of state of dark energy across the cosmological constant boundary is mildly favored, although the data are still far from being conclusive. As a theoretical setup we introduce a no-go theorem existing in quintom cosmology, and based on it we discuss the conditions for the equation of state of dark energy realizing the quintom scenario. The simplest quintom model can be achieved by introducing two scalar fields with one being quintessence and the other phantom. Based on the double-field quintom model we perform a detailed analysis of dark energy perturbations and we discuss their effects on current observations. This type of scenarios usually suffer from a manifest problem due to the existence of a ghost degree of freedom, and thus we review various alternative realizations of the quintom paradigm. The developments in particle physics and string theory provide potential clues indicating that a quintom scenario may be obtained from scalar systems with higher derivative terms, as well as from non-scalar systems. Additionally, we construct a quintom realization in the framework of braneworld cosmology, where the cosmic acceleration and the phantom divide crossing result from the combined effects of the field evolution on the brane and the competition between four and five dimensional gravity. Finally, we study the outsets and fates of a universe in quintom cosmology. In a scenario with null energy condition violation one may obtain a bouncing solution at early times and therefore avoid the Big Bang singularity. Furthermore, if this occurs periodically, we obtain a realization of an oscillating universe. Lastly, we comment on several open issues in quintom cosmology and their connection to future investigations.

We use observational data from Type Ia Supernovae (SN), Baryon Acoustic Oscillations (BAO), Cosmic Microwave Background (CMB) and observational Hubble data (OHD), and the Markov Chain Monte Carlo (MCMC) method, to constrain the cosmological scenario of holographic dark energy with varying gravitational constant. We consider both flat and non-flat background geometry, and we present the corresponding constraints and contour-plots of the model parameters. We conclude that the scenario is compatible with observations. In 1$\sigma$ we find $\Omega_{\Lambda0}=0.72^{+0.03}_{-0.03}$, $\Omega_{k0}=-0.0013^{+0.0130}_{-0.0040}$, $c=0.80^{+0.19}_{-0.14}$ and $\Delta_G\equiv G'/G=-0.0025^{+0.0080}_{-0.0050}$, while for the present value of the dark energy equation-of-state parameter we obtain $w_0=-1.04^{+0.15}_{-0.20}$.

We formulate Horava-Lifshitz cosmology with an additional scalar field that leads to an effective dark energy sector. We find that, due to the inherited features from the gravitational background, Horava-Lifshitz dark energy naturally presents very interesting behaviors, possessing a varying equation-of-state parameter, exhibiting phantom behavior and allowing for a realization of the phantom divide crossing. In addition, Horava-Lifshitz dark energy guarantees for a bounce at small scale factors and it may trigger the turnaround at large scale factors, leading naturally to cyclic cosmology.

We investigate the scalar and tensor perturbations in Horava gravity, with and without detailed balance, around a flat background. Once both types of perturbations are taken into account, it is revealed that the theory is plagued by ghost-like scalar instabilities in the range of parameters which would render it power-counting renormalizable, that cannot be overcome by simple tricks such as analytic continuation. Implementing a consistent flow between the UV and IR limits seems thus more challenging than initially presumed, regardless of whether the theory approaches General Relativity at low energies or not. Even in the phenomenologically viable parameter space, the tensor sector leads to additional potential problems, such as fine-tunings and super-luminal propagation.

We investigate cosmological scenarios with a non-minimal derivative coupling between the scalar field and the curvature, examining both the quintessence and the phantom cases in zero and constant potentials. In general, we find that the universe transits from one de Sitter solution to another, determined by the coupling parameter. Furthermore, according to the parameter choices and without the need for matter, we can obtain a Big Bang, an expanding universe with no beginning, a cosmological turnaround, an eternally contracting universe, a Big Crunch, a Big Rip avoidance and a cosmological bounce. This variety of behaviors reveals the capabilities of the present scenario.

We use observational data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis (BBN), to constrain the cosmological scenarios governed by Horava-Lifshitz gravity. We consider both the detailed and non-detailed balance versions of the gravitational sector, and we include the matter and radiation sectors. We conclude that the detailed-balance scenario cannot be ruled out from the observational point of view, however the corresponding likelihood contours impose tight constraints on the involved parameters. The scenario beyond detailed balance is compatible with observational data, and we present the corresponding stringent constraints and contour-plots of the parameters. Although this analysis indicates that Horava-Lifshitz cosmology can be compatible with observations, it does not enlighten the discussion about its possible conceptual and theoretical problems.

We survey the application of specific tools to distinguish amongst the wide variety of dark energy models that are nowadays under investigation. The first class of tools is more mathematical in character: the application of the theory of dynamical systems to select the better behaved models, with appropriate attractors in the past and future. The second class of tools is rather physical: the use of astrophysical observations to crack the degeneracy of classes of dark energy models. In this last case the observations related with structure formation are emphasized both in the linear and non-linear regimes. We exemplify several studies based on our research, such as quintom and quinstant dark energy ones. Quintom dark energy paradigm is a hybrid construction of quintessence and phantom fields, which does not suffer from fine-tuning problems associated to phantom field and additionally it preserves the scaling behavior of quintessence. Quintom dark energy is motivated on theoretical grounds as an explanation for the crossing of the phantom divide, i.e. the smooth crossing of the dark energy state equation parameter below the value -1. On the other hand, quinstant dark energy is considered to be formed by quintessence and a negative cosmological constant, the inclusion of this later component allows for a viable mechanism to halt acceleration. We comment that the quinstant dark energy scenario gives good predictions for structure formation in the linear regime, but fails to do that in the non-linear one, for redshifts larger than one. We comment that there might still be some degree of arbitrariness in the selection of the best dark energy models.

We present a novel dark energy candidate, based upon the existence and dynamics of Ghost D-branes in a warped compactification of type IIB string theory. Gp-branes cancel the combined BPS sectors of the Dp-branes, while they preserve the same supersymmetries. We show that this scenario can naturally lead to either quintessence or phantom-like behaviors, depending on the form of the involved potentials and brane tension.

We present a model of non-relativistic gravitational theory which is power-counting renormalizable in 3+1 dimensional spacetime. When applied to cosmology, the relativity-violation terms lead to a dark radiation component, which can give rise to a bounce if dark radiation possesses negative energy density. Additionally, we investigate a cyclic extension of the non-singular cosmology in which the universe undergoes contractions and expansions periodically. In both scenarios the background theory is well defined at the quantum level.

We perform a detailed phase-space analysis of Horava-Lifshitz cosmology, with and without the detailed-balance condition. Under detailed-balance we find that the universe can reach a bouncing-oscillatory state at late times, in which dark-energy, behaving as a simple cosmological constant, is dominant. In the case where the detailed-balance condition is relaxed, we find that the universe reaches an eternally expanding, dark-energy-dominated solution, with the oscillatory state preserving also a small probability. Although this analysis indicates that Horava-Lifshitz cosmology can be compatible with observations, it does not enlighten the discussion about its possible conceptual and theoretical problems.

We investigate, in a unified and general way, the thermodynamic properties of dark energy with an arbitrary, varying equation-of-state parameter w(a). We find that all quantities are well defined and regular for every w(a), including at the -1-crossing, with the temperature being negative in the phantom regime (w(a)<-1) and positive in the quintessence one (w(a)>-1). The density and entropy are always positive while the chemical potential can be arbitrary. At the -1-crossing, both temperature and chemical potential are zero. The temperature negativity can only be interpreted in the quantum framework. The regular behavior of all quantities at the -1-crossing, leads to the conclusion that such a crossing does not correspond to a phase transition, but rather to a smooth cross-over.

We investigate the holographic dark energy scenario with a varying gravitational constant, in flat and non-flat background geometry. We extract the exact differential equations determining the evolution of the dark energy density-parameter, which include $G$-variation correction terms. Performing a low-redshift expansion of the dark energy equation of state, we provide the involved parameters as functions of the current density parameters, of the holographic dark energy constant and of the $G$-variation.

We investigate quintessence and phantom dark energy scenarios, in which the scalar fields evolve in nearly flat potentials and are non-minimally coupled to gravity. We show that all such models converge to a common behavior and we provide the corresponding approximate analytical expressions for $w(\Omega_\phi)$ and $w(a)$. We find that non-minimal coupling leads to richer cosmological behavior comparing to its minimal counterpart. In addition, comparison with Baryon Acoustic Oscillation and latest Supernovae data reveals that agreement can be established more easily and with less strict constraints on the model parameters.

We investigate the restrictions on the equation-of-state parameter of phantom cosmology, due to the minimum quantum gravitational requirements. We find that for all the examined $w_\Lambda(z)$-parametrizations and for arbitrary phantom potentials and spatial curvature, the phantom equation-of-state parameter is not restricted at all. This is in radical contrast with the quintessence paradigm, and makes phantom cosmology more robust and capable of constituting the underlying mechanism for dark energy.

We investigate phantom models with power-law potentials and we extract the early-time, "tracker", solutions under the assumption of matter domination. Contrary to quintessence case, we find that energy positivity requires normal power-law potentials instead of inverse power-law ones, with the potential exponent being bounded by the quadratic form. In addition, we analytically present the general cosmological solution at intermediate times, that is at low redshifts, which is the period of the transition from matter to dark-energy domination. The comparison with the exact evolution, arising from numerical elaboration, shows that the tracker solution agrees with the later within 2% for redshifts z>1.5, while the intermediate solution is accurate within 2% up to $z\approx 0.5$.

We investigate several varying-mass dark-matter particle models in the framework of phantom cosmology. We examine whether there exist late-time cosmological solutions, corresponding to an accelerating universe and possessing dark energy and dark matter densities of the same order. Imposing exponential or power-law potentials and exponential or power-law mass dependence, we conclude that the coincidence problem cannot be solved or even alleviated. Thus, if dark energy is attributed to the phantom paradigm, varying-mass dark matter models cannot fulfill the basic requirement that led to their construction.

We perform a detailed phase-space analysis of various phantom cosmological models, where the dark energy sector interacts with the dark matter one. We examine whether there exist late-time scaling attractors, corresponding to an accelerating universe and possessing dark energy and dark matter densities of the same order. We find that all the examined models, although accepting stable late-time accelerated solutions, cannot alleviate the coincidence problem, unless one imposes a form of fine-tuning in the model parameters. It seems that interacting phantom cosmology cannot fulfill the basic requirement that led to its construction.

We investigate quintom evolution in power-law potentials. We extract the early-time, tracker solutions, under the assumption of matter domination. Additionally, we derive analytical solutions at intermediate times, that is at low redshifts, which is the period during the transition from matter to dark-energy domination. A comparison with exact evolution reveals that the tracker solutions are valid within 98% for z>1.5, while the intermediate ones are accurate within 98% up to z=0.5. Using these expressions we extract two new w-parametrizations, one in terms of the redshift and one in terms of the dark-energy density-parameter, and we present various quintom evolution sub-classes, including quintessence-like or phantom-like cases, realization of the -1-crossing, and non-monotonic w-evolution.

We examine dark energy models in which a quintessence or a phantom field, $\phi$, rolls near the vicinity of a local minimum or maximum, respectively, of its potential $V(\phi)$. Under the approximation that $(1/V)(dV/d\phi) \ll 1$, [although $(1/V)(d^2 V/d\phi^2)$ can be large], we derive a general expression for the equation of state parameter $w$ as a function of the scale factor for these models. The dynamics of the field depends on the value of $(1/V)(d^2 V/d\phi^2)$ near the extremum, which describes the potential curvature. For quintessence models, when $(1/V)(d^2 V/d\phi^2)<3/4$ at the potential minimum, the equation of state parameter $w(a)$ evolves monotonically, while for $(1/V)(d^2 V/d\phi^2)>3/4$, $w(a)$ has oscillatory behavior. For phantom fields, the dividing line between these two types of behavior is at $(1/V)(d^2 V/d\phi^2) = -3/4$. Our analytical expressions agree within 1% with the exact (numerically-derived) behavior, for all of the particular cases examined, for both quintessence and phantom fields. We present observational constraints on these models.

We investigate the physically accepted solutions of general Braneworld scenarios, scanning uniformly the associated parameter space. Without making any further assumptions we find that solutions which give "small" Hubble parameters on the physical brane, and therefore "small" effective cosmological constants on the 4D Universe, are far more probable than those with "large" ones. Eventually, their distribution tends to the $\delta$-function in the limit of continuous covering of the parameter space.

We examine quintom dark energy models, produced by the combined consideration of a canonical and a phantom field, with nearly flat potentials and dark energy equation-of-state parameter $w_{DE}$ close to -1. We find that all such models converge to a single expression for $w_{DE}(z)$, depending only on the initial field values and their derivatives. We show that this quintom paradigm allows for a description of the transition through -1 in the near cosmological past. In addition, we provide the necessary conditions for the determination of the direction of the -1-crossing.

We investigate general braneworld models, with a non-minimally coupled phantom bulk field and arbitrary brane and bulk matter contents. We show that the effective dark energy of the brane-universe acquires a dynamical nature, as a result of the non-minimal coupling which provides a mechanism for an indirect "bulk-brane interaction" through gravity. For late-time cosmological evolution and without resorting to special ansatzes or to specific areas of the parameter space, we show that the -1-crossing of its equation-of-state parameter is general and can be easily achieved. As an example we provide a simple, but sufficiently general, approximate analytical solution, that presents the crossing behavior.

We investigate canonical, phantom and quintom models, with the various fields being non-minimally coupled to gravity, in the framework of holographic dark energy. We classify them and we discuss their cosmological implications. In particular, we examine the present value of the dark energy equation-of-state parameter and the crossing through the phantom divide, and we extract the conditions for a future cosmological singularity. The combined scenarios are in agreement with observations and reveal interesting cosmological behaviors.

In the present work we investigate the cosmological implications of holographic dark energy density in the Gauss-Bonnet framework. By formulating independently the two cosmological scenarios, and by enforcing their simultaneous validity, we show that there is a correspondence between the holographic dark energy scenario in flat universe and the phantom dark energy model in the framework of Gauss-Bonnet theory with a potential. This correspondence leads consistently to an accelerating universe. However, in general one has not full freedom of constructing independently the two cosmological scenarios. Specific constraints must be imposed on the coupling with gravity and on the potential.

We investigate the full 5D dynamics of general braneworld models. Without making any further assumptions we show that cyclic behavior can arise naturally in a fraction of physically accepted solutions. The model does not require brane collisions, which in the stationary case remain fixed, and cyclicity takes place on the branes. We indicate that the cosmological constants play the central role for the realization of cyclic solutions and we show that its extremely small value on the observable universe makes the period of the cycles and the maximum scale factor astronomically large.

We investigate the phase-space structure of the quintom paradigm in the framework of a spatially flat, open, or closed isotropic and homogeneous universe. We examine the dynamical evolution under the assumption of late-time dark energy domination, without specifying the explicit quintom potential form. The obtained cosmological behavior is qualitatively different than that acquired from the single phantom model.

We apply the bulk holographic dark energy in general 5D two-brane models. We extract the Friedmann equation on the physical brane and we show that in the general moving-brane case the effective 4D holographic dark energy behaves as a quintom for a large parameter-space area of a simple solution subclass. We find that $w_\Lambda$ was larger than -1 in the past while its present value is $w_{\Lambda_0}\approx-1.05$, and the phantom bound $w_\Lambda=-1$ was crossed at $z_{p}\approx0.41$, a result in agreement with observations. Such a behavior arises naturally, without the inclusion of special fields or potential terms, but a fine-tuning between the 4D Planck mass and the brane tension has to be imposed.

We present a generalized version of holographic dark energy arguing that it must be considered in the maximally subspace of a cosmological model. In the context of brane cosmology it leads to a bulk holographic dark energy which transfers its holographic nature to the effective 4D dark energy. As an application we use a single-brane model and we show that in the low energy limit the behavior of the effective holographic dark energy coincides with that predicted by conventional 4D calculations. However, a finite bulk can lead to radically different results.

We demonstrate that energy can be transferred very rapidly from the bulk to the brane through parametric resonance. In a simple realization, we consider a massless bulk field that interacts only with fields localized on the brane. We assume an initial field configuration that has the form of a wave-packet moving towards the brane. During the reflection of the wave-packet by the brane the localized fields can be excited through parametric resonance. The mechanism is also applicable to bulk fields with a potential. The rapid energy transfer can have important cosmological and astrophysical implications.