Using entropy bounds to avoid the cosmological singularity and constrain cosmological particle production

In this work, we study the applications of entropy bounds in two toy cosmological models with particle production (annihilation), i.e., a radiation-dominated universe and a dust-dominated universe. We consider the co-moving volume and the volume covered by the particle horizon of a given observer as the thermodynamic systems satisfying entropy bounds. For the Bekenstein bound and the spherical entropy bound, it is found that the cosmological singularity can be avoided and cosmological particle production needs to be truncated in some special cases. Our study can be extended to other cosmological models with particle production.


Introduction
Black hole thermodynamics [1][2][3][4], which leads to the formulation of the holographic principle [5][6][7], may be the key to the understanding of quantum gravity.In 1972, Bekenstein found that black holes could have entropy [1], which triggered the research enthusiasm on black hole thermodynamics.In these studies, entropy is a crucial factor in the link between black holes and thermodynamics.However, for nearly a decade after Bekenstein found that black holes have entropy, few researchers considered the relationship between the entropy of black holes and the entropy of other gravitational systems.Until 1981, when Bekenstein studied the generalized second law (of thermodynamics) for a black hole [8], he for the first time argued that the generalized second law implies the entropy of any weakly gravitating matter system in asymptotically flat space should satisfy a bound S ≤ 2πkER/(ℏc), where E is the total mass-energy of the system and R is the radius of the smallest sphere that fits around the system [9].This entropy bound is also called the Bekenstein bound, which is independent of the gravitational theory.It is worth mentioning that Unruh and Wald did not agree with the original derivation of the Bekenstein bound [10,11].They stated that the entropy bound of a black hole is not needed for the validity of the generalized second law if there exists buoyancy force of the thermal atmosphere near the black hole horizon.Subsequent studies show that, according to physical processes, gravitational theories, the background space-time, etc., entropy bounds may have different forms.In 1995, Susskind argued that applying the generalized second law to a transformation that a system is converted to a black hole, one can get the spherical entropy bound S ≤ kA/(4l 2 p ), where A is a properly defined area enclosing the system [6,12].Subsequently, inspired by the work of Fischler and Susskind [13], Bousso proposed a covariant entropy bound (which is also called Bousso bound) S[L(B)] ≤ A(B)/4, where A is the area of the boundary B [14].The covariant entropy bound can be applied to any space-time including the strong gravitational system and satisfies general covariance, but it is only applicable to general relativity.In Ref. [15], the authors provided two ways to prove the covariant entropy bound and put forward a stronger entropy bound.In the year that the covariant entropy bound was proposed, Brustein and Veneziano proposed a causal entropy bound [16].The following year, Verlinde proposed the Bekenstein-Verlinde bound, the Bekenstein-Hawking bound and the Hubble bound [17].For more research on entropy bounds, one can refer to Refs.[18][19][20][21][22][23][24][25].
Although the concept of entropy bounds is the product of black hole research, it also has significant implications in cosmology.Some models of the Big Bang theory predict that there exists a cosmological singularity (initial singularity or Big Bang singularity) before the Big Bang, which contained all the energy and space-time of the universe.The Big Bang theory fits in well with cosmological observations and has been accepted by many physicists, but the initial singularity of the universe has been criticized.As a result, some cosmological models and theories have been proposed to explain or avoid the cosmological singularity, such as the cyclic model of the universe, multiverse, loop quantum gravity, etc.In 1989, Bekenstein found that the cosmological singularity is thermodynamically irrational [26] from the perspective of the entropy bound proposed by himself [9].Recently, Powell et al. proposed a re-examination of Bekenstein's approach in a radiation-dominated universe and also verified that the Bekenstein bound can be a feasible alternative to avoiding the cosmological singularity [27].In addition to the application of entropy bounds to the initial singularity of the universe, entropy bounds may also help us deduce the shape of the universe.For a closed universe, one can find that there exists a contradiction between the Fischler-Susskind bound and the positive curvature, which means that the Fischler-Susskind bound requires the shape of the universe to be non-closed [13].For more cosmological applications of various entropy bounds, one can refer to Refs.[25,[28][29][30][31][32][33].
It is an indisputable fact that the production and annihilation of particles occur continuously in our universe, so the cosmological model with particle production (annihilation) is more consistent with the real universe.The mechanisms and applications of particle production in cosmology have long attracted the attention of physicists (see a recent review [34] and references therein).In the 1960s, Parker established the first micromechanism of particle production in the context of cosmology with the quantum field theory in curved space-time [35][36][37][38].Then, particle production in different background space-time has been investigated successively [39][40][41][42].A few years ago, Harko found that, from a thermodynamic point of view, gravitational induced particle production can happen in non-minimal coupling theories [43,44], which provides a new mechanism for cosmological particle production.At present, most of the studies suggest that the production and annihilation of particles in an expanding universe have a profound theoretical foundation and they can not be ignored in the universe, especially in the early period.As for the applications of particle production in cosmology, they include multiple aspects, such as avoiding the cosmological singularity [45][46][47], explaining entropy production of the universe [45][46][47][48][49], accelerating the expansion of the universe [47,[50][51][52][53][54], and triggering inflation [55,56].
In this work, we study entropy bounds in the cosmological model with particle production (annihilation).We use entropy bounds to judge whether the cosmological singularity is thermodynamically rational in a radiation-dominated universe and a dust-dominated universe.Moreover, entropy bounds may be able to constrain cosmological particle production, which is a new application of entropy bounds in cosmology.The constraint on cosmological particle production is of significance for cosmology because it could provide us with truncations for some interactions in the context of cosmology.In the past, we usually used the (generalized) second law of thermodynamics to constrain cosmological particle production and the results were not unsatisfactory because it can only constrain the sign of the particle production rate [49,57].
The paper is organized as follows.Sec. 2 is a brief review on cosmological particle production and the corresponding entropy.In Sec. 3, we discuss the Bekenstein bound and the spherical entropy bound in a radiation-dominated universe with particle production.We focus on the cosmological singularity and the constraint on cosmological particle production inside the co-moving volume and the volume covered by the particle horizon of a given observer.Then, in Sec. 4, we study similar content in a dust-dominated universe with particle production.The last part, Sec. 5, is conclusions and discussions.

Cosmological particle production and entropy
In this section, we discuss cosmological particle production and entropy in the context of a homogeneous and isotropic universe, which can be described by the Friedmann-Lematîre-Robertson-Walker metric: For simplicity, we consider that the universe is spatially flat, i.e., k = 0. Here, a(t) is the scale factor.In this work, all the components of the universe are regarded as ideal fluids, so the energy-momentum tensor of all the ideal fluids is given as where ρ and p represent the total energy density and pressure of the ideal fluids, respectively.The four-velocity of the ideal fluids satisfies u ν u ν = −1.For general relativity, the Friedmann equations are given by Now, we focus on a spherical system with radius R in the universe.According to the purpose of research, the physical meaning of R can be multiple, such as the scale factor, the particle horizon, the apparent horizon, the radius of the visible universe, etc.Assuming that there are N particles in the system, one can define the particle production rate as Usually, if R is the scale factor and there is no interaction between these particles and other matter (or the background space-time), we have Γ = 0.For other cases of R, even if there is no interaction, Γ is generally nonzero due to the evolution of R. If the particle number in the system is non-conserved, the entropy of the system will be affected by the production (increase) or annihilation (decrease) of particles.We label the current particle number as N 0 and set the current radius of the system to R 0 .If the entropy of these particles is extensive ‡ in a homogeneous and isotropic gravitational system (such as the cosmological model we study), then the entropy of the system at time t can be written as ‡ In this work, we do not consider the non-extensive statistical entropy, such as Tsallis entropy [58].Some studies indicate that, compared to classical statistics, non-extensive statistics may be more applicable to a gravitational system [58][59][60].
where n t = N t /V t = N t /R 3 t is the particle number density and s(n t ) is the entropy per volume at time t.The radius R t of the system at time t is dependent on the way the universe is expanding.The particle number N t at time t is given by If the system is isolated, one can constrain the evolution of the system by the second law of thermodynamics: On the other hand, if the system possesses area entropy (such as the system covered by the apparent horizon [7,[61][62][63][64][65][66][67][68][69]), then the system can be constrained by the generalized second law of thermodynamics [1,2,8]: The area entropy S A of the system is usually proportional to the area of the system, but the coefficient is related to the gravitational theory.Owing to inequality (8) [or inequality (9)], the particle production rate Γ in the universe can not be arbitrary.Therefore, the (generalized) second law of thermodynamics is a means of constraining cosmological particle production (or the interaction between different substances in the context of cosmology).However, for a general radius R, the system is not isolated, so it may be not appropriate to use the (generalized) second law of thermodynamics to constrain cosmological particle production.In this work, we omit the area entropy of the system and try to constrain cosmological particle production with entropy bounds.Since entropy bounds do not require the corresponding system to be closed, we can study cosmological particle production in any system.On the other hand, we will also examine the effect of cosmological particle production on the cosmological singularity in the view of entropy bounds § .The entropy bounds employed in this work are the Bekenstein bound [9] and the spherical entropy bound [6,7,12].In order to obtain analytical solutions of the Friedmann equations, each of the toy cosmological models we consider only contains a species of matter.

Entropy bounds and particle production in a radiation-dominated universe
In 1989, Bekenstein studied the particle horizon of a given observer in a radiation-dominated universe and he found the cosmological singularity is thermodynamically impossible by considering the Bekenstein bound [26].In Bekenstein's Friedmann model, there is no cosmological particle production inside the co-moving volume, which may be inappropriate for the early universe because non-relativistic particles were not decoupled with other matter in the early universe.In this section, we extend Bekenstein's Friedmann model by introducing particle production into the radiation-dominated universe and discuss the relationship among the cosmological singularity, particle production and entropy bounds.As a toy model, we simplify the radiation-dominated universe as a universe containing only photons.The production of photons can be ascribed to the coupling between photons and the background spacetime [43,44] (or the running vacuum [49,70,71]).Since we do not know how to define the entropy of § In this work, we are not concerned about the mechanism by which matter was produced (annihilated) at the beginning of the universe.
the background space-time, the total entropy of the universe is considered to be only dependent on photons.We assume that the number of photons per unit volume is where c is the light speed, k is the Boltzmann constant, ζ(n) is the Riemann zeta function.If there is no coupling between photons and the background space-time (i.e., there is no particle production in the co-moving volume), the temperature T of photons should be proportional to a −1 .For a spherical system with radius R in the universe, with Eq. ( 5) and N = n R 3 , one can obtain the production rate of photons in the system is given as Since the entropy of photons per unit volume is given by the entropy of photons inside the spherical system with radius R is where N 0 is the current number of photons in the system.When R and Γ are given, we can calculate the entropy of photons in the system at any time.Next, we consider that R is the scale factor (i.e., the system is the co-moving volume) or the particle horizon of a given observer (i.e., the system is the volume covered by the particle horizon).As t approaches to zero, the production (annihilation) of photons may influence the cosmological singularity from the perspective of entropy bounds.Moreover, as t increases, entropy bounds may require that there exists a truncation for the production of photons.

Co-moving volume
If R is the scale factor, according to Eqs. (11) and (13), the entropy of photons inside the co-moving volume is given as If Γ = 0, then S is a constant.We first consider the Bekenstein bound, which requires the entropy of any weakly gravitating matter system to satisfy [9] where E is the total mass-energy including any rest mass and R is the radius of a sphere that can enclose the given system.For the photons inside the co-moving volume, we have E = π 2 k 4 15c 3 ℏ 3 a 3 T 4 , so the Bekenstein bound can be expressed as When Γ = 0, the energy density of photons inside the co-moving volume is given by ρ = ρ 0 a 4 (the Friedmann equations) or ρ = π 2 k 4 15c 3 ℏ 3 T 4 (thermodynamic state functions for black-body photons), so (16), then the Bekenstein bound can be simplified as where k ∼ 1.38 * 10 −23 J/K and ℏc ∼ 3.16 * 10 −26 J•m.If we set a 0 = 1 (m) and T 0 ∼ 2.7 K, then this inequality is tenable.Therefore, for a radiation-dominated universe without particle production, if its current temperature is the same as our real universe, applying the Bekenstein bound to the co-moving volume can not avoid the cosmological singularity.Moreover, based on inequality (17), one can find that if T 0 < 2 * 10 −4 K, the entropy of photons inside the co-moving volume will conflict with the Bekenstein bound.
When Γ ̸ = 0, a T will evolve over time.We can rewrite Eq. ( 11) as Then, one can obtain Substituting it into inequality ( 16) to cancel aT and N 0 yields where we can set a 0 = 1 (m) and T 0 ∼ 2.7 K. For 0 ≤ t < t 0 and Γ > 0, there must be a critical time t c corresponding to the equality sign of inequality (20).When t < t c , inequality (20) will be violated, which means that the Bekenstein bound requires the initial time of the universe to be nonzero.However, due to the existence of Γ, the scale factor does not satisfy a(t = 0) = 0, so the nonzero initial time is not equivalent to the nonzero initial volume.In the later discussion, we will see that, for Γ > 0, the universe naturally has no initial singularity according to the solution the scale factor.Therefore, in this case we do not need the Bekenstein bound to avoid the cosmological singularity.We will focus on the effect of the Bekenstein bound on the initial volume of the universe.Regarding this, we can see detailed examples later.
For 0 ≤ t < t 0 and Γ ≤ 0, since the right-hand side of inequality (20) increases with the decrease of t, the Bekenstein bound does not help to avoid the cosmological singularity.
For t 0 ≤ t and Γ ≥ 0, the right-hand side of inequality (20) increases with t.Based on the previous analysis, the Bekenstein bound is always true.Therefore, the photons inside the co-moving volume can be continuously produced and so the particle production rate Γ can not be limited by the Bekenstein bound.
For t 0 ≤ t and Γ < 0, as t increases, the entropy of photons inside the co-moving volume will violate the Bekenstein bound.Therefore, photons in the universe can not keep annihilating and the interaction between photons and the background space-time must be truncated at some point.It is worth noting that the result seems to defy our physical intuition.With the increase of the scale factor and the annihilation of photons, the entropy and energy of photons inside the co-moving volume will be reduced synchronously.However, the reduction of the energy is faster than that of the entropy, so the critical condition of the Bekenstein bound (i.e., S = 2πk ℏc ER) will appear and then be broken as t increases.It is different from the critical condition of the Bekenstein bound for compact objects, which only occurs in the case of the black hole and will not be violated as the mass of the black hole changes.The difference stems from the selection of the volume of the system.The former is the co-moving volume of the universe, which is affected by the background space-time.But, the latter is the volume of a sphere that can enclose the given system, which is related to the total mass-energy of the system.With regard to the issue, it is not the topic of this work, so we do not discuss it any further here.
Next, we consider the spherical entropy bound, which is given as Here, A is the area of the system and l p is the Planck length.For the photons inside the co-moving volume (A = a 2 ), the spherical entropy bound can be expressed as When Γ = 0, since the scale factor is monotonically increasing, once a 2 ≥ 8π 4 l 2 p 45ζ(3) N 0 , the spherical entropy bound will not be violated.Therefore, a = is the lower bound of the scale factor, which means that the spherical entropy bound could avoid the cosmological singularity in a radiation-dominated universe without particle production.
When Γ ̸ = 0, using Eq. ( 19) to eliminate Γ, we can obtain In order to judge whether this inequality is true, we need to figure out the relationship between a and T by solving Eq. ( 19).Here, we can discuss some issues qualitatively in the absence of the solutions of Eq. (19).Reviewing inequality (22), one can find that if Γ = 2 a da dt > 0, we have exp If a 0 = 1 (m) and T 0 ∼ 2.7 K, this inequality is true and independent of the evolution of the scale factor.We can set Γ c = 2 a da dt as a critical particle production rate.Note that Γ c is a function of time not a point.
for t > t 0 .As long as t is large enough, inequality (22) will be violated.Therefore, to meet the spherical entropy bound there should be a truncation for the production of photons at some point.Moreover, for Γ > Γ c > 0, as we mentioned earlier, there is naturally no cosmological singularity due to the solution of the scale factor, but the spherical entropy bound might modify the initial volume of the universe.
When Γ < Γ c and t > t 0 , we have aT a 0 T 0 3 < a 2 a 2 0 and then inequality (22) must be true.Therefore, the spherical entropy bound can not constrain the production (annihilation) rate of photons.But, if Γ < Γ c and t < t 0 , we have aT and then there could exist a lower bound for the scale factor, which can avoid the cosmological singularity.
In this section, we have not provided the solutions of Eq. ( 19), so there is only qualitative analysis.Next, we will discuss some similar issues in detail in the presence of analytical solutions.

Particle horizon
We will now take the volume covered by the particle horizon of a given observer as the thermodynamic system satisfying entropy bounds.For convenience, the production rate of photons inside the particle horizon will be characterized by the production rate of photons inside the co-moving volume.In an identical universe, the difference between the two particle production rates depends on the evolution of the universe.Therefore, the subsequent calculations are also characterized by the production rate of photons inside the co-moving volume.
We start from the solution of the energy density of photons in the presence of particle production.According to the laws of thermodynamics for open systems, the energy density of photons satisfies [45][46][47] Therefore, the solution of the energy density of photons can be expressed as For Γ = 0 and a 0 = 1, it degenerates into the standard solution ρ = ρ 0 a −4 .Since the universe is dominated by photons, we can ignore other matter when we solve the Friedmann equations ¶ .Thus, taking Eq. ( 27) into the Friedmann equations, we finally get where G 0 = 8 3 πG with G the Newtonian constant.When Γ = 0, we have a ∼ t 1/2 .As for Γ ̸ = 0, to get the analytical solution of the scale factor, we have to presuppose a specific form of Γ.We can analyze qualitatively the properties of the solution of the scale factor for different Γ's.If Γ = 0 results in a = b t 1/2 (where a(t 0 ) = a 0 and a(0) = 0), it can be expected that Γ > 0 corresponds to a > b t 1/2 (where a(t 0 ) = a 0 and a(0) > 0) and Γ < 0 corresponds to a < b t 1/2 (where a(t 0 ) = a 0 and a(0) < 0).Note that a(0) > 0 (i.e., Γ > 0) means that the beginning (t = 0) of the universe is not a singularity [45][46][47].Moreover, a(0) < 0 (i.e., Γ < 0) indicates that the singularity (i.e., a(t) = 0) of the universe appears at the non-beginning (t ̸ = 0) of the universe (see Fig. 1).In this section, we must be cautious about these two strange situations because the particle horizon may be unnormal, which should be avoided.Next, we discuss these issues in detail.
The particle horizon of a given observer in the radiation-dominated universe is given as ¶ Note that the reason why the number of photons inside the co-moving volume is not conserved is the interaction between photons and other matter (or the background space-time).Therefore, ignoring the matter (or the background space-time) interacting with photons is a rough approximation.However, as a toy model, we actually only need an analytical solution of the scale factor increasing monotonically over time.Incorporating the matter (or the background space-time) interacting with photons into the Friedmann equations will only make it more difficult to solve the equations analytically and is not helpful for the following research.
where t s ≥ 0 is the time at which the observer starts observing at r = 0. Since a(t) > 0, R H is a monotonically increasing function of t.For a = b t 1/2 , we have R H = 2b −1 c t 1/2 − t 1/2 s .In Ref. [26], based on the particle horizon of a given observer in a radiation-dominated universe and the Bekenstein bound, Bekenstein pointed out that the cosmological singularity is not thermodynamically possible (see Refs. [26,27] for more details).Now, let us analyze how Γ affects the cosmological singularity.For a system covered by the particle horizon of a given observer, the Bekenstein bound can be expressed as When Γ = 0, there is a lower bound for the scale of the universe determined by the Bekenstein bound [26,27].Since Γ influences the solution of the scale factor (which determines the particle horizon and the temperature of photons), it may modify the lower bound of the scale of the universe.
From the previous analysis, a positive-definite Γ will lead to a bigger scale factor (the scale factor given by Eq. ( 28) should satisfy a > b t 1/2 ) and a smaller particle horizon . However, a positive-definite Γ will also cause the temperature of photons to be higher.Therefore, we can not easily estimate the effect of Γ on the lower bound of the scale of the universe on the basis of inequality (30).With a specific Γ, we can analyze the issues further.
We set Γ = g t , where g > 0 (g < 0) represents particle production (annihilation).Then, taking it into Eq.( 28), the solution of the scale factor can be expressed as where b 1 = 2(G 0 ρ 0 ) 3a 2 0 3+2g t 0 + a 2 0 .We consider that g > 0, g = 0, and g < 0 correspond to d 1 > 0, d 1 = 0, and d 1 < 0, respectively.From Fig. 1, one can find that when Γ < 0, the value of the scale factor at t = 0 will be negative by extending the green dash-dotted line, which may lead to confusion when we calculate the particle horizon.In order to avoid such a nuisance, the lower bound of the integral in Eq. (29) should not be smaller than the minimum value of time (which ensures that the scale factor is non-negative).
Then, the particle horizon can be expressed as where 2 F 1 (a, b, c, z) is a hypergeometric function.For convenience, we can assume that, for any Γ, we always have F (t) = c a(t) dt and F (t s ) = 0, which will not alter our subsequent analysis and conclusion.Then, the second line in the above equation is vanishing.
With Γ = g t , ρ = π 2 k 4 15c 3 ℏ 3 T 4 , Eqs. ( 27) and ( 30), the temperature of photons is given as Reviewing inequality (30), to study the influence of Γ on the cosmological singularity, we only need to figure out the impact of the parameter g on T R H . Labelling all the positive coefficients in T R H as a ) marked with the green dash-dotted line, respectively.These three lines intersect at the point (t 0 , a 0 ).When Γ > 0 and t = 0, the scale factor is larger than zero.When Γ < 0, the minimum value of time is larger than zero.positive-definite parameter M 0 , T R H can be written as where M 0 = c . For a given g, T R H generally can be simplified.However, since the hypergeometric function is complicated, for a given g, it is more reasonable to calculate T R H starting from Eq. (31).Next, we consider the three cases of g adopted in Fig. 1.
is a constant (see the red dashed line in Fig. 2).In this case, one can prove that inequality (30) is tenable if a 0 = 1 (m) and T 0 ∼ 2.7 K, so the cosmological singularity can not be avoided.Note that t s = 0 just corresponds to a special observer.If there are other observers who deny the rationality of the cosmological singularity, the universe should not have the initial singularity.When Γ = 0 and t s > 0, one can find that T R increases with t and has a maximum value.When t → t s , we have T R H → 0, which inevitably causes inequality (30) to fail.Therefore, for a given t s , t must satisfy We set t c as the critical time corresponding to the equality sign of the above inequality.Then, a(t c ) characterizes the critical scale (i.e., lower bound of the scale) of the universe determined by the Bekenstein bound.If t s can approach to zero, the critical value (i.e., lower bound) of the scale factor can also approach to zero.However, if there is a minimum scale for time (such as the Planck-Wheeler time), then a(t c ) can not approach to zero [26,27].We will not dig and delve what value t s should take here.We just show that the Bekenstein bound could require a (nonzero) critical value for the scale factor in the radiation-dominated universe without particle production (Γ = 0), so the cosmological singularity can be avoided.When Γ ̸ = 0, in order to study the effect of Γ on the critical value of the scale factor, we first analyze the critical time for Γ ̸ = 0.For convenience, we set t s as the minimum value of time (see Fig. 1).Note that here the minimum value of time depends on the solution of the scale factor and has no connection with the Bekenstein bound.The "minimum value of time" determined by the Bekenstein bound is called critical time.We plot the schematic diagram of T R H in Fig. 2. From the left panel in Fig. 2, we can find that when t ≪ t 0 , if there exists production of photons, T R H (see the black solid line) will be smaller than the one in the case of Γ = 0 (see the red dashed line).According to inequality (30), if the critical time for Γ = 0 is t c , then the critical time for Γ > 0 would be t c1 > t c .Moreover, from Fig. 1, the scale factor for Γ > 0 is larger than that for Γ = 0, so t c1 > t c means that the critical value of the scale factor for Γ > 0 is larger than that for Γ = 0.In other words, the production of photons causes the critical scale of the universe to be larger.As we mentioned earlier, for Γ > 0, the solution of the scale factor naturally removes the cosmological singularity (i.e., a(t = 0) > 0, see Fig. 1).Here, the lower bound of the scale factor determined by the Bekenstein bound is unquestionably larger than a(t = 0).As for Γ < 0, it is difficult to tell what the critical time is, because it could be t c2 < t c or may not exist.The former means that the critical value of the scale factor for Γ < 0 is smaller than that for Γ = 0, and the latter means the cosmological singularity can appear.2 ) marked with the green dash-dotted line, respectively.
Next, we analyze the constraint of the Bekenstein bound on the production of photons.One can find, from the right panel in Fig. 2, that with the increase of t, whatever the value of Γ is, T R H is always increasing.Only when Γ = 0 and t s = 0, T R H is a constant, for which the Bekenstein bound is tenable.Recalling inequality (30), the larger T R H , the easier it is satisfied.Since T R H increases with t, the Bekenstein bound usually will not be violated in the late universe.Therefore, when we consider the entropy of photons in a system covered by the particle horizon in the late radiation-dominated universe, the Bekenstein bound can not provide effective constraint on the production of photons (i.e., the Bekenstein bound is easily satisfied).The interaction between photons and other matter (or the background space-time) does not need to be truncated.Now, we consider the spherical entropy bound for a system covered by the particle horizon of a given observer.We still assume that a 0 = 1 (m), T 0 ∼ 2.7 K, and t s is the minimum value of time.With Eq. ( 27) and ρ = π 2 k 4 15c 3 ℏ 3 T 4 , the spherical entropy bound can be written as where M 0 = c , inequality (36) can be simplified as Since t > t s ≥ 0, for a given nonzero t s , when t → t s , the left-hand side of the above inequality will approach to infinity.As t increases, it reaches the minimum at t = 9 4 t s and then monotonically increases to infinity as t → ∞.In this case, it is not easy to determine the lower bound of time (critical time) according to inequality (37).However, when t s = 0, we can find that inequality (37) requires a nonzero lower bound for time (see the red dashed line in Fig. 3).As long as there exists an observer requiring a nonzero lower bound for time, the lower bound of time must be nonzero.Therefore, considering the spherical entropy bound in a system covered by the particle horizon in the radiation-dominated universe, the cosmological singularity can be avoided.
When Γ ̸ = 0, we can still set Γ = g t .Based on the previous calculations, if g = 1 2 , inequality (36) can be reduced as where R H is given by Eq. (32) with t s = 0 (the minimum value of time, see Fig. 1).If g = − 1 2 , inequality ( 36) is reduced as where R H is given by Eq. ( 32) with t s = 3 −3/2 (the minimum value of time, see Fig. 1).
In order to compare the critical times (the lower bounds of time) for the three cases, we can set = 1 and t 0 = 1.Then, we label uniformly the left-hand sides of these three inequalities as a function f (t).The evolutions of f (t) over time for the three cases are plotted in Fig. 3. From Fig. 3, one can find that when there exists production of photons (Γ > 0, see the black solid line), f (t) decreases monotonically with time, so the corresponding critical time does not exist.In other words, even if t → 0, inequality (38) still holds.Therefore, in this case the spherical entropy bound can not modify the critical scale of the universe (which depends on the solution of the scale factor).On the other hand, since f (t) decreases monotonically with time, there should be a truncation for the production of photons (see the black solid line).When there exists annihilation of photons (Γ < 0, see the green dash-dotted line), f (t) increases monotonically with time.We find that the critical time for such case is smaller than the one obtained by Γ = 0.According to Fig. 1, the scale factor for Γ < 0 is smaller than that for Γ = 0, so the critical scale of the universe for Γ < 0 must be smaller than that for Γ = 0.Moreover, since f (t) increases monotonically with time, the spherical entropy bound can not constrain the annihilation of photons in the radiation-dominated universe.It is worth mentioning that these results are based on t s being the minimum value of time.As for other t s 's, there may be different results, which we will not continue to discuss in this work.) marked with the green dash-dotted line, respectively.

Entropy bounds and particle production in a dust-dominated universe
In the late 1980s, Prigogine proposed that cosmological particle production can avoid the initial singularity and solve the entropy problem [45][46][47].In his cosmological model, the entropy of the universe is expressed as the product of the particle number and the specific entropy (the average entropy of a single particle).If there exists particle production in the universe, the entropy inside the co-moving volume will keep increasing and so the entropy problem can be solved.In this section, we consider a similar toy cosmological model, i.e., a dust-dominated universe.The entropy of dust in any system can be expressed as where σ(t) is the specific entropy of dust and Γ is the production rate of dust in the system.Next, we study the effect of the production of dust on the cosmological singularity in light of entropy bounds and the constraint of entropy bounds on the production of dust.We still focus on the co-moving volume and the volume covered by the particle horizon of a given observer.

Co-moving volume
For the co-moving volume, we have N 0 = n 0 a 3 0 in Eq. ( 40) and Γ is the production rate of dust inside the co-moving volume.Comparing Eq. ( 40) with Eq. ( 14), one can find that if σ(t) is a constant, the entropy evolution of the dust-dominated universe is similar to that of the radiation-dominated universe.Therefore, an evolving σ(t) has more general properties in characterizing the entropy evolution of the universe.Moreover, the solutions of the Friedmann equations for the dust-dominated universe are different from the ones for the radiation-dominated universe, which could also lead to different conclusions for similar situations.
We first consider the Bekenstein bound, which, for the dust inside the co-moving volume, is given as For convenience, the potential energy of dust can be integrated into the rest mass of a single particle.In this case, the energy of dust inside the co-moving volume can be written as where m is the equivalent mass of a single particle.The lower bound of m corresponds to the rest mass of a single particle.Since the temperature of the universe is changing continuously, m also varies with time.Theoretically, as a → 0, m could be infinite.So the Bekenstein bound is reduced as If σ(t) is a constant, the Bekenstein bound must be robust in the late dust-dominated universe due to the nonzero lower bound of m.Therefore, the production (annihilation) of dust in the late universe will not be limited by the Bekenstein bound.In addition, whether there is a lower bound for the scale factor relies on the evolutions of m, σ(t) and the scale factor itself.If inequality (43) does not hold when a is less than a certain threshold, then the threshold is the lower bound of the scale factor determined by the Bekenstein bound, which avoids the cosmological singularity.But, if inequality (43) holds for any a, the cosmological singularity can not be avoided.
For a general σ(t), since it can be independent of the scale factor, it is convenient to absorb m into σ(t) as a new parameter.Next, we set σ m (t) = σ(t)/m signifying the entropy of dust per unit mass, so inequality (43) can be re-expressed as Note that although there is no apparent Γ in the inequality, the solution of the scale factor depends on Γ.Therefore, Γ and σ m (t) jointly determine the cosmological singularity, and Γ can be constrained by the Bekenstein bound for a given σ m (t).With an analytical solution of the scale factor, we can discuss these issues in detail.
Since we have absorbed the potential energy of dust into its rest mass, the (equivalent) pressure of dust can be approximately equal to zero and the (equivalent) energy density of dust satisfies The solution of the above equation can be expressed as Taking it into the Friedmann equations, one can get the solution of the scale factor: Reviewing inequality (44), in order to obtain the upper bound (which is related to the constraint on the production of dust) and the lower bound (which is related to the cosmological singularity) of the scale factor, we need to presuppose specific forms of σ m (t) and Γ.We take σ m (t) = (p t) 2 3 (p > 0 to guarantee that the entropy of the system is always increasing) and Γ = g t as an example to illustrate the issues in detail.In this case, inequality (44) can be expressed as For convenience, we can set 2πkc ℏ = 1, 3 2 (G 0 ρ 0 ) 1 2 = 1, a 0 = 1, and t 0 = 1.Then, the Bekenstein bound is reduced to We can set the right-hand side of the above inequality as a new function G(t), which is plotted in Fig. 4 with some representative values of the parameters (g, p).From Fig. 4, we can judge whether the universe has the initial singularity for different situations and obtain the constraint of the Bekenstein bound on the production of dust.For different values of the parameters (g, p), the results are totally different.
When g = 0 (Γ = 0), as long as p ≤ 1, the Bekenstein bound is always valid (see inequality (49) and the black dotted line), so it can not avoid the cosmological singularity.For p > 1, the Bekenstein bound does not hold at any time (see inequality (49) and the red dotted line), so this situation can not happen inside the co-moving volume of the dust-dominated universe.
When g = 2 (Γ > 0) and p = 3 (see the red solid line), since the Bekenstein bound holds at the beginning (t = 0) of the universe, it can not avoid the cosmological singularity.Note that for Γ > 0 the scale factor is naturally nonzero at the beginning (t = 0) of the universe.So, there is still no initial singularity and the Bekenstein bound can not modify the critical scale of the universe.But over time, the production of dust will break the Bekenstein bound, so it requires a truncation for the interaction between dust and other matter (or the background space-time).When g = 2 (Γ > 0) and p = 3 5 (see the black solid line), the Bekenstein bound is tenable.Therefore, it can not avoid the cosmological singularity and constrain the production of dust.
When g = −3 (Γ < 0) and p = 3 (see the red dashed line), the Bekenstein bound is invalid at any time, so it is an unreasonable case.When g = −3 (Γ < 0) and p = 3 5 (see the black dashed line), there exists a nonzero critical time.However, a nonzero critical time does not signify a nonzero critical scale of the universe.For Γ < 0, according to the solution of the scale factor (see Fig. 1), only when the nonzero critical time is larger than the minimum value of time, the cosmological singularity can be avoided.One can take the nonzero critical time into the solution of the scale factor and then it is found that the scale factor is negative, which means that the Bekenstein bound can not prevent the appearance of the cosmological singularity.With the increase of t, there will be a new intersection between the black dashed line and the horizontal axis, which requires the annihilation of dust to be cut off in the late universe.
Here we can briefly summarize the above discussion as follows.When Γ ≤ 0, the Bekenstein bound can not avoid the cosmological singularity, but there could exist a truncation for the annihilation of dust.When Γ > 0, there is no cosmological singularity due to the solution of the scale factor, and the Bekenstein bound can not affect the critical value of the scale factor (the critical scale of the universe).As for whether there is a truncation for the production of dust, it depends on the values of the parameters (g, p).All of these results are based on σ m (t) = (p t) Fig. 4: Plot of G(t).There are six sets of the parameters (g, p): g = 2 (Γ > 0) and p = 3 marked with the red solid line; g = 2 (Γ > 0) and p = 3 5 marked with the black solid line; g = 0 (Γ = 0) and p = 3 marked with the red dotted line; g = 0 (Γ = 0) and p = 3  5 marked with the black dotted line; g = −3 (Γ < 0) and p = 3 marked with the red dashed line; g = −3 (Γ < 0) and p = 3  5 marked with the black dashed line.
If we use the spherical entropy bound to discuss the cosmological singularity and the constraint on the production of dust inside the co-moving volume, we need the entropy of dust inside the co-moving volume to satisfy We can still set σ(t) = (p t) 2 3 (p > 0) and Γ = g t .Note that the spherical entropy bound does not involve the equivalent mass of a single particle, so we do not have to define the entropy of dust per unit mass.The spherical entropy bound can be expressed as For convenience, we can set k Similarly, we set the right-hand side of the above inequality as a new function G 1 (t), which is plotted in Fig. 5 with some representative values of the parameters (g, p).Comparing Fig. 5 with Fig. 4, one can find that when Γ > 0 (see the solid lines), the results are similar.Therefore, we will not rehash the discussion on the case.When Γ = 0 (see the dotted lines), there exists a nonzero critical time (0 < p 1/2 ≤ t c ), which can be obtained directly from inequality (52) with g = 0. Since the evolution of the scale factor starts from a(t = 0) = 0 for Γ = 0, the cosmological singularity can be avoided, which is different from the previous case (see the dotted lines in Fig. 4).Apparently, as the parameter p grows, the critical time will be larger and so the critical scale of the universe will be larger.When Γ < 0 (see the dashed lines), there exists a nonzero critical time.But, the nonzero critical time usually can not avoid the cosmological singularity because a(t = 0) < 0 for Γ < 0. Only when the critical time is larger than the minimum value of time, the cosmological singularity can be avoided.In addition, there does not exist a truncation for the annihilation of dust.Again, these results are based on σ(t) = (p t) Fig. 5: Plot of G 1 (t).There are six sets of the parameters (g, p): g = 2 (Γ > 0) and p = 3 marked with the red solid line; g = 2 (Γ > 0) and p = 3 5 marked with the black solid line; g = 0 (Γ = 0) and p = 3 marked with the red dotted line; g 0 (Γ = 0) and p = 3  5 marked with the black dotted line; g = −3 (Γ < 0) and p = 3 marked with the red dashed line; g = −3 (Γ < 0) and p = 3  5 marked with the black dashed line.

Particle horizon
At last, we consider the volume covered by the particle horizon of a given observer in the dustdominated universe.The solution of the scale factor is also given by Eq. (47).We still set Γ = g t , c = 1, 3 2 (G 0 ρ 0 ) 1 2 = 1, a 0 = 1, and t 0 = 1.So, the particle horizon can be expressed as where 2 F 1 (a, b, c, z) is a hypergeometric function.According to our previous assumption, for any Γ, if F (t) = c a(t) dt and F (t s ) = 0, then the Bekenstein bound (replacing the scale factor in inequality ( 44) with R H ) can be expressed as 0 ≤ (g + 2)t g Here, we have set 2πkc ℏ = 1 and σ m (t) = (p t) 2 3 (p > 0).The right-hand side of the above inequality is defined as a new function G 2 (t), which is plotted in Fig. 6 with some representative values of the parameters (g, p).In all cases, t s is the corresponding minimum value of time.
From Fig. 6, we can find that when Γ > 0 (see the solid lines), there are two possible consequences.For a smaller p (such as g = 2 and p = 3 5 , the black solid line), there exists a truncation for the production of dust due to the Bekenstein bound (see the black solid line in the right panel), and the cosmological singularity can be avoided by the Bekenstein bound (see the black solid line in the left panel).In this case, the Bekenstein bound will make the critical scale of the universe bigger.For a larger p (such as g = 2 and p = 3, the red solid line), the Bekenstein bound can not be satisfied, so this situation can not happen inside the volume covered by the particle horizon.When Γ = 0 (see the dotted lines), there is only one possible consequence.The entropy of dust inside the volume covered by the particle horizon is consistent with the Bekenstein bound at the beginning (t = 0) of the universe, and then it will break the Bekenstein bound at some point (i.e., the special critical time).Therefore, in this case, the Bekenstein bound can not avoid the cosmological singularity but it provides an upper bound for time (i.e., the special critical time).In this work, we can not offer a definitive explanation for the upper bound of time.When Γ < 0 (see the dashed lines), it is found that the critical time is larger than the minimum value (t s = 4 9 ) of time (at which the scale factor is vanishing).Therefore, the cosmological singularity can be avoided.Comparing the critical time for p = 3 (see the red dashed line in the right panel) and the one for p = 3 5 (see the black dashed line in the left panel), one can find that as p increases, the critical scale of the universe will be larger.In addition, we can find that the annihilation of dust can not be truncated in the late universe.
For the spherical entropy bound, the entropy of dust inside the volume covered by the particle horizon of a given observer needs to satisfy In this case, we do not need to define the entropy of dust per unit mass.Note that since Γ is the particle production rate of dust inside the co-moving volume, there is an extra dimensionless factor  Fig. 6: Plot of G 2 (t).There are six sets of the parameters (g, p): g = 2 (Γ > 0) and p = 3 marked with the red solid line; g = 2 (Γ > 0) and p = 3 5 marked with the black solid line; g = 0 (Γ = 0) and p = 3 marked with the red dotted line; g = 0 (Γ = 0) and p = 3 5 marked with the black dotted line; g = 0 (Γ = 0); g = −3 (Γ < 0) and p = 3 marked with the red dashed line; g = −3 (Γ < 0) and p = 3 5 marked with the black dashed line.
The right-hand side of the above inequality is defined as a new function G 3 (t), which is plotted in Fig. 7 with some representative values of the parameters (g, p).In all cases, t s is still the corresponding minimum value of time.
From the upper left panel in Fig. 7, it can be found that when Γ > 0, there exists a critical time and a truncation for the production of dust due to the spherical entropy bound.But the critical time is later than the time corresponding to the truncation, which seems strange.Here, we can deem that the truncation is non-physical and the universe began at the critical time.In this case, the universe has no initial singularity and the production of dust can not be cut off in the late universe.But, if we accept the existence of the truncation, then the universe originated from the minimum value (t = 0) of time, but has no initial singularity due to the solution of the scale factor.Therefore, the spherical entropy bound can not modify the critical scale of the universe.The production of dust can only occur in the early universe and beyond the critical time.According to the upper right panel and lower panel in Fig. 7, when Γ ≤ 0, the cosmological singularity can be avoided by the spherical entropy bound.And there is no truncation for the annihilation of dust.
Finally, we emphasize again that these results are based on the specific assumptions about the model and the parameters.The results may be completely different for other models and parameters.

Conclusions and Discussions
The Bekenstein bound has been proposed more than 30 years, after which there emerge multiple definitions of entropy bounds.Although most of these entropy bounds are established on black hole research, their applications in cosmology are also widely studied, especially the topics related to the initial singularity and entropy of the universe.The entropy of the universe caused by particle production is usually constrained by the (generalized) second law of thermodynamics and the law of thermal equilibrium, but the results in practice are not satisfactory [49,57].In this work, we use entropy bounds (the Bekenstein bound and the spherical entropy bound) to constrain the entropy of particles in two toy cosmological models with particle production, thereby constraining the production of the corresponding particles.We also study the effect of particle production on the cosmological singularity in light of entropy bounds.For the two toy cosmological models we study, the cosmological singularity can be avoided by considering entropy bounds in some special cases, such as the system covered by the particle horizon in the radiation-dominated universe with Γ > 0 constrained by the Bekenstein bound (see Fig. 2).In addition, the production (annihilation) of particles needs to be truncated in some special cases due to entropy bounds, such as the co-moving volume in the dust-dominated universe with Γ < 0 and p = 3 5 constrained by the Bekenstein bound (see Fig. 4).There are also some cases that always do not satisfy certain entropy bound, and thus they can not happen in the corresponding cosmological model, such as the system covered by the particle horizon in the dust-dominated universe with Γ > 0 and p = 3 constrained by the Bekenstein bound (see Fig. 6).Moreover, some cases always satisfy certain entropy bound, which means that the corresponding entropy bound does not play much of a role, such as the co-moving volume in the dust-dominated universe with Γ > 0 and p = 3 5 constrained by the spherical entropy bound (see Fig. 5).In general, for different situations there are various possible results about the cosmological singularity and the truncation of particle production (annihilation).It should be emphasized that these results are mainly dependent on the selections of the system in the universe and the entropy bound.
Since most entropy bounds are obtained from black hole research, they may not apply to the universe.Therefore, it is worth discussing whether we can employ directly the entropy bound obtained from black hole research to cosmology.Seeking an entropy bound which is really practical in cosmology is the focus of our next research, which can also avoid the emergence of various possible results for different entropy bounds.To sum up, our methods are basically appropriate for all cosmological models with particle production, but it is still a preliminary attempt.How to implement the researches in the real universe and how to choose the system in the universe and the entropy bound need to be further studied in the future.

2 3
(p > 0) and Γ = g t .For other forms of σ m (t) and Γ, the results could be completely different.Truncation for the production of dustTruncation for the annihilation of dust

R 3 Ha 3
on the left-hand side of the above inequality.Similarly, we set σ(t) = (p t)

Fig. 7 :
Fig.7: Plot of G 3 There are six sets of the parameters (g, p).The upper left panel corresponds to g = 2 (Γ > 0), where p = 3 is marked with the red solid line and p =3  5 is marked with the black solid line.The upper right panel corresponds to g = 0 (Γ = 0), where p = 3 is marked with the red dotted line and p =3  5 is marked with the black dotted line.The lower panel corresponds to g = −3 (Γ < 0), where p = 3 is marked with the red dashed line and p =3  5 is marked with the black dashed line.