    More than a century ago, Newtonian gravity failed to explain the anomalies in the orbits of Mercury. In 1915, Einstein introduced his theory of General Relativity (GR), which not only described the answer to the physics problems of that time, including the issue of Mercury but also its ability in its interesting predictions like the gravitational lensing, gravitational waves, black holes, etc.
However, despite all of these successes, in recent years this theory, failed to be consistent with observation and was not able to describe the dark energy and dark matter that are confirmed by observations. Also, GR was not renormalizable and, therefore, could not be quantized by using the conventional methods in quantum field theory. It was shown that if one demands the renormalization to be satisfied at one-loop, the Einstein-Hilbert action should be supplemented by higher-order curvature terms [1]}. Later on, Einstein's theory of GR had many forms of modifications such as: \(f(R)\)  gravity with \(R\)  as the Ricci scalar [2]}, [3]}, [4]}, \(f(R, T)\)  gravity, where \(T\)  is the trace of energy-momentum tensor [5]}, [6]}, \(f(G)\)  gravity with \(G\)  the Gauss-Bonnet invariant [7]} etc. All of these modified theories have received much attention to investigating the accelerated expansion of our universe, flat rotation curves of galaxies, wormhole behavior [8]}, [9]}, [10]}, [11]}. One of the familiar extensions of GR is to include terms such as
\(f(R)=...+\dfrac{\beta _{2}}{R^{2}}+\dfrac{\beta _{1}}{R}+R-\Lambda +\alpha _{2} R^2+\alpha _{3} R^3+...\) 
