We study the scenario of Kaniadakis horizon-entropy cosmology, which arises from the application of the gravity-thermodynamics conjecture using the Kaniadakis modified entropy. The resulting modified Friedmann equations contain extra terms that constitute an effective dark energy sector. We use data from cosmic chronometers, Type Ia supernova, H ii galaxies, strong lensing systems, and baryon acoustic oscillation observations, and we apply a Bayesian Markov chain Monte Carlo analysis to construct the likelihood contours for the model parameters. We find that the Kaniadakis parameter is constrained around 0, namely around the value where the standard Bekenstein–Hawking is recovered. Concerning the normalized Hubble parameter, we find |$h=0.708^{+0.012}_{-0.011}$|⁠, a result that is independently verified by applying the |$\mathbf {\mathbb {H}}0(z)$| diagnostic and, thus, we conclude that the scenario at hand can alleviate the H₀ tension problem. Regarding the transition redshift, the reconstruction of the cosmographic parameters gives |$z_{\rm T}=0.715^{+0.042}_{-0.041}$|⁠. Furthermore, we apply the Akaike, Bayesian, and deviance information criteria, and we find that in most data sets the scenario is statistical equivalent to Λ cold dark matter one. Moreover, we examine the big bang nucleosynthesis, and we show that the scenario satisfies the corresponding requirements. Additionally, we perform a phase-space analysis, and we show that the Universe past attractor is the matter-dominated epoch, while at late times the Universe results in the dark-energy-dominated solution. Finally, we show that Kaniadakis horizon-entropy cosmology accepts heteroclinic sequences, but it cannot exhibit bounce and turnaround solutions.

Recently, the research community has been exploring fractional calculus to address problems related to cosmology; in this approach, the gravitational action integral is altered, leading to a modified Friedmann equation, then the resulting theory is compared against observational data. In this context, dynamical systems can be used along with an analysis the phase spaces for different values of the fractional order of the derivative and their different matter contents. The equilibrium points are classified, providing a range for the order of the fractional derivative in order to investigate whether the cosmological history can be reconstructed and a late-time accelerating power-law solution obtained for the scale factor. In this paper, we discuss the physical interpretation of the corresponding cosmological solutions with particular emphasis on the influence of the fractional order of the derivative in a theory of gravity that includes a scalar field minimally coupled to gravity. The presented results improve and extend those obtained previously, further demonstrating that fractional calculus can play a relevant role in cosmology.

We study the dynamics of the field equations in a four-dimensional isotropic and homogeneous spatially flat Friedmann-Lemaître-Robertson-Walker geometry in the context of Einstein-Gauss-Bonnet theory with a matter source and a scalar field coupled to the Gauss-Bonnet scalar. In this theory, the Gauss-Bonnet term contributes to the field equations. The mass of the scalar field depends on the potential function and the Gauss-Bonnet term. For the scalar field potential, we consider the exponential function and the coupling function between the scalar field and the Gauss-Bonnet scalar is considered to be the linear function. Moreover, the scalar field can have a phantom behavior. We consider a set of dimensionless variables and write the field equations into a system or algebraic-differential equations. For the latter, we investigate the equilibrium points and their stability properties. We use compactified variables to perform a global analysis of the asymptotic dynamics. This gravitational theory can explain the Universe’s recent and past acceleration phases. Therefore, it can be used as a toy model for studying inflation or as a dark energy candidate.

Scalar field cosmologies with a generalized harmonic potential and a matter fluid with a barotropic Equation of State (EoS) with barotropic index $$\gamma $$γ for locally rotationally symmetric (LRS) Bianchi III metric and open Friedmann–Lemaître–Robertson–Walker (FLRW) metric are investigated. Methods from the theory of averaging of nonlinear dynamical systems are used to prove that time-dependent systems and their corresponding time-averaged versions have the same late-time dynamics. Therefore, simple time-averaged systems determine the future asymptotic behavior. Depending on values of barotropic index $$\gamma $$γ late-time attractors of physical interests for LRS Bianchi III metric are Bianchi III flat spacetime, matter dominated FLRW universe (mimicking de Sitter, quintessence or zero acceleration solutions) and matter-curvature scaling solution. For open FLRW metric late-time attractors are a matter dominated FLRW universe and Milne solution. With this approach, oscillations entering nonlinear system through Klein–Gordon (KG) equation can be controlled and smoothed out as the Hubble factor H – acting as a time-dependent perturbation parameter – tends monotonically to zero. Numerical simulations are presented as evidence of such behaviour.

Scalar field cosmologies with a generalized harmonic potential and a matter fluid with a barotropic equation of state (EoS) with barotropic index $$\gamma $$γ for the locally rotationally symmetric (LRS) Bianchi I and flat Friedmann–Lemaître–Robertson–Walker (FLRW) metrics are investigated. Methods from the theory of averaging of nonlinear dynamical systems are used to prove that time-dependent systems and their corresponding time-averaged versions have the same late-time dynamics. Therefore, the simplest time-averaged system determines the future asymptotic behavior. Depending on the values of $$\gamma $$γ , the late-time attractors of physical interests are flat quintessence dominated FLRW universe and Einstein-de Sitter solution. With this approach, the oscillations entering the system through the Klein–Gordon (KG) equation can be controlled and smoothed out as the Hubble parameter H – acting as time-dependent perturbation parameter – tends monotonically to zero. Numerical simulations are presented as evidence of such behavior.