HTML
-
Accelerated expansion of the universe has been confirmed via observations by the Type Ia supernovae [1, 2], cosmic microwave background radiation [3, 4] and baryon acoustic oscillations [5, 6] since the 1990s. To explain this phenomenon, the concept of dark energy was introduced, and various models were proposed. From the perspective of the Einstein field equation, there are two main approaces to describe dark energy. One is to modify the matter sector of the Einstein field equation. These proposed models include the lambda cold dark matter (
$ \Lambda\mathrm{CDM} $ ) model [7],$ \Lambda(t)\,\mathrm{CDM} $ (i.e. the dark energy varies with time t) model [8, 9], particle creation model [10, 11], etc. Another is to modify the geometrical sector of the Einstein field equation. These theories are also referred to as the modified gravity, for example Lanczos-Lovelock gravity [12],$ f(R) $ gravity [13, 14], etc. Irrespective of the approach, however, one of the main goals of these models is to describe the current accelerated expansion of the universe [15, 16]. Among them, the particle creation model has some of the following advantages [11, 17]: (i) A description of the non-equilibrium thermodynamics exists, as the process of particle creation is an irreversible process accompanied by entropy production [18, 19]. (ii) The particle creation mechanism unifies the dark energy and dark matter, as it contains only a single free parameter, namely the particle creation rate. The one-parameter model is preferred by statistical Bayesian analysis [20]. In Ref. [21], the authors found that the bulk viscous model (the one-parameter model) has some advantages over the$ \Lambda\mathrm{CDM} $ model through the Bayesian analysis of Supernovae Type Ia data. In particular, the particle creation model can be described by the gravitationally induced particle creation mechanism. The particle number was shown not to represent a constant of the motion, even though the gravitational field is not quantized considering the covariant quantized free field equations of elementary particles in an expanding universe [22-24]. In addition to the above microscopic description, Prigogine et al. provided a macroscopic description of the particle creation mechanism induced by the gravitational field [10].In the past 30 years, significant progress has been made in the study of gravity from the viewpoint of thermodynamics. In 1995, Jacobson showed that the Einstein field equation can be derived based on the Clausius relation and entropy-area relation [25]. Subsequently, Clausius and his collaborators investigated the properties of non-equilibrium thermodynamics of spacetime by introducing the term of entropy generation [26]. In contrast, Padmanabhan derived the important thermodynamic relation
$ S = 2\beta E $ [27] and the first law of thermodynamics [28], which are consistent with the conventional form of thermodynamics from a similar but different thermodynamic perspective. Further, there is a simple correspondence between the surface$ L_{\rm sur} $ and bulk$ L_{\rm bulk} $ terms when the Lagrangian of gravity L is decomposed into$ L = L_{\rm sur}+L_{\rm bulk} $ [29, 30], which shows that the gravitational action is holographic, as the same information is encoded in both terms. These results indicate a deep connection between gravitational dynamics and horizon thermodynamics.In particular, Padmanabhan proposed that the spacetime is composed of microscopic degrees of freedom, and the de Sitter universe satisfies the holographic equipartition (an especial holographic principle) [31, 32]. Further, he thought that our universe is asymptotically de Sitter and derived the evolution equations of the universe based on the difference between the number of degrees of freedom on the surface and those in the bulk. This model was subsequently extensively studied [33-41]. Therefore, there is a solid physical foundation to study the evolution of the universe from the perspective of holographic equipartition.
In this paper, we propose a cosmological scenario that describes the evolution history of the universe based on particle creation and holographic equipartition. Throughout the evolution of the universe, we assume that the universe consistently creates particles, and holographic equipartition is always satisfied. In our scenario, the creation rate of the radiation is
$ \Gamma = \alpha H^{2} $ in the early universe, and the creation rate of the pressureless matter is$ \Gamma = \alpha H $ in the present and late universe, where$ \alpha $ is a positive parameter, and H is the Hubble parameter. The entire evolution history of the universe may be explained as follows. The universe starts from an unstable de Sitter universe ($ a(t)\propto {\rm e}^{\alpha t/3} $ ) and evolves into a standard radiation stage ($ a(t)\propto t^{1/2} $ ) because of the creation of radiation. With the expansion of the universe, the pressureless matter, whose creation rate is$ \Gamma = \alpha H $ , begins to dominate the universe. The negative creation pressure of the matter accelerates the expansion of the universe and drives the present accelerated universe ($ a(t)\propto t^{\delta}\,(\delta>1) $ ) to the de Sitter universe ($ a(t)\propto {\rm e}^{Ht} $ ).This paper is organized as follows. In Section 2, we introduce non-equilibrium thermodynamics, which describes particle creation. In Section 3, we obtain the energy of the universe enclosed by the Hubble horizon using the laws of energy and holographic equipartition. Then, we obtain the laws of evolution for the early universe as well as the present and late universe by choosing suitable creation rates of particles in Section 4. In Section 5, we investigate thermodynamical properties for this model. The summary is provided in Section 6. We employ units
$ {c = \hbar = 1} $ . -
The first law of thermodynamics in a closed system that only performs volume work can be expressed as
$ \begin{array}{l} {\rm d}E = \delta Q-p{\rm d}V, \end{array} $
(1) where
$ \delta Q $ is the amount of heat exchanged by the system; p is the pressure;$ {\rm d}E $ is the change of internal energy, and$ {\rm d}V $ is the change in volume. In a reversible thermodynamic process, the amount of heat exchanged by the system can be expressed as$ \delta Q = T{\rm d}S $ , where T is the temperature, and$ {\rm d}S $ is the entropy change.For an adiabatic open system, where the particle number is not conserved due to particle creation, the first law of thermodynamics is written as [10]
$ \begin{split} {\rm d}(\rho V)+p{\rm d}V = \frac{h}{n}{\rm d}(nV), \end{split} $
(2) where n is the particle number per unit volume;
$ h = \rho+p $ is the enthalpy, and$ \rho $ is the energy density.To calculate the creation pressure attributed to particle creation, we consider the universe as a sphere with a radius
$ a(t) $ . Thus, the volume is$ V = \frac{4\pi}{3}a^{3}(t) $ , and the right term of Eq. (2) is expressed as$ \begin{split} \frac{h}{n}{\rm d}(nV) = (\rho+p)V\left(3H+\frac{\dot{n}}{n}\right){\rm d}t, \end{split} $
(3) where
$ H = \dfrac{\dot{a}(t)}{a(t)} $ is the Hubble parameter. Inserting Eq. (3) and the relation$ \dot{n}+3Hn = n\Gamma $ [10, 42-44] into Eq. (2), where$ \Gamma $ is the particle creation rate, we obtain$ \begin{array}{l} \dot{\rho}+3H(\rho+p+p_{\rm e}) = 0, \end{array} $
(4) where the effect of particle creation is expressed as an extra pressure
$ \begin{split} p_{\rm e} = -\frac{(\rho+p)\Gamma}{3H}. \end{split} $
(5) Thus, we can deem that Eq. (4) is the modified continuity equation attributed to the particle creation.
Particle creation in the universe generates entropy, whose change is given by
$ {\rm d}S = \frac{s}{n}{\rm d}(nV) $ , where$ s = S/V $ [10], which provides an explanation for the origin of cosmological entropy. Furthermore, the negative pressure caused by particle creation can explain the current accelerated expansion. In fact, bulk viscous fluids are introduced to study the accelerated expansion of the universe also due to the negative pressure effect of viscous fluids [45-48].Although the particle creation model is widely employed in cosmology to investigate the accelerated expansion of the universe, and has made great progress therein, a question arises naturally on whether the effect of particle creation can be observed. Lima et al. [49, 50] performed kinematic tests for the flat, closed, and open FRW universes with matter creation by investigating the lookback time, age of the universe, luminosity distance, angular diameter and galaxy number counts redshift relations, and finally obtained results that agreed with the observation. In particular, research by Pigozzo et al. [51] indicated a cosmological, late-time dark matter creation at 95% confidence level by analyzing the linear matter power spectrum, distance measurements from Type Ia supernovae, and the position of the first peak in the anisotropy spectrum of the cosmic microwave background.
-
Research in recent decades has uncovered an important fact, namely that spacetime can be heated up like matter. Unruh discovered that an observer can measure a temperature
$ T = \frac{\kappa}{2\pi} $ when they accelerate through the inertial vacuum with a proper acceleration$ \kappa $ in 1976 [52]. Furthermore, the temperature measured by the accelerated observer in the vacuum is as real, as that measured by an inertial thermometer in ordinary matter [53]. Hence, it can be concluded that the spacetime, just as ordinary matter, is made up of microscopic degrees of freedom. In the description of the evolution of the universe based on the emergent perspective of gravity, the energy equipartition law$ \begin{split} E = \frac{1}{2}N_{\rm bulk}kT \end{split} $
(6) holds, where E is the energy in volume V,
$ N_{\rm bulk} $ is the number of degrees of freedom in volume V, and T is the local acceleration temperature.In Ref. [31], Padmanabhan proposed that our universe is asymptotically, rather than exactly, de Sitter, and the expansion of the universe is being driven by the difference (
$ N_{\rm sur}-N_{\rm bulk} $ ), where$ N_{\rm sur} $ is the number of surface degrees of freedom on the horizon. However, in this study, we assume that the universe always obeys holographic equipartition because of particle creation. Therefore, the number of bulk degrees of freedom enclosed by the horizon is$ \begin{array}{l} N_{\rm bulk} = N_{\rm sur}. \end{array} $
(7) In addition, the number of surface degrees of freedom can be expressed as [54]
$ \begin{split} N_{\rm sur} = \frac{A}{L_{p}^{2}}, \end{split} $
(8) where A is the area of the horizon. Hence, the total energy in the volume V is
$ \begin{split} E = \frac{A}{2L_{p}^{2}}kT. \end{split} $
(9) Considering a universe enclosed by the Hubble horizon, which is a sphere with the radius
$ H^{-1} $ , the volume, temperature, and area of the Hubble horizon can be written as$ V = \frac{4\pi}{3H^{3}}, $ $ kT = \frac{H}{2\pi} $ , and$ A = \frac{4\pi}{H^{2}} $ , respectively. Inserting these physical quantities into Eq. (9), we obtain$ \begin{split} E = \frac{1}{HL_{p}^{2}}. \end{split} $
(10) This is an explicit relation between the energy enclosed by the Hubble horizon and the Hubble parameter. The relation of Eq. (10) implies that the universe enclosed by the Hubble horizon has a negative specific heat, as
$ E\propto T^{-1} $ . This is also the expected result, even in Newtonian gravitating systems [55].There are two relevant important energy contents in the flat universe, namely the Misner-Sharp energy and Komar energy. Thus, a question naturally arises as to whether the energy E in the energy equipartition law represents the Misner-Sharp energy or the Komar energy. We followingly investigate the two energy conditions.
-
The Misner-Sharp energy is [56]
$ \begin{array}{l} E = \int T_{\mu\nu}u^{\mu}u^{\nu}{\rm d}V, \end{array} $
(11) where
$ u^{\mu} = \delta_{0}^{\mu} $ is the four-velocity, and$ T_{\mu\nu} = (\rho,\: p+p_{\rm e},\: $ $p+p_{\rm e},\: p+p_{\rm e}) $ is the energy-momentum tensor. Then, the energy inside the Hubble horizon is$ \begin{split} E = \int_{0}^{H^{-1}}4\pi r^{2}\rho {\rm d}r = \frac{4\pi\rho}{3H^{3}}, \end{split} $
(12) where we have assumed that the energy density
$ \rho $ is homogeneous. Comparing Eq. (10) with Eq. (12), we obtain$ \begin{split} \rho = \frac{3H^{2}}{4\pi L_{p}^{2}}. \end{split} $
(13) Thus, we obtain one of the evolution equations of the universe from the energy equipartition law. The equipartition law is an equation of state, while Eq. (13) can be deemed as a dynamic evolution equation of the universe. Thus, the evolution of the universe can be derived from an equation of state, namely the equipartition law.
-
The Komar energy is defined as [57, 58]
$ \begin{array}{l} E = \int(2T_{\mu\nu}-Tg_{\mu\nu})u^{\mu}u^{\nu}{\rm d}V, \end{array} $
(14) where T is the trace of the energy-momentum source
$ T_{\mu\nu} $ . Then, the energy inside the Hubble horizon can be reduced to$ \begin{split} E = \int_{0}^{H^{-1}}4\pi r^{2}|\rho+3(p+p_{\rm e})|{\rm d}r = \frac{4\pi}{3H^{3}}|\rho+3(p+p_{\rm e})|. \end{split} $
(15) Thus, Eq. (10) changes to
$ \begin{split} |\rho+3(p+p_{\rm e})| = \frac{3H^{2}}{4\pi L_{p}^{2}}. \end{split} $
(16) This demonstrates that the change of the Hubble parameter is related to the creation pressure of particles when the Komar energy is chosen. Interestingly, Padmanabhan suggested that the Komor energy is the active gravitational mass-energy [27] and derived the standard Friedmann equation taking into consideration the Komar energy [31, 32]. Therefore, we investigate the evolution of the universe using the Komar energy as the energy of the universe enclosed by the Hubble horizon in this study. Evidently, we can also use the Misner-Sharp energy as the energy of the universe and obtain the evolution laws of the universe consistent with astronomical observations. However, this is outside of the scope of the present study.
Conventionally, the cosmic fluid is deemed as an ideal fluid, whose equation of state is
$ \begin{array}{l} p = \omega\rho, \end{array} $
(17) where
$ \omega $ is constant.Here, we would like to argue the validity of Eq. (8) from the viewpoint of gravitational dynamics. Suppose that the gravity can be quantized and has a minimum quantum of area on the order of
$ L_{p}^{2} $ [54], then the horizon with the area A can be divided into$ N = \dfrac{A}{c_{1}L_{p}^{2}} $ cells, where$ c_{1} $ is a numerical factor. Further, if there are$ c_{2} $ microscopic states for the every cell, every cell has$ c_{2} $ degrees of freedom. According to the energy equipartition law$ (6) $ , we obtain the total energy$ \begin{split} E^{'} = \frac{c_{2}}{2}NkT = \frac{c_{2}}{c_{1}}\frac{A}{2L_{p}^{2}}kT. \end{split} $
(18) From Eq. (18), we find that Eq. (9) can be recovered by choosing the relation
$ \frac{c_{2}}{c_{1}} = 1 $ . In contrast, Eq. (18) can be reduced to$ \rho = \frac{3H^{2}}{8\pi L_{p}^{2}} $ , which is the standard Friedmann equation for the flat FRW spacetime, if we choose$ \frac{c_{2}}{c_{1}} = \frac{1}{2} $ and the Misner-Sharp energy of Eq. (12). Hence, the energy equipartition law and Einstein equation are equivalent to some extent when we take the Misner-Sharp energy as the energy of the universe and choose$ \frac{c_{2}}{c_{1}} = \frac{1}{2} $ . Thus, the Misner-Sharp energy seems to be a good choice for the energy of the universe. However, we employ the Komar energy as the energy of the universe in this paper. The reasons are as follows. First, the Komar energy is the true gravitational energy, as shown by Padmanabhan. Second, the accelerated expansion of the present universe cannot be described by the Einstein equation without a cosmological constant, whereas it can be explained by a modified gravitational field equation. In contrast, Eq. (16) can be derived by a modified gravitational field equation. Third, Eq. (8), rather than the Einstein equation, is the ansatz from the perspective of thermodynamics. The above discussions show only the validity of Eq. (8) from the viewpoint of gravitational dynamics and do not imply that only the Minser-Sharp energy can explain the validity of Eq. (8). The reason we use the Misner-Sharp energy to argue the validity of Eq. (8) is that it is convenient to compare with the standard evolution equation derived from the general relativity. Indeed, the modified Friedmann equation can be obtainedby choosing the Komar energy.A question arises naturally: How does the holographic equipartition hold true in non-equilibrium? The holographic equipartition implies that the number of degrees of freedom on the surface equals that in the bulk. With the expansion of the universe, the Hubble horizon is increasing. The number of surface degrees of freedom is likewise increasing, as it is proportional to the area of Hubble horizon. Furthermore, the number of bulk degrees of freedom also rises because of the creation of particles. Therefore, the holographic equipartition is satisfied in the universe, where the particles are being created.
-
In this section, we provide two specific cases as examples to illustrate the validity of the present model. Furthermore, we analyze the specific physical process of the evolution of the universe from the two cases.
In every epoch, we assume that the rate of particle creation has the following formula
$ \begin{split} \Gamma = \gamma\left(\frac{H}{H_{0}}\right)^{n-1}H = \alpha H^{n}, \end{split} $
(19) where
$ \gamma $ is a positive constant of the order one, and n is a non-negative integer,$ H_{0} $ is a quantity with the same dimension as H, and$ \alpha = \frac{\gamma}{H_{0}^{n-1}} $ . The purpose of introducing$ H_{0} $ is to ensure that the dimension of$ \Gamma $ is consistent with that of H. Combining Eq. (4), Eq. (5), and Eq. (17), we obtain the specific modified continuity equation$ \begin{split} \dot{\rho}+3(1+\omega)H\rho\left(1-\frac{\alpha}{3}H^{n-1}\right) = 0. \end{split} $
(20) Further, the dynamic evolution equation of the universe (16) can be reduced to
$ \begin{split} \mid\alpha(1+\omega)H^{n-1}-(1+3\omega)\mid\rho = \frac{3H^{2}}{4\pi L_{p}^{2}}. \end{split} $
(21) Thus Eq. (20) and Eq. (21) constitute the fundamental equations of the evolution of the universe. As long as the two equations are combined, some basic physical quantities of the evolution of the universe, such as the scale factor
$ a(t) $ and energy density$ \rho $ , can be solved. -
A natural choice is
$ n = 2 $ , namely$ \Gamma/H = \alpha H $ . According to astronomical observations and cosmological theories, we know that the Hubble parameter H has been decreasing since the very early stage of the universe except for the period of inflation. In this case, the particle creation rate is particularly high in the early universe and decreases rapidly over time, hence choosing such a particle creation rate to study the evolution of the early universe is reasonable. In this subsection,$ \omega $ is set to$ \frac{1}{3} $ , as the early universe is dominated by radiation.With the above assumption, Eq. (20) and Eq. (21) are reduced to
$ \begin{split} \dot{\rho}+4H\rho-\frac{4\alpha}{3}H^{2}\rho = 0 \end{split} $
(22) and
$ \begin{split} \left|\frac{4\alpha}{3}H-2\right|\rho = \frac{3H^{2}}{4\pi L_{p}^{2}}. \end{split} $
(23) Combining the above equations, we obtain the result
$ \begin{array}{l} (\alpha H-3)(-6H^{3}+4\alpha H^{4}-3H\dot{H}) = 0. \end{array} $
(24) This equation has two solutions
$ \begin{array}{l} H = \alpha/3 \end{array} $
(25) and
$ \begin{split} t = \frac{\alpha}{3}\ln\left|\frac{4\alpha H-6}{H}\right|+\frac{1}{2H}, \end{split} $
(26) where the integration constant is chosen as zero.
The solution in Eq. (25) implies that there is an inflation solution
$ a(t)\propto {\rm e}^{\frac{\alpha t}{3}} $ in the early universe. In contrast, the evolution law of the universe is$ H = \frac{1}{2t} $ when t is large from Eq. (26), which can be reduced to the standard evolution law$ a(t)\propto t^{1/2} $ in the radiation dominated stage. These solutions can be explained as follows. The universe starts from an unstable de Sitter space$ (a(t)\propto {\rm e}^{\frac{\alpha t}{3}} $ ), after which it evolves to the standard radiation phase$ (a(t)\propto t^{1/2}) $ . During the evolution of the early universe, the creation rate of the radiation is consistently$ \Gamma = \alpha H^{2} $ . -
Another natural choice is
$ n = 1 $ , namely$ \Gamma/H = \alpha $ . With this choice, Eq. (20) and Eq. (21) are reduced to$ \begin{array}{l} \dot{\rho}+(1+\omega)\left(3-\alpha\right)H\rho = 0 \end{array} $
(27) and
$ \begin{split} \mid(\alpha-1)+(\alpha-3)\omega\mid\rho = \frac{3H^{2}}{4\pi L_{p}^{2}}. \end{split} $
(28) From Eq. (27) and Eq. (28), we obtain the results
$ \rho = {\rm constant} $ and$ H^{2} = 8\pi L_{p}^{2}\rho/3 $ if the parameter$ \omega $ is set to$ -1 $ . The results are consistent with the evolution laws of the universe when the universe is dominated by vacuum energy in general relativity. In contrast, the same results are obtained if$ \alpha = 3 $ . This implies that the evolution law of the universe is$ a(t)\propto {\rm e}^{Ht} $ and independent of the nature of matter when the particle creation rate$ \Gamma $ equals$ 3H $ . Therefore, we can describe the exponentially expanding de Sitter universe without introducing vacuum energy, which is equivalent to the negative pressure matter when$ \alpha = 3 $ .When
$ \omega\neq-1 $ and$ \alpha\neq3 $ , by combining Eq. (27) and Eq. (28), we obtain$ \begin{split} \dot{H} = -\frac{1}{\delta} H^{2}. \end{split} $
(29) Further, we obtain
$ \begin{array}{l} a(t) = t^{\delta}, \end{array} $
(30) where
$ \delta = \frac{2}{(1+\omega)(3-\alpha)} $ . The parameter$ \delta $ satisfies the inequality$ \delta>1 $ , as the expansion of the universe is speeding up according to astronomical observations. If the universe is dominated by pressureless matter at present, namely$ \omega = 0 $ , we obtain the range of the parameter$ \alpha $ as$ 1<\alpha<3 $ . From this result, we can describe the accelerated expansion of the universe without introducing dark energy. For example, the authors of Ref. [59] pointed out that the rate of expansion, which is consistent with supernova observations, is currently$ a(t) = t^{2} $ . This rate of expansion can be obtained when$ \alpha = 2 $ for the present universe, dominated by pressureless matter. Hence, we can describe the accelerated expansion of the current universe and the evolution of the de Sitter universe at a later time using these settings.Thus far, we analyzed the specific physical processes of the evolution of the universe from the cases of
$ \Gamma/H = \alpha H $ and$ \Gamma/H = \alpha $ based on particle creation and holographic equipartition. These results are consistent with the conclusions of the analysis of$ \Lambda $ dark energy model, which is generally accepted. Therefore, the current model can efficiently explain the entire evolution of the universe without introducing dark energy.At the end of this section, we discuss the evolution mechanism of
$ \Gamma $ . We observe that$ H^{2} $ evolves to H naturally if the creation rate of particles is chosen as$ \Gamma = \gamma(\frac{H^{2}}{H_{0}}+H) $ in this section. At the early universe, the term$ \gamma\frac{H^{2}}{H_{0}} $ dominates. The Hubble parameter H decreases as the universe evolves, the second term$ \gamma H $ begins to dominate. Thus, the behaviour of$ \Gamma $ from$ H^{2} $ to H changes naturally at the matter-radiation equality. -
We investigate the thermodynamical properties of the universe based on the holographic equipartition in this model. To observe the energy equipartition, we calculate the physical quantity
$ \begin{split} \frac{1}{2}\beta E = \frac{1}{2}\frac{2\pi}{kH}\frac{1}{HL_{p}^{2}} = \frac{\pi}{H^{2}L_{p}^{2}}. \end{split} $
(31) Hence, the relation
$ \begin{split} S = \frac{1}{2}\beta E \end{split} $
(32) holds in this model, as
$ S = \frac{A}{4L_{p}^{2}} = \frac{\pi}{H^{2}L_{p}^{2}} $ . This is an important relation of thermodynamics that has been shown by Padmanabhan in the cases of static spacetime in general relativity [27] and a wider class of gravity theories, like the Lanczos-Lovelock gravity [60]. This relation is indeed the energy equipartition law in the bulk when the holographic equipartition$ N_{\rm sur} = N_{\rm bulk} $ is satisfied, as$ E = \frac{2}{\beta}S = \frac{kT}{2}\frac{A}{L_{p}^{2}} = \frac{1}{2}N_{\rm sur}kT = \frac{1}{2}N_{\rm bulk}kT $ . Hence, the ansatz that the holographic equipartition remains satisfied throughout the evolution of the universe is consistent with the conclusion that the relation$ S = \frac{1}{2}\beta E $ depicts the energy equipartition in static spacetime. Moreover, the application of the energy equipartition law is generalized to dynamic spacetime. Consequently, the validity of relation presented in Eq. (32) also indicates that it is a reasonable ansatz that the universe always obeys the holographic equipartition.Furthermore, differentiating Eq. (10), we obtain the change of total energy enclosed by the Hubble horizon during a time interval
$ {\rm d}t $ $ \begin{split} {\rm d}E = -\frac{\dot{H}}{H^{2}L_{p}^{2}}{\rm d}t. \end{split} $
(33) On the other hand, we obtain
$ \begin{split} T{\rm d}S = \frac{H}{2\pi}\left(-\frac{2\pi\dot{H}}{H^{3}L_{p}^{2}}\right){\rm d}t = -\frac{\dot{H}}{H^{2}L_{p}^{2}}{\rm d}t, \end{split} $
(34) where the temperature
$ T = \frac{H}{2\pi} $ and the entropy-area relation$ S = \frac{A}{4L_{p}^{2}} $ are used. Comparing Eq. (33) with Eq. (34), we obtain the relation$ \begin{array}{l} {\rm d}E = T{\rm d}S. \end{array} $
(35) This relation depicts the energy conservation relation for the universe enclosed by the Hubble horizon. However, the first law of thermodynamics is expressed as
$ -{\rm d}E = T{\rm d}S $ in Ref. [61] (the relation is also used in some references, e.g., Refs. [44, 62-66]). We explain the difference between the notion$ {\rm d}E $ used in Ref. [61] and that used in the present paper. In Ref. [61],$ -{\rm d}E $ is the heat flux crossing the horizon during time$ {\rm d}t $ . However, in this case$ {\rm d}E $ is the change of the energy inside the Hubble horizon. Moreover,$ -{\rm d}E $ is defined by$ -{\rm d}E = 4\pi R^{2}T_{\mu\nu}k^{\mu}k^{\nu}{\rm d}t $ , where$ k^{\mu} $ is the future directed ingoing null vector field in Ref. [61]. As$ {\rm d}E $ is mainly composed of two parts, one related to the particle creation and the other to the particles across the Hubble horizon, the change of entropy in this study is mainly caused by the change in the number of particles. Further, the energy E in the bulk can be also considered as the energy on the horizon due to the holographic equipartition. Further, the entropy change on the horizon$ {\rm d}S $ is the irreversible entropy change attributed to particle creation. Therefore, non-equilibrium thermodynamics accompanying the entropy production can be expressed in the form of Eq.$ (35) $ of the equilibrium thermodynamics on the Hubble horizon.