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.

Our paper introduces a new theoretical framework called the Fractional Einstein–Gauss–Bonnet scalar field cosmology, which has important physical implications. Using fractional calculus to modify the gravitational action integral, we derived a modified Friedmann equation and a modified Klein–Gordon equation. Our research reveals non-trivial solutions associated with exponential potential, exponential couplings to the Gauss–Bonnet term, and a logarithmic scalar field, which are dependent on two cosmological parameters, m and ${\alpha }_0=t_0{H}_0$ and the fractional derivative order $\mu$. By employing linear stability theory, we reveal the phase space structure and analyze the dynamic effects of the Gauss–Bonnet couplings. The scaling behavior at some equilibrium points reveals that the geometric corrections in the coupling to the Gauss–Bonnet scalar can mimic the behavior of the dark sector in modified gravity. Using data from cosmic chronometers, type Ia supernovae, supermassive Black Hole Shadows, and strong gravitational lensing, we estimated the values of m and ${\alpha }_0$, indicating that the solution is consistent with an accelerated expansion at late times with the values ${\alpha }_0=1.38\pm 0.05$, $m=1.44\pm 0.05$, and $\mu =1.48\pm 0.17$ (consistent with ${\mathrm{\Omega }}_{m,0}=0.311\pm 0.016$ and $h=0.712\pm 0.007$), resulting in an age of the Universe $t_0=19.0\pm 0.7$ [Gyr] at 1$\sigma$ CL. Ultimately, we obtained late-time accelerating power-law solutions supported by the most recent cosmological data, and we proposed an alternative explanation for the origin of cosmic acceleration other than $\mathrm{\Lambda }$CDM. Our results generalize and significantly improve previous achievements in the literature, highlighting the practical implications of fractional calculus in cosmology.

Fractional differential equations model processes with memory effects, providing a realistic perspective on complex systems. We examine time-delayed differential equations, discussing first-order and fractional Caputo time-delayed differential equations. We derive their characteristic equations and solve them using the Laplace transform. We derive a modified evolution equation for the Hubble parameter incorporating a viscosity term modeled as a function of the delayed Hubble parameter within Eckart’s theory. We extend this equation using the last-step method of fractional calculus, resulting in Caputo’s time-delayed fractional differential equation. This equation accounts for the finite response times of cosmic fluids, resulting in a comprehensive model of the Universe’s behavior. We then solve this equation analytically. Due to the complexity of the analytical solution, we also provide a numerical representation. Our solution reaches the de Sitter equilibrium point. Additionally, we present some generalizations.

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.