Cosmological model from the holographic equipartition law with a modified R\'{e}nyi entropy

Cosmological equations were recently derived by Padmanabhan from the expansion of cosmic space due to the difference between the degrees of freedom on the surface and in the bulk in a region of space. In this study, a modified R\'{e}nyi entropy is applied to Padmanabhan's `holographic equipartition law', by regarding the Bekenstein--Hawking entropy as a nonextensive Tsallis entropy and using a logarithmic formula of the original R\'{e}nyi entropy. Consequently, the acceleration equation including an extra driving term (such as a time-varying cosmological term) can be derived in a homogeneous, isotropic, and spatially flat universe. When a specific condition is mathematically satisfied, the extra driving term is found to be constant-like as if it is a cosmological constant. Interestingly, the order of the constant-like term is naturally consistent with the order of the cosmological constant measured by observations, because the specific condition constrains the value of the constant-like term.


I. INTRODUCTION
To explain the accelerated expansion of the late universe, ΛCDM (lambda cold dark matter) models assume a cosmological constant Λ related to an additional energy component called dark energy [1][2][3]. However, Λ measured by observations is many orders of magnitude smaller than the theoretical value estimated by quantum field theory [4]. This discrepancy is the so-called cosmological constant problem. To resolve this theoretical difficulty, numerous cosmological models have been proposed [5,6], such as CCDM (creation of CDM) models [7][8][9] and Λ(t)CDM models, which assume a time-varying cosmological term Λ(t) [10][11][12][13][14][15][16][17][18][19]. In CCDM models, a constant term is obtained from a dissipation process based on gravitationally induced particle creation [7][8][9], while in Λ(t)CDM models, a constant term is obtained from an integral constant of the renormalization group equation for the vacuum energy density [12].
For these models, thermodynamic scenarios have attracted attention, in which the Bekenstein-Hawking entropy (which is proportional to the surface area of the event horizon) [20,21] and the holographic principle (which refers to the information of the bulk stored on the horizon) [22] play important roles. For example, using the holographic principle it has been proposed that gravity is itself an entropic force derived from changes in the Bekenstein-Hawking entropy [23][24][25][26]. Based on this concept, the cosmological equations have been extensively examined in a homogeneous and isotropic universe [27,28], although the cosmological constant has not been discussed. In an alternative treatment, Easson et al. [29] proposed an entropic cosmology that assumes the usually neglected surface terms on the horizon of the * E-mail: komatsu@se.kanazawa-u.ac.jp universe [30][31][32][33][34][35][36][37][38][39]. In entropic cosmology, an extra driving term to explain the accelerated expansion is derived from entropic forces on the horizon of the universe. An area entropy (the Bekenstein-Hawking entropy), a volume entropy (the Tsallis-Cirto entropy) [40], a quartic entropy [34,35], and a general form of entropy [38,39] have been applied to entropic cosmology. However, the entropic-force term (related to dark energy and Λ) is generally considered to be a tuning parameter, which makes it difficult to include in discussion of the cosmological constant problem.
Padmanabhan recently provided a new insight into the origin of spacetime dynamics using another thermodynamic scenario called the 'holographic equipartition law' [41]. Based on the holographic equipartition law with Bekenstein-Hawking entropy, cosmological equations in a flat universe can be derived from the expansion of cosmic space due to the difference between the degrees of freedom on the surface and in the bulk [41]. The emergence of the cosmological equations (i.e., the Friedmann and acceleration equations) has been examined from various viewpoints, such as a non-flat universe and quantumcorrected entropy [42][43][44][45][46]. However, dark energy and Λ have not yet been discussed fundamentally, though they have been considered in several studies [44][45][46]. This is likely because dark energy and Λ (which are related to an extra driving term) have been assumed to explain the accelerating universe. If the extra driving term is naturally derived from the holographic equipartition law, it is possible to study various cosmological models based on holographic equipartition law. However, such a driving term should not be derived using Bekenstein-Hawking entropy.
In this paper, the modified Rényi entropy is applied to the holographic equipartition law in a homogeneous, isotropic, and spatially flat universe. It is expected that an extra driving term can be derived from the holographic equipartition law with the modified Rényi entropy. Therefore, the present study should help in developing cosmological models based on the holographic equipartition law. In addition, the extra driving term is expected to be constant-like under specific conditions. The specific condition and constant-like term may provide new insights into the cosmological constant problem. (The present study focuses on the derivation of the extra driving term and the specific condition for the constantlike term. Accordingly, the inflation of the early universe and density perturbations related to structure formations are not examined.) The remainder of the article is organized as follows. In Sec. II, Λ(t)CDM models are briefly reviewed for a typical formulation of the cosmological equations. In Sec. III, entropies on the horizon are discussed. The Bekenstein-Hawking entropy is reviewed in Sec. III A, while a modified Rényi is introduced in Sec. III B. The holographic equipartition law is discussed in Sec. IV. In Sec. V, the modified Rényi entropy is applied to the holographic equipartition law, to derive the acceleration equation that includes an extra driving term. The extra driving term is then discussed in Sec. VI, focusing on the specific condition required for obtaining a constant-like term. Finally, in Sec. VII, the conclusions of the study are presented.
A homogeneous, isotropic, and spatially flat universe is considered, and the scale factor a(t) is examined at time t in the Friedmann-Lemaître-Robertson-Walker metric [35,38]. In the Λ(t)CDM model, the Friedmann equation is given as and the acceleration equation is where the Hubble parameter H(t) is defined by G, c, ρ(t), and p(t) are the gravitational constant, the speed of light, the mass density of cosmological fluids, and the pressure of cosmological fluids, respectively [35,38], and Λ(t) is a time-varying cosmological term. Based on Eqs. (1) and (2), the continuity equation [35] is given byρ .
The right-hand side of this continuity equation is usually non-zero, except for the case Λ(t)=Λ. Accordingly, the Λ(t)CDM model can be interpreted as a kind of energy exchange cosmology in which the transfer of energy between two fluids is assumed [58]. When Λ(t)=Λ, the Friedmann, acceleration, and continuity equations are identical to those for the standard ΛCDM model. In this paper, cosmological models based on the holographic equipartition law are assumed to be a particular case of Λ(t)CDM models, although the theoretical backgrounds are different.

III. ENTROPY ON THE HORIZON
In the holographic equipartition law, the horizon of the universe is assumed to have an associated entropy [41]. The Bekenstein-Hawking entropy is generally used, replacing the horizon of a black hole by the horizon of the universe. In Sec. III A, the Bekenstein-Hawking entropy is briefly reviewed. In Sec. III B, a novel type of Rényi entropy proposed by Biró and Czinner [50] is introduced for the entropy on the horizon of the universe.

A. The Bekenstein-Hawking entropy
The Bekenstein-Hawking entropy S BH [20] is given as where k B and are the Boltzmann constant and the reduced Planck constant, respectively [33,34]. The reduced Planck constant is defined by ≡ h/(2π), where h is the Planck constant. A H is the surface area of the sphere with the Hubble horizon (radius) r H , given by Substituting A H = 4πr 2 H into Eq. (5) and using Eq. (6), we obtain [33,34] or equivalently where K is a positive constant given by and L p is the Planck length, written as As shown in Eqs. (5) and (8), the Bekenstein-Hawking entropy S BH on the Hubble horizon is proportional to H −2 (and A H ) and is related to the Planck length.
In this paper, a spatially flat universe, k = 0, is considered, where k is a curvature constant. For a spatially non-flat universe (k = 0), the apparent horizon given by r A = c/ H 2 + (k/a 2 ) is used instead of the Hubble horizon, see, e.g., Ref. [43].

B. Modified Rényi entropy
Nonextensive entropies for black hole thermodynamics have been extensively investigated [40,50,51,54]. Recently, Biró and Czinner [50] suggested a novel type of Rényi entropy on black-hole horizons, by regarding the Bekenstein-Hawking entropy as a nonextensive Tsallis entropy and using a logarithmic formula. In the present study, the modified Rényi entropy is used for the entropy on the Hubble horizon. Therefore, the modified Rényi entropy [50,51] is briefly reviewed here.
In nonextensive thermodynamics, the Tsallis entropy S T for a set of W discrete states is defined as where p i is a probability distribution and q is the socalled nonextensive parameter [40,49]. In addition, the original Rényi entropy [40,47] is defined as When q = 1, both S T and S org R recover the Boltzmann-Gibbs entropy given by The original Rényi entropy [40,47] can be written as where, for simplicity, 1 − q has been replaced by λ: In this paper, λ is considered to be non-negative, as discussed later. A novel type of Rényi entropy has been proposed and examined in Refs. [50,51], in which not only is the logarithmic formula of the original Rényi entropy used but the Bekenstein-Hawking entropy S BH is considered to be a nonextensive Tsallis entropy S T . Using Eq. (14) and replacing S T by S BH , we obtain a modified Rényi entropy S R [51] given by When λ = 0 (i.e., q = 1), S R becomes S BH . We call S R the modified Rényi entropy. In the present study, the modified Rényi entropy is applied to the holographic equipartition law. To this end, the holographic equipartition law is reviewed in the next section. Note that the physical origin of the modified Rényi entropy is likely unclear at the present time. However, it is worthwhile to examine cosmological models based on the holographic equipartition law from various viewpoints. Therefore, in this paper, the modified Rényi entropy is applied to the holographic equipartition law.

IV. HOLOGRAPHIC EQUIPARTITION LAW
Padmanabhan [41] derived the Friedmann and acceleration equations in a flat universe from the expansion of cosmic space due to the difference between the degrees of freedom on the surface and in the bulk in a region of space. In this section, the Padmanabhan idea called the 'holographic equipartition law' is reviewed, based on his work [41] and other related works [42][43][44][45][46]. (The surface terms assumed in entropic cosmology [29] are not considered in the holographic equipartition law. In Ref. [46], the holographic equipartition law is examined, including the surface terms.) In an infinitesimal interval dt of cosmic time, the increase dV of the cosmic volume can be expressed as where N sur is the number of degrees of freedom on the spherical surface of Hubble radius r H and N bulk is the number of degrees of freedom in the bulk [41]. L p is the Planck length given by Eq. (10) and ǫ is a parameter discussed later. Equation (17) represents the so-called holographic equipartition law. Note that the right-hand side of Eq. (17) includes c, because c is not set to be 1 in the present paper. Using r H = c/H given by Eq. (6), the Hubble volume V can be written as The number of degrees of freedom in the bulk is assumed to obey the equipartition law of energy [41]: where the Komar energy |E| contained inside the Hubble volume V is given by and ǫ is a parameter defined as [41,42] ǫ ≡ +1 (ρc 2 + 3p < 0), From this definition, |E| is confirmed to be nonnegative. Note that ρc 2 + 3p < 0 corresponds to an accelerating universe (e.g., dark energy dominated universe), while ρc 2 + 3p > 0 corresponds to a decelerating universe (e.g., matter and radiation dominated universe). The temperature T on the horizon is written as This temperature is used for calculating N bulk from Eq. [59]. In contrast, the number of degrees of freedom on the spherical surface is given by where S H is the entropy on the Hubble horizon. When S H = S BH , Eq. (23) is equivalent to that in Ref. [41]. In the next section, S H is replaced by the modified Rényi entropy S R given by Eq. (16). We now derive the cosmological equations from the holographic equipartition law. To this end, we first calculate the left-hand side of Eq. (17). Substituting Eq. (18) into Eq. (17), the left-hand side of Eq. (17) is given by In this calculation, r has been set to be r H before the time derivative is calculated [41]. Next, we calculate N bulk included in the right-hand side of Eq. (17). According to Ref. [41], ρc 2 + 3p < 0 is selected and, therefore, ǫ = +1 from Eq. (21). Note that the following discussion is not affected even if ρc 2 + 3p > 0 is selected [42]. Substituting Eqs. (20) and (22) into Eq. (19), using Eq. (18), and rearranging, we obtain In addition, substituting ǫ = +1 and Eqs. (10), (23), (24), and (25) into Eq. (17) gives and solving this with regard toḢ, we havė where K is given by Eq. (9). As shown in Eq. (2),ä/a is written asä/a =Ḣ + H 2 . Substituting Eq. (27) into this equation, we obtain The above equation is the acceleration equation derived from the holographic equipartition law, where S H is the entropy on the Hubble horizon. Interestingly, the second term on the right-hand side, i.e., H 2 (1 − S H H 2 /K), appears to be an extra driving term. When the Bekenstein-Hawking entropy is used, i.e., S H = S BH , the second term H 2 (1 − S H H 2 /K) is zero because S BH = K/H 2 given by Eq. (7). Consequently, the acceleration equation can be written as The obtained cosmological equation is considered to be a particular case of Λ(t)CDM models, as mentioned in Sec. II. Accordingly, Eq. (29) indicates that Λ(t) in Eq.
In the above, the second term on the right-hand side of Eq. (28) is zero because the Bekenstein-Hawking entropy S BH is used for the entropy S H on the horizon. However, if a different S H is assumed, the second term (corresponding to an extra driving term) is expected to be non-zero. In the next section, we examine this expectation, applying the modified Rényi entropy S R to the holographic equipartition law. Note that a similar driving term has been studied, e.g., in Ref. [46], using quantum corrected entropy.

V. ACCELERATION EQUATION FROM THE HOLOGRAPHIC EQUIPARTITION LAW WITH A MODIFIED RÉNYI ENTROPY
In this section, the modified Rényi entropy S R is applied to the holographic equipartition law, instead of the Bekenstein-Hawking entropy. The modified Rényi entropy given by Eq. (16) can be written as where S R becomes S BH when λ = 0. The acceleration equation derived from the holographic equipartition law given by Eq. (28) is where S H is the entropy on the Hubble horizon and K given by Eq. (9) is written as Regarding S R as S H and substituting Eq. (32) into Eq. (33), we havë where, using S BH = K/H 2 given by Eq. (7), the extra driving term f (H) is written as and λ is a non-negative constant, i.e., λ ≥ 0. The acceleration equation given by Eq. (35) (39) is plotted. The two approximation curves for 1/2 and 1/X given by Eq. (42) are also shown.
The positive dimensionless parameter X is defined as where λ is considered to be positive, i.e., λ > 0. Dividing Eq. (37) by λK, and applying Eq. (38), we have the normalized extra-driving term given by To observe the properties of f (H), the normalized term is plotted in Fig. 1. Consequently, f (H)/(λK) gradually approaches 1/2 with decreasing X, while it approaches a common curve with increasing X. Therefore, we focus on two specific conditions, X ≤ 1 and X ≫ 1. Under the two conditions, the approximate formulas are given mathematically by High-order terms have been neglected in Eq. (40). Substituting Eqs. (40) and (41), respectively, into Eq. (39), we have From Eqs. (38) and (42), f (H) can be written as The constant-like and H 2 -like terms are respectively obtained from each condition, as shown in Eq. (43). The constant-like term, λK/2, can be interpreted as a kind of cosmological constant, although the specific condition λK/H 2 ≤ 1 is required. In the next section, we discuss the constant-like term and the specific condition. As a matter of fact, the extra driving term (which behaves like effective dark energy) can be derived even if dark energy is not assumed, i.e., even when p = 0. Alternatively, the modified Rényi entropy is assumed in the present paper. In this sense, the original holographic equipartition law is modified. This modification is expected to provide new insights into cosmological models based on the holographic equipartition law.
As shown in Eq. (43), H 2 -like terms are derived under a specific condition. However, the driving term is not likely to be related to the inflation of the early universe, because higher exponents terms such as H 4 terms should be required for the inflation [38]. In Λ(t)CDM models, the higher exponents terms have been closely examined in Ref. [13]. In entropic cosmology, quantum corrected entropy is introduced for the higher exponents terms, see, e.g., Ref. [46] and the second paper of Ref. [29].

VI. CONSTANT-LIKE TERM UNDER A SPECIFIC CONDITION
The acceleration equation can be derived from the holographic equipartition law with the modified Rényi entropy, as examined in the previous section. Interestingly, an extra driving term included in the equation is found to be constant-like under a specific condition. In this section, we discuss the constant-like term under the specific condition. To this end, the present cosmological model is reviewed, assuming that it is a particular case of Λ(t)CDM models. Substituting Λ(t)/3 = f (H) into Eqs.
(1), (2), and (4), we obtain the Friedmann, acceleration, and continuity equations written as where the extra driving term f (H) is given by and K is As shown in Eq. (46), the right-hand side of the continuity equation is not zero generally. This non-zero righthand side may be interpreted as the interchange of energy between the bulk (the universe) and the surface (the horizon of the universe) [31,38]. (The previous works of Solà et al. [14,15] imply that the value of the nonzero right-hand side of the continuity equation should be small in Λ(t)CDM models [38].) In the following, the right-hand side of the continuity equation is zero, because f (H) is considered to be constant under a specific condition. General solutions for the present model, which includes a logarithmic term, are separately discussed in Appendix A. This paper focuses on the constant-like term under a specific condition. From Eq. (43), the constant-like term and the specific condition can be written as where H 0 is the Hubble parameter at the present time and λK/(2H 2 0 ) is a dimensionless constant (corresponding to α 4 in Ref. [35]). Hereafter, we call the constantlike term the constant term. Under the specific condition given by Eq. (49), the Friedmann, acceleration, and continuity equations are written as ρ + 3ȧ a ρ + p c 2 = 0.
The three equations are essentially the same as those in the standard ΛCDM model. Setting p = 0 for nonrelativistic matter, the background evolution of the universe is analytically given by where a 0 represents the scale factor at the present time and Ω h is a holographic parameter defined by Using Eq. (54) andH ≡ H/H 0 , Eq. (38) can be written as In the present model, the constant term λK/2 and the holographic parameter Ω h have been restricted by the inequality given in Eq. (49). Accordingly, the constraint on both λK/2 and Ω h can be discussed, without tuning. In the history of an expanding universe, H 0 is expected to be the minimum value of H, because r H increases with time and H is given by H = c/r H from Eq. (6). Therefore, when H = H 0 , the most severe constraint on λK can be obtained from the inequality of Eq. (49). The constraint is written as In the present model, Eq. (56) is mathematically required to obtain the constant term. From Eq. (56), the order of the constant term, λK/2, is approximately given by or equivalently, from Eqs. (54) and (56), the order of Ω h can be approximately written as We now compare the constant term λK/2 in the present model and the cosmological constant term Λ/3 in the standard ΛCDM model. In the ΛCDM model, the density parameter Ω Λ for the cosmological constant Λ is given by Numerous cosmological observations indicate that the order of Ω Λ is expected to be 1, e.g., Ω Λ = 0.692 from the Planck 2015 results [3]. Accordingly, the order of Ω Λ can be approximately written as Using Eqs. (59) and (60), we have Equations (57) and (61) imply that the order of the constant term λK/2 is consistent with the order of Λ/3. That is, interestingly, the constant term in the present model is naturally consistent with the order of Λ measured by cosmological observations as if it is Λ. Similarly, from Eqs. (58) and (60), the order of Ω h is likely consistent with the order of Ω Λ . Of course, in the present model, λK H 2 0 is required to obtain the constant term, assuming the modified Rényi entropy instead of the Bekenstein-Hawking entropy. This extension may be beyond the original holographic equipartition law. However, the constant term and the specific condition considered here may provide new insights into the cosmological constant problem. Equation (56) can be written as λ L 2 p /(πk B r 2 H0 ) using Eq. (48) and H 0 = c/r H0 , where r H0 is the Hubble horizon at the present time. This constraint indicates that λ is an extremely small positive value. In other words, the modified Rényi entropy is approximately equivalent to the Bekenstein-Hawking entropy although it slightly deviates from the Bekenstein-Hawking entropy. Accordingly, the present model may imply that the cosmological constant is related to a small deviation from the Bekenstein-Hawking entropy. As time passes, the constraint is expected to be more severe because r H0 increases with time. Consequently, the extra driving term in the present model should gradually deviate from the constant value, even if the severe condition is satisfied at the present time. That is, the present model can be distinguished from the ΛCDM model, as discussed in Appendix A.
In this paper, we have focused on the derivation of an extra driving term from the holographic equipartition law with the modified Rényi entropy. The obtained term f (H) given by Eq. (47) is a logarithmic formula. However, under two specific conditions, it can be systematically written as f (H) = C 0 H 2 0 + C 1 H 2 , where C 0 and C 1 are dimensionless constants. In the Λ(t)CDM model, the above formula has been examined in detail [15][16][17]32]. For example, Gómez-Valent et al. [16] and Basilakos et al. [32] have shown that it is not the H 2 ,Ḣ, and H terms but rather the constant term that plays an important role in cosmological fluctuations related to structure formations [38]. In addition, recently, Λ(t)CDM models, which include power series of H, have been found to be more suitable than the standard ΛCDM model [14]. In particular, the simple combination of the constant and H 2 terms, f (H) = C 0 H 2 0 + C 1 H 2 , is likely favored; see e.g., the works of Solà et al. [15], Gómez-Valent et al. [16], and Lima et al. [17], which indicate that C 0 is dominant and C 1 is expected to be small [14,15]. The smallness of C 1 has not yet been explained by the holographic approach [38], though it may be explained by a deeper understanding of the present model. For example, the logarithmic formula discussed in the present model slightly deviates from the constant value when λK/H 2 ≈ 1, as shown in Fig. 1. Of course, this small deviation is not proportional to H 2 . However, the small deviation may be related to the smallness of C 1 , if the smallness can be interpreted as a deviation from a constant term. This task is left for future research.
Keep in mind that the present result depends on the choice of entropy. In this paper, the modified Rényi entropy [50,51] is employed, in which the Bekenstein-Hawking entropy is considered to be a nonextensive Tsallis entropy. That is, the Bekenstein-Hawking entropy is assumed to satisfy the Tsallis composition rule. In addition, a logarithmic formula of the original Rényi entropy is used for the modified Rényi entropy [50,51], where the logarithmic formula is obtained from the Tsallis composition rule. For details of the modified Rényi entropy, see Refs. [50,51]. The physical origin of the modified Rényi entropy is likely unclear at the present time. However, it is worthwhile to examine cosmological models based on the holographic equipartition law from various viewpoints. Note that the choice of temperature does not affect main results in the present study, because the extra driving term discussed here is related to N sur = 4S H /k B from Eq. (23) and is independent of the temperature [59].

VII. CONCLUSIONS
Recently, a novel type of Rényi entropy on a black-hole horizon was proposed by Biró and Czinner [50], in which a logarithmic formula of the original Rényi entropy was used and the Bekenstein-Hawking entropy was regarded as a nonextensive Tsallis entropy. In this study, the modified Rényi entropy has been applied to the holographic equipartition law proposed by Padmanabhan [41], to investigate cosmological models based on the holographic equipartition law. Consequently, the acceleration equation, which includes an extra driving term, can be derived in a homogeneous, isotropic, and spatially flat universe. The extra driving term is a logarithmic formula which can explain the accelerated expansion of the late universe, as for effective dark energy.
Under two specific conditions, the extra driving term is found to be the constant-like and H 2 -like terms, respectively. In particular, when a specific condition is satisfied, the extra driving term is found to be constantlike, i.e., λK/2, as if it is a cosmological constant Λ. In the present model, λK H 2 0 is required to obtain the constant-like term from the most severe constraint. In other words, the constant-like term must be extremely small when the extra driving term behaves as if it is constant. The present model may imply that the cosmological constant is related to a small deviation from the Bekenstein-Hawking entropy. In addition, interestingly, the order of the constant-like term is naturally consistent with the order of the cosmological constant measured by observations, because the specific condition constrains the value of the constant-like term. The present model may provide new insights into the cosmological constant problem.
This paper focuses on the derivation of the extra driving term and the specific condition for the constantlike term. Accordingly, general solutions for the present model that includes a logarithmic driving term have been separately studied in Appendix A. To solve the cosmological equations, the present model is assumed to be a particular case of Λ(t)CDM models. From the obtained solution, the background evolutions of the universe in the present model are found to agree well with those in the standard ΛCDM model when both 2Ω h /H 2 and 2Ω h are small. For lower redshift, the present model gradually deviates from the ΛCDM model, due to the logarithmic term. Therefore, the present model can be distinguished from the standard ΛCDM model.

,
(A7) and the normalized Hubble parameterH is defined as Similarly, the normalized scale factorã is defined as Substituting H =HH 0 and a =ãa 0 into Eq. (A5), and arranging the resultant equation, we have In addition, a parameter N is defined by N ≡ lnã, and therefore dN = dã a .
Using Eq. (A11), Eq. (A10) can be written as When Ω h is constant, Eq. (A12) can be integrated as This solution is given by where C is an integral constant and 'li' represents the logarithmic integral [60] defined as From Eqs. (A8) and (A9), the present values ofH andã are 1. SubstitutingH = 1 andã = 1 into Eq. (A14), the integral constant C can be written as

Background evolutions of the late universe
In this subsection, the background evolutions of the late universe in the present model are briefly examined. For this purpose, Ω h is considered to be a free parameter. Note that, when Eq. (A17) for Ω h = 0 is calculated numerically, Ω h is set to be 10 −4 .
To observe the properties of the present model, the luminosity distance is examined. The luminosity distance [61] is generally given by where the integrating variable y, the function F (y), and the redshift z are given by The relationship betweenã andH for each Ω h is obtained from Eq. (A17). Using this relationship, the luminosity distance is calculated from Eqs. (A18) and (A19). For typical results, the luminosity distances for Ω h = 0, 0.4, and 0.8 are plotted in Fig. 2. The luminosity distance for the present model tends to be more consistent with supernova data points with increasing Ω h . This result indicates that the present model can describe an accelerated expansion of the late universe. However, the influence of the increase in Ω h on d L is likely to be smaller than the influence of the increase in Ω Λ in the ΛCDM model on d L (the result for the ΛCDM model is not shown). To compare the two models, a temporal deceleration parameter is examined.
The temporal deceleration parameter q is defined by The temporal deceleration parameter q for the present model can be calculated from Eq. (A22). Note that the relationship betweenã andH is obtained from Eq. (A17). For the ΛCDM model, substituting f (H) = Λ/3 into Eq. (A21), and using Ω Λ = Λ/(3H 2 0 ) andH = H/H 0 , we have whereH 2 is given byH 2 = (1 − Ω Λ )ã −3 + Ω Λ . In a flat universe, the density parameter Ω m for matter is given  by Ω m = 1 − Ω Λ , neglecting the density parameter for the radiation. Figure 3 shows the dependence of the temporal deceleration parameter q on the redshift z. As expected, q for Ω h = 0 and 0.2 is consistent with that for Ω Λ = 0 and 0.2, respectively. However, q for the present model gradually deviates from that for the ΛCDM model, with increasing Ω h , especially for low redshift. In fact, even for Ω h = 0.2, q gradually deviates from that for Ω Λ = 0.2 with decreasing z, although the two results agree well for high redshift. As shown in Fig. 4, 2Ω h /H 2 increases with decreasing z, becauseH decreases with decreasing z. That is, the deviation from the ΛCDM model for low redshift (shown in Fig. 3) is related to the increase in 2Ω h /H 2 . In addition, for Ω h = 0.6 and 0.8, q varies from positive to negative with decreasing z, as shown in Fig. 3. Therefore, the present model can describe a decelerated and accelerated expansion of the late universe.
In this Appendix, the general solution for the present model has been examined, assuming that the present model is a particular case of Λ(t)CDM models. In addition, background evolutions of the late universe in the present model have been discussed. Consequently, the present model is found to agree well with the standard ΛCDM model, when both 2Ω h /H 2 and 2Ω h are small. For lower redshift, the present model gradually deviates from the ΛCDM model, due to the properties of the logarithmic driving term. Accordingly, the present model examined here can be distinguished from the ΛCDM model.