Emergent Universe in the Braneworld Scenario

According to Padmanabhan's proposal, the difference between the surface degrees of freedom and the bulk degrees of freedom in a region of space may result in the acceleration of Universe expansion through the relation $\Delta V/\Delta t = N_{\rm sur}-N_{\rm bulk}$ where $N_{\rm bulk}$ and $N_{\rm sur}$ are referred to the degrees of freedom related to the matter and energy content inside the bulk and surface area, respectively \cite{Pad1}. In this paper, we study the dynamical effect of the extrinsic geometrical embedding of an arbitrary four dimensional brane in a higher dimensional bulk space and investigate the corresponding degrees of freedom. Considering the modification of Friedmann equations arising from a general braneworld scenario, we obtain a correction term in Padmanabhan's relation, denoting the number of degrees of freedom related to the extrinsic geometry of the brane embedded in higher dimensional spacetime as $\Delta V /\Delta t=N_{\rm sur}-N_{\rm bulk}-N_{\rm extr}$ where $N_{\rm extr}$ is referred to the degree of freedom related to the extrinsic geometry of the brane while $N_{\rm sur}$ and $N_{\rm bulk}$ are as well as before. Finally, we study the validity of the first and second laws of thermodynamics for this general braneworld scenario in the state of thermal equilibrium and in the presence of confined matter fields to the brane with the induced geometric matter fields.


Introduction
Recent researches support the idea that the gravitational field equations can be derived in the same way that the equations of an emergent phenomena like fluid mechanics or elasticity are obtained [1,2]. According to emergent gravity paradigm the gravitational field equations can be derived from thermodynamic principle [3,4]. In this way, Padmanabhan has treated the Einstein field equations as emergent, while the existence of a spacetime manifold, its metric and curvature have been assumed [5]. In a cosmological context, it has been argued that the accelerated expansion of the Universe can be derived from the difference between the surface and bulk degrees of freedom through the relation ∆V /∆t = N sur −N bulk , in which N bulk and N sur are referred to the degrees of freedom related to matter and energy content (or dark matter and dark energy) inside the bulk and surface area, respectively [6]. In order to explain the present accelerated expansion of the Universe, which is in agreement with different data set or observational data [8], different models have been proposed. One of these models is the dark energy model which admits that the universe is dominated by a dark fluid with negative pressure. However, there are several dark energy models such as dynamical dark energy [9], quintessence [10] and k-essence [11], for a review the reader is referred to [12]. Also, the LCDM model or the concordance model is a particular case of dark energy that is parameterized by a cosmological constant Λ with the equation of state parameter equal minus one, i.e p = −ρ. The strong energy condition, i.e ρ+3p > 0, is violated by the dark energy because of demanding for the positive acceleration of the Universe through the second Friedmann equationä a = − 4πG 3 (ρ + 3p). Another approach lies in the framework of modified gravity theories which describe the present acceleration of the Universe such as f (R) gravity [13], f (T ) gravity [14], Hořava-Lifshitz gravity [15], Gauss-Bonnet gravity [16], Weyl gravity [17], Lovelock gravity [18], massive gravity [19] and braneworld scenarios [20]. In these models, there is no need for the introduction of an adhoc component usually called as dark energy with unusual features. In these models, some additional terms are considered in the gravitational Lagrangian which lead to the modification of the gravitational theory resulting in an effective dark energy sector with a geometrical origin. One can also find such models in the low energy limit of heterotic string theory [21]. All of these models admit a series of conditions coming from various laws of physics such as thermodynamics laws [22] or astrophysical data.
To explain the structure of spacetime and its relation with thermodynamics of the system, one can refer to four laws of black hole mechanics which are derived from the classical Einstein field equations. These four laws are analogous to those of thermodynamics [23]. Discovery of the quantum Hawking radiation [24] turns out that this analogy is an identity. By deriving the Einstein field equation from the relation of entropy and horizon area together with the thermodynamic law of Q = T dS which connects the heat Q, the entropy S, and the temperature T , Jacobson showed that the classical general relativity is a kind of thermodynamics where the surface gravity is a temperature [25]. The generalized second law of thermodynamics is specially investigated in different modified gravity models. For example, we can refer to the investigations devoted to the study of generalized second law (GSL) of thermodynamics in f (T ) gravity models in which two types of horizons, are used to check the validity of the generalized second law of thermodynamics with corrected entropies [26]. One can also find that in the state of thermal equilibrium, in Kaluza-Klein universe which is composed of dark matter and dark energy, the validity of the laws of thermodynamics are true [27]. The investigations on the deep connection between gravity and thermodynamics have been widely considered in the cosmological context where it has been shown that in the form of the first law of thermodynamics on apparent horizon, the differential form of the Friedmann equation in the FRW universe can be written [28,29,30,31,32,33,34,35,36]. The GSL in an accelerating universe related to the apparent horizon has been considered in [37,38,39]. It was discussed in [38], that in contrast to the case of the apparent horizon, the general second law of thermodynamics breaks down in the case of a universe enveloped by the event horizon with the usual definitions of entropy and temperature. This study reveals that from the thermodynamical point of view, in an accelerating universe with spatial curvature, the apparent horizon is a physical boundary. Also, the general expression of temperature at apparent horizon of FRW universe, allows one to show that the GSL holds in Einstein, Gauss-Bonnet and more general Lovelock gravity [40]. Also, the GSL of thermodynamics in the framework of braneworld scenarios is studied in [41]. One can find other studies on the GSL of thermodynamics in [42,43,44].
In particle physics, the warped product geometries, well known as Randall-Sundrum models, are very important [45]. In these models, it is imagined that our real world is a higher-dimensional universe described by a warped geometry with Z 2 symmetry. The standard gauge interactions are confined to the four dimensional brane embedded in a higher dimensional bulk space where gravitons are propagating through the extra dimensions. More specifically, our universe is assumed to be a five-dimensional anti-de Sitter space where the elementary particles, except for the gravitons, are localized on a (3 + 1)-dimensional brane or branes.
In this paper, we consider a general braneworld model which provides a geometrical origin for dark energy or accelerating expansion of the Universe [46]. Considering the modification of Friedmann equations resulted from this general braneworld scenario, we obtain a correction term on Padmanabhan's relation. This paper is organized as follows: In section 2, we introduced general geometrical setup of the braneworld. In section 3, this braneworld model is studied under Israeal-Darmois-Lanczos junction condition, which provides the Z 2 symmetry, and the corresponding number of degrees of freedom related to the extrinsic geometry of such a brane model is obtained. In section 4, we find the correction term to the Padmanabhan's relation in our general braneworld model which does not have any specific junction condition. In section 5, we explored the thermodynamics of such a general brane model. At last, in section 6, we presented our concluding remarks.

General Geometrical Setup of the Braneworld
The effective Einstein-Hilbert action for the 4D spacetime (M 4 , g) embedded in an nD bulk space (M n , G) can be written where α * , R and L m are respectively as gravitational coupling constant in the bulk space, bulk Ricci scalar and the Lagrangian density of the matter fields confined to the brane. Variation of the action (1) with respect to the bulk metric G AB (A, B = 0, ..., n − 1) leads to the following field equations for the bulk space where S AB is and T AB = −2 δLm δgAB + g AB L m is the energy-momentum tensor of the matter fields confined to the four dimensional brane through the action of the confining potential V. The confining potential V satisfies three general conditions: (I) It has a deep minimum on the original non-perturbed brane (we will discuss on the original non-perturbed and perturbed geometry in the following), (II) It depends only on extra coordinates, and (III) It preserves the gauge group related to the subgroup of the isometry group of the bulk space [48]. Using the confining potential V the matter fields are exactly localized on the brane and one obtains where µ, ν = 0, ..., 3 and a, b = 4, ..., n − 1 labels the number of four dimensional brane and bulk extra dimensions, respectively and T µν is the confined matter source on the four dimensional brane. This is the so called "confinement hypothesis". Now, it is worth to have a brief discussion on the bulk and brane energy scales and their corresponding gravitational coupling constants. The bulk space gravitational coupling constant α * is where G n is equivalently known as the bulk gravitational constant and M n is the fundamental energy scale of the bulk space. In a usual four dimensional spacetime, we have G 4 = G = M −2 P l where G is the Newtonian gravitational constant. In the static weak field limit of the Einstein field equations, one obtains n-dimensional Poisson equation for the gravitational potential which admits the following solution where by supposing that the length scale of the extra dimensions is denoted by L, this potential behaves as V (r) ∼ r 3−n , on scales with size r L, and depends on the number of dimensions of spacetime n. On the other hand, for the scales larger than L, the potential V (r) behaves as V (r) ∼ L 4−n r −1 [49]. For n = 4, we recover the Newtonian four dimensional gravitational potential V (r) ∼ r −1 . This means that the Newtonian gravitational constant G or usual Planck scale M P l are effective coupling constants, which describe gravity on the scales much larger than the length scale of extra dimensions, and are proportional to the bulk fundamental energy scale M n via where L n−4 denotes the volume of the extra dimensional space. This relation indicates that for the extra dimensional volume which is about the Planck scale, i.e. L ∼ M −1 P l , we have M n ∼ M P l . But for the volume which is significantly above the Planck scale, we find that the fundamental energy scale of the bulk space M n is much smaller than the four dimensional effective energy scale M P l ∼ 10 19 GeV.
In order to obtain the effective Einstein field equation induced on the brane, we consider the following geometrical setup. Suppose that the 4D background manifold M 4 is isometrically embedded in a n dimensional bulk M n by a differential map Y A : M 4 −→ M n such that where is the basis of the bulk (brane),N A a are (n − 4) normal unit vectors orthogonal to the brane and g ab = ǫδ ab in which ǫ = ±1 corresponds to two possible signatures for each extra dimension of the bulk space. Perturbation of the background M 4 manifold in a sufficiently small neighborhood of the brane along an arbitrary transverse direction ξ a is given by where L ξ represents the Lie derivative along ξ a denoting the non-compact extra dimensions. The presence of tangent component of the vector ξ along the brane can cause some difficulties because it can induce some undesirable coordinate gauges. But, it was shown that in the theory of geometric perturbations, it is quite possible to choose this vector to be orthogonal to the background manifolds [50]. Then, choosing the extra dimensions ξ a to be orthogonal to the brane ensures us about the gauge independency [51] and having perturbations of the geometrical embedding along the orthogonal extra directions N A a . Thus, the local coordinates of the perturbed brane will be Equation (9) implies that since the vectors N A depend only on the local coordinates x µ , N A = N A (x µ ), they will not propagate along the extra dimensions which can be shown as The above assumptions lead to the embedding equations of the perturbed geometry as where by setting N A a = δ A a , the metric of the bulk space G AB in the Gaussian frame and in the vicinity of M 4 takes the form of Then, the line element of the bulk space will have the following form where is the metric of the perturbed brane whileḡ µν is the metric of original non-perturbed brane (the first fundamental form) andK is the extrinsic curvature of the original brane (the second fundamental form). In what follows, we will use the notation which represents the twisting vector fields (the normal fundamental form). For any fixed extra dimension ξ a , we have a new perturbed brane and can define an extrinsic curvature similar to the original one bỹ Note that the definitions (13), (15) and (18) require the extrinsic curvature of the perturbed brane to bẽ In a geometric setup, the presence of gauge fields A µa tilts the embedded family of sub-manifolds with respect to the normal vector N A . According to our construction, although the original brane is orthogonal to the normal vector N A , equations (12) imply that this will not true for the deformed geometry. Hence, we change the embedding coordinates in the following form where the coordinates X A describes new family of embedded manifolds whose members always will be orthogonal to the normal vector N A . In this coordinate system, the embedding equations of the perturbed brane will be similar to the original one, described by the relations given in equation (8), so that the coordinates Y A are replaced by the new coordinates X A . This geometrical embedding of the new local coordinates will be suitable for obtaining the induced Einstein field equations on the brane. In this coordinates, the extrinsic curvature of a perturbed brane is given by which is known as the generalized York's relation and shows the propagation of the extrinsic curvature because of the propagation of the metric in the direction of extra dimensions in the bulk space. The Gauss-Codazzi equations for the components of the Riemann tensor of the bulk space in the embedding vielbein where R ABCD and R αβγδ are the Riemann tensors of the bulk and the perturbed brane, respectively [52]. Then, one can find the Ricci tensor by contracting the Gauss equation (22) as where a next contraction will give the Ricci scalar as where use has been made of the notations K • K ≡ K aµν K aµν and K a ≡ g µν K aµν . Consequently, using equations (24) and (25), the relation between Einstein tensors of the bulk and brane can be obtained as where G AB , G µν are the Einstein tensors of the bulk and brane respectively, and where K a = g µν K aµν and K • K = K aµν K aµν . From the definition of Q µν , it is an independent conserved geometrical quantity, i.e. ∇ µ Q µν = 0 [46]. Using the decomposition of the Riemann tensor of the bulk space into the Weyl curvature tensor, the Ricci tensor and the scalar curvature as we obtain the induced 4D Einstein equation on the brane as where E µν = g ab C ABCD X A ,µ N B a N C b X D ,ν is the electric part of the Weyl tensor C ABCD of the bulk space . From the brane point of view, the electric part of the Weyl tensor describes a traceless matter, denoted by dark radiation or Weyl matter. For a constant curvature bulk space, we have E µν = 0.
Then, the induced Einstein equation, in a constant curvature and Ricci flat bulk (E µν = R AB = 0) will take the following form where T µν is the confined matter source on the brane and Q µν is a pure geometrical energy-momentum source. We also assume that the spacetime on the brane is isotropic and homogeneous and so we have Friedmann-Robertson-Walker (FRW) metric on the brane, where a(t) is the cosmic scale factor, k = +1, −1 and 0 corresponds to the closed, open and flat universes, respectively, and dΩ 2 = dθ 2 + sin 2 θdφ 2 . The confined matter source on the brane T µν can be considered in the perfect fluid form in a co-moving coordinates as where u α = δ 0 α is the 4-velocity vector of the fluid, ρ and p are energy density and isotropic pressure, respectively. For the metric (31), the components of the extrinsic curvature tensor can be obtained by using the Codazzi equation as where dot denotes the derivative with respect to the cosmic time t, and b = b(t) is an arbitrary function of time [46,53]. Then, by defining the parameters h(t) =ḃ/b and H(t) =ȧ/a, the components of Q µν represented by (27) take the form of Similar to the confined matter field source on the brane T µν , the geometric energy-momentum tensor Q µν can be identified as Q µν = (ρ extr + p extr )u µ u ν + p extr g µν , where the ρ extr and p extr denote the "extrinsic geometric energy density" and "extrinsic geometric pressure", respectively (the suffix "extr" stands for "extrinsic") [46]. Then, using Eqs. (34) and (35) we obtain Using Eqs. (30), (34) and (36) and separating the space and time components we arrive aẗ Eliminating theä terms gives the following modified Friedmann equation on the brane which possesses a modified term arising from the extrinsic geometry of the brane in the bulk space. In the next sections, we will study this braneworld modification to the Friedmann equations with more details.

The Brane Model with Junction Conditions
Using the Israeal-Darmois-Lanczos junction condition which exactly provides the Z 2 symmetry 1 (mirror symmetry) [47] (see [46] for a brief review), one can obtain the extrinsic curvature tensor component of the original non-perturbed brane in terms of the confined matter sources on brane as k 11 = b(t) = −α 2 * ρa 2 [46]. Then, the modified Friedmann equations (37), (38) and (39) will take the forms which shows the ρ 2 dependent cosmology [54]. Now, we intend to obtain modification of the basic law governing the emergence of space due to the difference between the degrees of freedom in the framework of this model. Using relationä/a =Ḣ + H 2 , Eq. (41), can be written asḢ Multiplying Eq. (42) Assuming V = 4πH −3 /3 as the volume of the sphere on the brane with Hubble radius H −1 , we have On the other hand, according to Padmanabhan's idea, the number of degrees of freedom on the spherical surface of Hubble radius H −1 is given by [6] where L p is the Planck length and A = 4πH −2 represents the area of the Hubble horizon. Using the area law S = A/4L 2 p , as the saturation of Bekenstein limit [7], we can write 2 N sur = 4S.
1 The Z 2 symmetry means that when you approach the brane from one side and go through it in the bulk, you face with the same bulk having reversed normal unit vector to the brane, i.e N a → −N a . Indeed, in the presence of Z 2 symmetry, the original non-perturbed brane located at ξ a = 0 acts as a mirror for all objects that feels the extra dimensions. The Z 2 symmetry governs for any perturbation of the original brane leading to a mirror perturbation on the other side of the brane [50]. 2 The Bekenstein limit is an upper limit on the entropy or information that can be contained within a given finite region of space which has a finite amount of energy. It implies that the information necessary to perfectly describe a system, must be finite if the region of space and the energy is finite.

Also, the bulk degrees of freedom obey the equipartition law of energy
where E, k B and T are the energy inside of the bulk, the Boltzmann constant and the temperature of the bulk, respectively. In the following, we use the units of k B = c = = G = L p = 1 for simplicity. We also assume the temperature associated with the Hubble horizon as the Hawking temperature T = H/2π, and the energy contained inside the Hubble volume in Planck units V = 4π/3H 3 as the Komar energy The novel idea of Padmanabhan is that the cosmic expansion, conceptually equivalent to the emergence of space, is being driven towards holographic equipartition, and the basic law governing the emergence of space must relate the emergence of space to the difference between the number of degrees of freedom in the holographic surface and the ones in the emerged bulk [6]. Using equations (47) and (48), the bulk degrees of freedom may be obtained as where it is assumed that ρ + 3p < 0. Thus, Eq. (44) can be written as where N extr = 3 2πǫ appears as the number of degrees of freedom corresponding to the extrinsic geometry of the embedded brane in a higher dimensional spacetime. Indeed, there are three modes of degrees of freedom, the surface degrees of freedom, the bulk degrees of freedom and the ones that are related to the extrinsic geometry of the embedded brane. Since N extr represents the number of degrees of freedom, it must be positive. Equation (51) shows that the positiveness of N extr demands for ǫ = +1, representing a spacelike extra dimension. Therefore, it turns out that by applying the Israeal-Darmois-Lanczos junction condition, the timelike extra dimension will be ruled out. This indicates that unlike the other braneworld scenarios where there is no essential requirement for ǫ being positive, in the present scenario the positiveness of ǫ is imposed by the equation (51). Generally, the braneworlds with different extrinsic geometries have different cosmological evolutions. It is seen that in this scheme, by applying the Israeal-Darmois-Lanczos junction condition, the number of degrees of freedom depends on the bulk space energy scale α * , the signature of the extra dimensions ǫ = +1 and the confined matter density ρ as well as the volume V and horizon temperature T . At the end of this section, we remark that the presence of the quadratic energy density in the Friedmann equations which was initially anticipated as a possible solution to the observed accelerated expansion of the Universe, was shown to be incompatible with the big-bang nucleosynthesis [55]. Also, it is shown that this quadratic ρ term can constrain the high energy inflationary regimes in comparison with the observational SDSS/2DF/WMAP data [56]. In order to reconcile the above mentioned braneworld scenario with Z 2 symmetry or Israeal-Darmois-Lanczos condition with the observational data, one may propose that this scenario should be modified, see [57]. Therefore, in the following section we will study the case of a general braneworld embedding procedure without any simplifying junction condition or Z 2 symmetry.

The General Braneworld Model without any Specific Junction Condition
We consider the geometric quantity (35) with the barotropic equation of state where ω extr is the geometric equation of state parameter and generally can be a function of time. Using equations (36) and (52), we obtain the following equation for b(t) where ω extr is an unknown function. In general, solving the above equation is impossible unless the functional form of ω extr is given. Let us consider the simple case where ω extr = constant. In this case Eq. (53) can be solved immediately as where a 0 = a(t 0 ) is the scale factor of the Universe at the present time and b 0 is an integration constant representing the curvature warp of the Universe at the present time. Substituting the solution (54) into equations (34) gives the components of the geometric quantity in terms of b 0 , a 0 and a(t) as and consequently using equations (36) we get Then, using equations (56) and (35), the induced Einstein equation on the brane (30) gives us the following equation for the confined energy density Note that we have not included the cosmological constant because it is possible to construct a geometrical origin for the dark energy in a general geometrical embedding scheme with a brane possessing an extrinsic curvature, to recover the acceleration of the Universe [46,58]. The generalization to the case that the cosmological constant is not zero, is trivial. Similarly, the confined isotropic pressure component can be obtained from equations (30), (35) and (56) as Combining these equations leads to the following equation Using the same procedure as in the previous section, we obtain where the number of degrees of freedom related to the extrinsic geometry of spacetime has the general form For the case of a general geometric embedding scheme, the number of degrees of freedom related to the geometric embedding state of the brane in a higher dimensional bulk beside the scale factor of the Universe a, depend on the signature of extra dimensions ǫ, the volume V and horizon temperature T and warp factor of the Universe b 0 as well as the equation of state parameter of the geometric fluid ω extr . If the curvature warp of the universe b 0 vanishes, all of the extrinsic curvature components will also vanish, and the braneworld will behave just as a trivial plane. In this case, the number of degrees of freedom corresponding to the extrinsic curvature vanishes and we recover the original Padmanabhan's relation dV dt = N sur − N bulk . Another interesting result is that when the geometric equation of state parameter becomes ω extr = −1/3, the corresponding degrees of freedom also vanishes. Moreover, there are two possibilities for satisfying the positivity of N extr . The first possibility is ω extr > −1/3 with ǫ = +1 which indicates a spacelike extra dimension. For this case, all of the known energy conditions such as weak, null, strong and dominant energy conditions are satisfied for the geometric fluid. The second possibility is ω extr < −1/3 with ǫ = −1, which accounts for a timelike extra dimension. In this case, the energy conditions may be violated by the geometric fluid. It is worth mentioning that, according to (56), satisfying the weak energy condition requires ǫ = +1 which is the same result coming from the positiveness of N extr .
Our Universe, is not pure de Sitter but we know that it evolves toward an asymptotically de Sitter phase. For the purpose of reaching holographic equipartition we need to have dV /dt → 0 in the equation (60) which leads to N sur = N bulk + N extr . In order to understand the prominent feature of N extr it is better to look at equation (60) without this term. Following the discussion of Padmanabhan, one can consider that N bulk consists of two terms, one related to the normal matter with ρ + 3p > 0, and the other one related to the dark energy with ρ + 3p < 0 [6]. So, it is possible to divide the degrees of freedom of bulk into two terms, one coming from the degrees of freedom of dark energy leading to acceleration and the other one coming from the degrees of freedom of normal matter leading to deceleration. Then, equation (60) takes the form of dV dt = N sur + N m − N de . Thus, it is seen that a universe without a dark energy component has no hope of reaching the holographic equipartition [5].
We can remark the important results of the present study as follows.
• We can avoid of the term N de , namely dark energy or cosmological constant, which has been proposed by Padmanabhan. In our general setup, dark energy has a completely geometrical origin [46]. In fact, the geometrical component denoted by N extr plays the role of N de proposed by Padmanabhan. Similarly, we can understand equation (60) in a better way if we separate out the matter component resulting in deceleration from the geometrical component resulting in acceleration. For the sake of simplicity, we will assume that the Universe has just two components, the normal matter with ρ + 3p > 0 and an effective (geometric) matter with ρ extr + 3p extr < 0 . By our consideration, equation (60) can be expressed in an equivalent form as where N sur , N m , N extr are positive with It is seen that the holographic equipartition condition with asymptotically vanishing emergence of space (dV /dt → 0) can be satisfied only if the Universe possesses N extr . Equivalently, the existence of a geometric term (due to the embedding of the brane) is required for the asymptotic holographic equipartition which leads the cosmos to find its equilibrium. In the presence of N extr , the emergence of space will lead to N extr dominating over N m with the universe experiencing accelerated expansion due to dark energy. Asymptotically, N extr will approach N sur and dV /dt → 0 in a de Sitter universe.
• We can keep the term N de as the dark energy degrees of freedom in the bulk and consider the new term N extr as a geometric contribution to the term N de . In this regard, we can write the following equation The holographic equipartition condition is satisfied if N sur + N m = (N de + N extr ). Moreover, we have In the presence of (N de + N extr ), the emergence of space will lead to (N de + N extr ) dominating over N m with the universe experiencing accelerated expansion due to dark energy and geometric embedding. Asymptotically, (N de + N extr ) will approach N sur and dV /dt → 0 in a de Sitter universe.

Thermodynamics of a General Braneworld Scenario
The first law of thermodynamics for apparent horizon 3 reads as [61] − dE = T dS, where T is the time-dependent temperature of a thermal heat bath as is perceive by an observer at r = 0, dE is the change in the mass of the matter present on the observer's side of the horizon, and dS is the increases in the horizon entropy. In this section, we confirm the validity of the first law of thermodynamics for apparent horizon in the presence of the additional terms due to the extrinsic curvature of a braneworld model and then we study the second law of thermodynamics with the assumption that in the apparent horizon the space-time has thermodynamical behaviour. We also show that the second law of thermodynamics for apparent horizon is always satisfied for an expanding universe. The entropy of system is obtained as the sum of the surface entropy and the internal entropy. The internal entropy includes the entropy related to the ordinary matter fields localized on the brane and the geometric entropy corresponding to the induced geometric matter. In this section, we consider a general embedding scheme without any specific junction condition.
For a flat universe k = 0, we consider a perfect fluid form for the geometric fluid Q µν as in equation (35) and using equations (56), (57), (59) and conservation equation we obtain three following equationṡ ρ +ρ extr + 3H(ρ + p + ρ extr + p extr ) = 0, As is seen from the above equations, there are two kind of matter sources for these equations. The first one is the normal matter ρ confined on the brane and the second one is the induced geometric matter ρ extr . In order to obtain the entropy expression associated with the normal and geometric matter, we consider the variation of their corresponding energies, dE Am and dE Aextr , which are achievable by the energy crossing formula on the apparent horizon as [61] − dE Am = 4πR 2 T µν K µ K ν dt = 4πR 2 (ρ + p)dt, where we have used the subscript A to denote for the quantities on the "Apparent horizon". In order to achieve the total energy crossing formula, we should consider equations (71) and (72) together with dE A = dE Am + dE Aextr . Then, we obtain − dE A = 4πR 2 (T µν + Q µν )K µ K ν dt = 4πR 2 (ρ + p + ρ extr + p extr )dt.
Then, using the equations (74) and (75), we obtain Inserting the equation (76) into the equation (73) and considering H −1 instead of R A , the energy crossing term will take the form of We note that the radius of apparent horizon is R A = (H 2 + k a 2 ) −1 2 where for the flat spatial geometry k = 0, the radius of the apparent horizon R A coincides with hubble horizon R A = 1 H . Thus, by taking into account the area law of the entropy, i.e S A = A/4 = 4πH −2 /4, we find which denotes the total surface crossing entropy of the system. Thus, one can confirm the validity of the first law of thermodynamics using equations (77) and (78) as where T A = (2πR A ) −1 . Using (78), we obtainṠ On the other hand, we have the internal entropy S I for the system which is related to the volume inside the horizon (we have used the subscript I to denotes for the quantities "Inside apparent horizon"). For the internal entropy, we have where P = p + p extr and E I = (ρ + ρ extr )V. Variation of (81) results iṅ where we have divided it to the entropy corresponding to the normal matter S m and the geometric matter S extr . The extrinsic geometric entropy shows its effect through the induced geometric fluid on the brane by adding a new term to the total internal entropy. In the above equation, T I is the temperature of the thermal system inside the horizon. We consider a thermal system which is bounded by an apparent horizon that has reached equilibrium with its internal volume. This assumption allows us to put T I = T A [62]. By putting V = 4π 3 H −3 and using equations (75) and (76) in the equation (82), we obtaiṅ Then, the total derivative of the entropy by adding equation (5) to equation (80), becomeṡ According to the second law of thermodynamics, entropy of the thermodynamical systems can never decrease. So, the derivative of the entropy with respect to time is always greater than zero, i.e.,Ṡ t 0, so we havė For the Universe which expands, H > 0, the above equation is always true representing that the second law of thermodynamics always holds. Then, the entropy of the Universe always increases and it depends on the normal matter ρ and the induced geometric matter ρ extr , and by looking at equation (70), we see that it is independent of the pressure profiles p and p extr . Therefore, in our braneworld model, similar to the Kaluza-Klein model withṠ t = 21π 2 8G5H 2 0, the entropy of the Universe on the apparent horizon always increases [27].

Conclusion
In this paper, we have addressed the following question: what is the dynamical effect of the extrinsic geometrical embedding of an arbitrary four dimensional brane in a higher dimensional bulk space, from Padmanabhan's point of view on the emergent Universe? We have shown that other than the surface degrees of freedom and the bulk degrees of freedom, there are new degrees of freedom related to the extrinsic geometry of the brane embedded in a higher dimensional bulk space which may play a basic role in the cosmological evolution. Based on this scenario, we have corrected the Padmanabhan's relation as ∆V /∆t = N sur − N bulk − N extr where N extr accounts for the new degrees of freedom corresponding to the extrinsic geometry of the brane. This term has a contribution to the Padmanabhan's relation such that it plays the role of a dark energy. In this regard, we have separately investigated the braneworld scenarios with and without specific junction conditions. Moreover, we have shown that for the case of braneworlds with Israeal-Darmois-Lanczos junction condition, the number of degrees of freedom is determined by the bulk space energy scale α * , the signature of the extra dimensions ǫ = +1 and the confined matter density ρ as well as the volume V and horizon temperature T . In this case, because of positivity of the number of degrees of freedom, N extr , the possibility of having timelike extra dimension is ruled out.
For the case of a general geometric embedding scheme, the number of degrees of freedom related to the geometric embedding state of the brane in a higher dimensional bulk space, beside the scale factor of the Universe a, depends on the signature of extra dimensions ǫ = ±1, the volume V , the horizon temperature T and the warp factor of Universe b 0 , as well as the equation of state parameter of the geometric fluid ω extr . If the curvature warp of the Universe b 0 vanishes, all of the extrinsic curvature components vanish too, and the braneworld behaves just as a trivial plane. In this case, the corresponding number of degrees of freedom also vanishes and we recover the original Padmanabhan's relation dV dt = N sur − N bulk . Also, when the geometric equation of state parameter will be ω extr = −1/3, the corresponding degrees of freedom also vanishes in this approach. Moreover, the positivity of of the number of degrees of freedom requires ω extr > −1/3 for a spacelike extra dimension where ω extr < −1/3 for a timelike extra dimension. For the first case, all of the known energy conditions as weak, null, strong and dominant energy conditions are satisfied by the geometric fluid while they may be violated by the geometric fluid for the second case. Then, we investigated the thermodynamical aspects of this general braneworld model. We confirmed the validity of the first law of thermodynamics on apparent horizon in the presence of the additional terms due to the extrinsic curvature of a braneworld model and then we studied the second law of thermodynamics with the assumption that on the apparent horizon the space-time has thermodynamical behaviour. We found that the second law of thermodynamics is always satisfied for an expanding universe. It should be noticed that the presence of thermal equilibrium fluctuations and quantum fluctuations can contribute as the entropy correction, and consequently the number of degrees of freedom will be corrected, correspondingly. This work is under our current study and will be reported in the near future.