1. Introduction
Einstein’s theory of gravity predicts the accelerating expansion of the late-time universe, which has
been confirmed by the cosmological observations from type Ia Supernovae (SNIa) [1, 2]. The expansion
of the universe can be explained by the simplest model of the cosmological constant that currently
gives the best fit with the current observational data from Cosmic Microwave Background (CMB)
radiation, the measurement of Type Ia Supernovae (SNIa), Baryon Acoustic Oscillation (BAO), and
H0measurement [3, 4]. However, some issues remain for the cosmological constant or ΛCDM model
as encountered from the observations [3, 5, 6] and a seriously theoretical inconsistency [7, 8], and for
a review, one may refer to [9, 10]. In these regards, several modified theories of gravity have been
proposed to refine or provide alternatives. Without exhaustion, we list down some of the theories which
deviate from the standard cosmological model: scalar-tensor theory [11, 12, 13, 14, 15], vector-tensor
[16], tensor-tensor theory [17, 18, 19], massive gravity [20, 21, 22], higher-order derivative theory, in
particular, the f(R)-gravity for cosmology [23, 24, 25, 26], and higher-dimensional and string-motivated
theory [27, 28, 29, 30, 42, 43]. For a nice review, one is referred to [31]. Also, there are many dark energy
models under these frameworks such as Chevallier Polarski Linder (CPL) model [32], Holographic Dark
Energy model and its constraint [33, 34, 35, 36, 37], Generalized holographic cosmology via AdS/CFT
constraint with low redshift [38], and Chaplygin gas model [39, 40] to name a few. One can refer
to [41] for model comparison using observational data from Planck 2015.
It is well-known that Lovelock’s theorem is a natural generalization to general relativity, which
satisfies the diffeomorphism and local Lorentz invariance in D-dimensional space-time and, in partic-
ular, leads to second-order field equations [42, 44, 45]. A special case of this is the second order in
Lanczos-Lovelock gravity Lagrangian density so-called Einstein-Gauss-Bonnet gravity. Consequently,
one possible modification of General Relativity is the D≥5 Gauss-Bonnet (GB) gravity which sat-
isfies the properties in Lovelock’s theorem, including the ghost-free [46], natural generalization with
Einstein, and cosmological terms [47]. For its study of higher dimensions in cosmology, one can refer
to [48, 49, 50]. However, the GB term in 4Dis topologically invariant, so it does not modify Ein-
stein’s theory since it does not contribute to gravitational dynamics. Nevertheless, the GB term can
contribute to the 4Ddynamical equations if one introduces a coupling function of a scalar field to the
GB term [51, 52]. Due to the extra scalar degree of freedom, this consideration leads to many gravity
studies, including cosmic acceleration in inflationary models of the early universe, and tested against
the observational data [53, 54, 55, 56, 57, 58, 59].
Recently, the 4DEGB model proposed in [60] introduces a scaling coupling constant α→α/(D−
4), where Dis the space-time dimensions, and considers the limit D→4. This model gives rise to
nontrivial contributions to the gravitational dynamics due to the extra contribution in the equations
of motion (EoM) from the GB term as the consequence of the scaling coupling constant without
introducing extra degrees of freedom. Notice that in this model the space-time was assumed to be
continuous [61]. To be less exhaustive, there is an interesting work in cosmology adopting this model
with the observational constraints which can resolve the coincidence problem [62] and also indicate
that the re-scaling coupling constant of the model still needs the help of the cosmological constant to
explain the accelerated expansion. For a review of the model, one may refer to [63].
Due to its richness and peculiar consequences, comments and arguments on this novel 4DEGB
theory have been raised [64, 65, 66, 67, 68, 69, 70], and the diffeomorphism invariant regularization
is also considered in [71, 72, 73, 74] to point out that the theory advertised in [60] is the subclass
of Horndeski theory. However, the arguments from Horndeski theory still rely on a 2 + 1 degrees
of freedom where higher dimensions are compactified and where in [60] the space-time is based on
the D-dimensional direct product and D→4 limit, and further investigation into this regularized
tensor-scalar theory leading to a strong-coupling scalar field is studied in [73, 75] and without the
strong coupling field in [76]. Another investigation into the regularized 4DEGB theory is that
looking at diffeomorphism invariant property reveals the inconsistency in the regularized 4DEGB
theory [67, 69, 70]. Nevertheless, the consistent study of D→4 EGB theory can be achieved up
to spatially diffeomorphism invariance by using Hamiltonian formalism so-called minimally modified
gravity theory; otherwise, an extra degree of freedom is required [77, 78]. Therefore, we will proceed
by introducing an additional degree of freedom but still exploiting the regularization scheme which
1