Universal Thermodynamics in the Context of Dynamical Black Hole

The present work is a brief review of the development of dynamical black holes from the geometric point view. Furthermore, in this context, universal thermodynamics in the FLRW model has been analyzed using the notion of the Kodama vector. Finally, some general conclusions have been drawn.


Introduction
A black hole is a region of space-time from which all future directed null geodesics fail to reach the future null infinity I + .More specifically, the black hole region B of the space-time manifold M is the set of all events P that do not belong to the causal past of future null infinity, i.e., B = M − J − (I + ). ( Here, J − (I + ) denotes the causal past of I + , i.e., it is the set of all points that causally precede I + .The boundary of the black hole region is termed as the event horizon (H), H = ∂B = ∂(J − (I + )). ( A cross-section of the horizon is a 2D surface H(Σ) obtained by intersecting the event horizon with a space-like hypersurface Σ.As event the horizon is a causal boundary, it must be a null hypersurface generated by null geodesics that have no future end points.
In the black hole region, there are trapped surfaces that are closed 2-surfaces (S) such that both ingoing and outgoing congruences of null geodesics are orthogonal to S, and the expansion scalar is negative everywhere on S. Therefore, the 3D boundary of the space-time region that contains trapped surfaces is the trapping horizon.Further, a closed two-surface on which one of the congruences has zero expansion is termed as a marginally-trapped surface or apparent horizon.For a stationary black hole, the apparent horizon coincides with the event horizon, while in the dynamical situation, the apparent horizon always lies inside the black hole region unless the null energy condition (NEC) is violated.
According to Hawking and Penrose, the existence of trapped surfaces inside a black hole implies the formation of a singularity provided the strong energy condition (SEC) is violated [1].Geometrically, the existence of singularity implies geodesic (time-like or null) incompleteness.
A stationary black hole that is static must be spherically symmetric [2][3][4][5][6].On the other hand, the stationary space-time of a rotating black hole is necessarily axially symmetric, and it admits two Killing vectors t α and φ α , i.e., the Lie derivatives of the metric tensor along the vectors vanish, more specifically L t g αβ = 0 and L φ g αβ = 0, respectively.Then, the linear combination ξ α = t α + Ω H φ α is null on the event horizon (Ω H is the black hole's angular velocity).In fact, ξ α is tangent to the horizon's null generators.The surface gravity κ of the black hole is defined by the relation [7]: on the horizon.Note that κ is constant along the horizon's null generators, i.e., uniform over the entire horizon.Now, we shall consider a congruence of time-like geodesics with four velocity u α .The divergence of the four velocity is called the expansion scalar, i.e., θ = u α ;α .The rotation tensor is given by: where The shear tensor is given by . Now, using the facts (i) due to hypersurface orthogonal nature of congruence of null generators, the rotation tensor ω αβ ≡ 0, and (ii) due to the stationarity of the event horizon, both θ and dθ dλ are both (λ is the proper time), the Raychaudhuri equation for null geodesics [8,9]: implies (using Einstein field equations): (a) T αβ ξ α ξ β = 0 on the event horizon.Therefore, matter cannot flow across the event horizon, i.e., the generators would get focused and the black hole would not be stationary.(b) σ αβ = 0, i.e., the null generators of the event horizon have a vanishing shear tensor.Usually, a black hole is characterized by its event horizon.However, the very definition of event horizon demands that it could never be observed by any means in finite time; its existence can only be located in infinite future.According to Hayward, the event horizon does not have any physical effect [10], and it may pass through an observer at any finite time without any imprint.Consequently, the notion of trapped surfaces may be interesting in the context of a black hole, as they may hint that a black hole has been formed.
In 4D space-time, let us consider 2D surfaces that are space-like and closed (i.e., topologically, they are two-spheres).Now, at each point on this 2D surface, one may consider two null directions normal to the surface, i.e., there are two distinct future directed families of null geodesics emerging from the surface; along the outer and inner directions.If it is so happens that both the families of null geodesics orthogonal to the surface converge, then the two-surface is said to be trapped.Geometrically, it implies that both the null geodesics (orthogonal to the surface) emitting from every point of the surface will decrease in area.More specifically, these two-surfaces are termed as future trapped surfaces.The convergence of outgoing null geodesics implies that there is a strong gravitational field, tending to focus on the geodesics, and something is unusual due to the presence of a trapped surface.
Let ζ be a closed space-like two-surface (topologically two-sphere) and contained in a 3D space-like hypersurface Σ. Suppose s is a unit space-like vector normal to ζ, but lying on Σ. Similarly n, is a future-directed unit time-like normal field to Σ.Then, K ± = n ± s are the outgoing and ingoing future-directed null normal fields with normalization: K + K − = −2.Then, the expansion scalars θ ± of the congruence of null geodesics having tangent vector fields K ± are defined as: where γ AB is the first fundamental form, i.e., the induced metric on the two-surface ζ and K AB is the second fundamental form, i.e., extrinsic curvature of the two-surface ζ.The above scalars measure to what extent the congruence of null geodesics expand, namely θ > 0(< 0) implies the divergence (convergence) of the congruence of null geodesics.Thus, the trapping nature of ζ is characterized by the signs of θ ± .It should be noted that the values of θ ± depend on the choice of the null normals, which cannot be uniquely defined; rather, they depend on the choice of the hypersurface Σ.However, the signs of θ ± do not depend on the choice of the null normals K ± [11].Thus, the compact spatial two-surface ζ with future pointing null directions K ± (normal to ζ) is said to be a future (past) trapped surface if θ ± | ζ < 0 (θ ± | ζ > 0).However, future trapped surfaces are not quite characterized as a black hole.Therefore, one has to consider marginal surfaces that are just trapped to characterize the boundary of a black hole.
A marginal surface is a spatial two-surface S on which one of the null expansions vanishes (say θ + | S = 0).A trapping horizon is the closure H of a three-surface H foliated by marginal surfaces on which θ − | H = 0 and L − θ + | H = 0.The trapping horizon, as well as marginal surfaces are said to be: Therefore, for a black hole, the most relevant horizon is the future outer trapping horizon having properties:
Therefore, the future outer trapping horizon characterizes a quasi-local black hole horizon.Hawking's area theorem, a fundamental feature of a classical black hole bounded by an event horizon, can be extended to a quasi-local black hole.In fact, the future condition θ − | H < 0 guarantees the non-decrease of the horizon area [11].
As the geometric structure of null and spatial hypersurfaces is distinct, it is interesting to study the stationary and dynamical regimes of a quasi-local black hole, characterized by a future outer trapping horizon.In stationary regimes, the geometric structures constructed on a null hypersurface that is foliated by closed outer marginally-trapped surfaces is provided by isolated horizons.In fact, different levels of stationarity for a black hole horizon (in a dynamical environment) are characterized by an isolated horizon.It should be noted that an isolated horizon does not interact with its surroundings, i.e., there is no exchange of information with the surroundings.The event horizons of static Schwarzschild and Reissner-Nordström black holes are isolated horizons.
The dynamical horizon characterizes the non-stationary regime of quasi-local black hole.A space-like hypersurface foliated by closed future marginally-trapped surfaces is termed a dynamical horizon.
The Killing horizon, which is purely local in nature, is defined by the vanishing of the norm of a Killing vector field.The null hypersurface, which is the connected component of the set of points where a Killing field is null, is the Killing horizon.Note that a Killing horizon is not necessarily a black hole horizon, and also, it is not an isolated horizon.Subsequently, Hayward [10][11][12][13] initiated the study of black hole thermodynamics for dynamical black holes.He introduced a one-form equation known as the unified first law (UFL), and he had shown its equivalence with the Einstein field equations.It was also shown that by projecting the UFL along the tangent to the trapping horizon, the first law of black hole physics can be obtained for the dynamical black hole.
The paper is organized as follows: In Section 2, we discuss some basic ideas about dynamical black holes.Section 3 deals with the thermodynamics in the homogeneous and isotropic FLRW model.The derivation of the Clausius relation from UFL is presented in Section 4. A new definition of the surface gravity is discussed in Section 5. Section 6 is devoted to concluding remarks.

Preliminaries on Dynamical Black Holes
The discovery of Hawking that the surface gravity determines the effective temperature of a quantum field on a stationary black hole background created a nice analogy between the laws of thermodynamics and certain properties of black holes.As a result, one has to take into account black hole dynamics (i.e., dynamically-evolving black holes) instead of black hole statics.Usually, a black hole is characterized by the event horizon for which the global assumption of asymptotic flatness is essential, and this may not be the case for the real universe.
Hayward, in a series of papers [10][11][12][13], introduced the notion of a dynamical, quasi-local definition of black holes in terms of trapping horizons, which are hypersurfaces foliated by marginal surfaces.He then derived a second law of black hole dynamics by expressing the increase of the area element along the trapping horizon, and consequently, black hole dynamics was formulated as the dynamics of trapping horizons.Hayward introduced a new definition of surface gravity for non-stationary black holes to formulate a corresponding first law of black hole dynamics.Instead of using the Kodama vector, he introduced the surface gravity similar to the stationary definition with the Killing vector on the Killing horizon.In this work, we have considered the 3 + 1 space-time in the framework of Einstein's gravity.
There is another debated issue in general relativity, namely the idea of energy, as there is no well-accepted definition of gravitational energy.However, for spherical symmetry, Misner-Sharp energy has all the relevant physical properties, and it represents active gravitational energy.This energy is defined as [14]: where # stands for the contravariant dual (index raising) with respect to the space-time metric g and the dot denotes contraction.Here, the Newtonian gravitational constant is chosen to be unity, and spatial metrics are taken to be positive definite by sign convention.One can construct two invariant quantities from the energy tensor T as: and the vector where the trace is taken over the 2D space normal to the spheres of symmetry.Here, w is termed as the work function, and ψ is called the energy supply vector.Then, the gradient of the Misner-Sharp energy is related to the Einstein equations as: where denotes the covariant dual.The 2D space normal to the spherical symmetry is usually characterized by one spatial and one temporal direction, which cannot be chosen uniquely.However, using double-null coordinates (ξ + , ξ − ), each symmetric sphere has two preferred normal directions (the null directions) ∂ ∂ξ ± .The line element for 4D space-time in the double-null form can be written as: where dΩ 2 2 is the metric on unit two-sphere and r ≥ 0 and f are a function of the null coordinates.In this coordinate system, the different components of the UFL are given by: where the following Einstein field equations are used to derive the last equation.

•
First law of black hole statics: In stationary space-time ∃, a Killing vector ξ that generates asymptotic time translations and ξ satisfies the Killing equation: with ⊗ denoting the symmetric tensor product.A Killing horizon is a null hypersurface for which the Killing vector is tangent to the null generators, and it may be considered as the boundary of a stationary black hole (or a white hole/cosmological regions).Now, by definition, the vector ξ • (∇ ∧ ξ ) is parallel to ξ on a Killing horizon, and the proportionality constant is termed as surface gravity, i.e., on a Killing horizon with ∧ denoting the antisymmetric tensor product.
As an example, one may choose the Reissner-Nordström solution of the Einstein-Maxwell equations.The two parameters e and m characterize the charge and asymptotic energy m.For m ≥ |e|, it describes a stationary black hole with surface gravity κ = √ m 2 −e 2 r 2 on the Killing horizon: Therefore, on the horizon, one gets the first law of black hole statics: with 'd' the differential in the state space (m, σ).Due to the 'm' term, the above first law cannot be generalized to space-times that are not asymptotically flat.By introducing the quasi-local energy E for the present Reissner-Nordström solution as: the above first law can be rewritten as: This is a physically realistic form of the first law of black hole statics with electromagnetism.In analogy to the internal energy of a black hole, E is the energy on the horizon (rather than at infinity) or, equivalently, E is the energy of the black hole, not the energy of the space-time and is defined as the active gravitational energy.
In black hole dynamics, there is no time-like Killing vector to define surface gravity.Kodama [15] introduced a vector, which for spherically symmetry is defined as: with curl defined on the 2D space normal to the spheres of symmetry.For the double null co-ordinate system, the Kodama vector takes the form: The Kodama vector satisfies (i) k • ∇r = 0 and (ii) k • k = 2E r − 1.The Kodama vector is spatial, null or temporal for trapped, marginal or untrapped spheres, respectively.Therefore, similar to a stationary black hole where the boundary is the hypersurface with the null Killing vector, the horizon of a dynamic black hole is the hypersurface where the Kodama vector is null, i.e., the trapping horizon.Thus, it is possible to have a preferred flow of time for a general spherically-symmetric space-time by the Kodama vector (which coincides with the Killing vector for stationary space-time), and active gravitational energy is the conserved charge associated with the Kodama vector; consequently, it is speculated that dynamic gravitational entropy might be defined in analogy to the method of Wald [16].Now, it is natural to define dynamic surface gravity by replacing the Killing horizon with the trapping horizon and the Killing vector by the Kodama vector.Therefore, the surface gravity is defined as: on a trapping horizon ∂ ± r = 0 or equivalently a purely geometrical definition: where the divergence and gradient are taken over the normal 2D plane of the spheres of symmetry.Then, one may speculate that a dynamic spherically-symmetric quantum field has a local Hawking temperature hκ 2π , but it remains an open question.Subsequently, in the context of black hole dynamics with the above generalized definition of surface gravity, the first law of black hole dynamics can be written as: where prime stands for differentiation along the trapping horizon: f = z • ∇ f with z a tangent vector to the trapping horizon.Here, the first term on the r.h.s. is the usual term for the black hole dynamics, while the second term denotes the work done along the tangent to the trapping horizon.
In black hole statics, the zeroth law for a stationary black hole states that the surface gravity κ is constant.In black hole dynamics, the thermal equilibrium is characterized by the constancy of: where S is the entropy of the system.Therefore, gravithermal equilibrium is defined as the thermal equilibrium in a stationary space-time, and a weak type unified zeroth law is defined as both κ and β are constants.

Universal Thermodynamics in the FLRW Model
In universal thermodynamics, it is interesting to study the FLRW space-time from the point of view of a dynamical black hole.It is found that for this homogeneous and isotropic space-time, there is no outer trapping horizon; rather, there is an inner trapping horizon, which coincides with the apparent horizon [10][11][12][13][17][18][19].It is found [17] that the Clausius relation can be obtained by projecting the UFL along the tangent to the trapping horizon.Due to absence of the time-like Killing vector field for non-stationary models, one can introduce the divergence-free Kodama vector [15] field for any time-dependent spherically-symmetric space-time.This Kodama vector acts as an asymptotically time-like Killing vector field.Furthermore, there is a conserved current associated with this vector field [15,20].In the following, the universal thermodynamics for the FLRW model will be discussed in the context of the Kodama vector [21].
We start with the homogeneous and isotropic four-dimensional non-stationary FLRW space-time: where a(t) is the scale factor of the expanding universe, k is the curvature and R = a(t)r is the area radius.This metric can also be written as: where We can write this metric in terms of null co-ordinates as: where R is a function of l and m, respectively.Suppose the space-time is time orientable, and we consider the future pointing null vectors ∂ ∂l , ∂ ∂m , which are given by: The expansion of the geodesic of congruences is given by: Depending on θ l and θm, we can have three kinds of surfaces: i.e., for a trapped surface, both null expansions converge.In other words, in a trapped surface, both the incoming and outgoing light rays converge to the surface of the sphere [22].For an untrapped surface, one of the expansions of the null congruence converges, and one of expansions of null congruence diverges.Explicitly, the incoming light rays converge, and the outgoing light rays diverge.In the case of a marginal surface, at least oneof the expansions of null congruences vanishes.Now, for the trapped surface, two situations may occur: both θ l , θ m > 0 and both θ l , θ m < 0. If both θ l , θ m > 0, i.e., both 1 then the surface is called the past trapped surface, and if both θ l , θ m < 0, i.e., both 1 then the surface is called the future trapped surface.For a marginal surface with θ l = 0, we may have the following cases: • bifurcating: θ m = 0, • outer: ∂θ l ∂m < 0, • inner: ∂θ l ∂m > 0, • degenerate: ∂θ l ∂m = 0.A hypersurface foliated by a non-degenerate marginal hypersurface is called a trapping horizon.For the marginal surface of the future outer trapping horizon of black holes, we have: while for the marginal surface of the past outer trapping horizon of white holes, we have: For the FLRW model, the horizon is taken in analogy to the trapping horizon of a dynamical black hole.Since gravity is negative, we have to deal with the inner horizon.
The trapping horizon is defined as (∂ l R) R=R T = 0, i.e., where R A is the geometric radius of the apparent horizon and H is the Hubble parameter.We can write the internal energy inside a sphere of area radius R totally in terms of the area radius, i.e., the internal energy is totally a geometric quantity.Now, the internal energy can be written as [14]: For the FLRW model, the internal energy can be written as: i.e., for the FLRW model, the internal energy can be written in terms of the radius of the apparent horizon and the area radius.According to Hayward [10][11][12][13], the UFL can be written as: i.e., the change in internal energy can be written as the sum of the work function and the work done for dV amount of volume change.Here, A = 4πR 2 is the surface area of a sphere of area radius R and V = 4 3 πR 3 is the volume of that sphere; and the work density, is regarded as the work done by a change of the horizon, while the energy-supply term is given by [23,24]: ψ = Ψ a dx i denotes the flow of energy through the horizon of area radius R. The Kodama vector for the FLRW model is given by [15,25]: where ab is the two-dimensional Levi-Civita tensor in the radial temporal plane.For the FLRW model, the Levi-Civita tensor is given by: and the Kodama vector for the FLRW model is explicitly given by: The Kodama vector serves as the time-like Killing vector for dynamical black holes and FLRW space-time (as there is no time-like Killing vector in non-stationary FLRW space-time).In de Sitter space, it is very similar to the Killing vector ∂ ∂t .For the present FLRW model, the work density is given by: and the energy supply term is given by: The work done for dV amount of volume change is given by: and the energy flow across the surface area A of a sphere of area radius R is given by: Therefore, the change in internal energy E is given by: Now, we shall obtain the Friedmann equation directly from the UFL by projecting it along the Kodama vector.
Projecting the internal energy change along the Kodama vector, we obtain: Projecting the heat flow across the area A along the Kodama vector, we obtain: WdV, K i = 0.
Therefore, the projection of UFL along the Kodama vector gives: i.e., we obtain the second Friedmann equation.
We shall now show the first Friedmann equation by projecting it along a vector U µ orthogonal to the Kodama vector in the two-dimensional radial temporal plane.The above-mentioned vector is given by: The vector has the following properties: (i) Depending on the components, the vector may be space-like, null or time-like.(ii) The divergence of this vector vanishes (i.e., ∇ µ U µ = 0).Therefore, we can find a current associated with that vector given by ξ µ = G µν U ν .From the Bianchi identity, we can say the vector ξ µ is conserved, i.e., the divergence of the vector vanishes, i.e., in mathematical form: ∇ µ ξ µ = 0.
Projecting the change in internal energy along the vector U µ gives: Projecting the heat flow along the vector U µ gives: and the projection of the work done by dV amount of volume along the vector U µ gives: Therefore, projecting the UFL along the vector orthogonal to the Kodama vector gives: • projecting along the vector 2 HR 3 , 0 and 0, − 2 3R 2 a Projecting the UFL along the vector 2 HR 3 , 0 , we obtain: Projecting the UFL along the vector 0, − 2 3R 2 a , we obtain: • Projecting along the vector U µ and Ūµ Let us consider two linearly independent vectors U µ = (V 0 , V 1 ) and Ūµ = (V 0 , −V 1 ) in FLRW space-time where V 0 ,V 1 are arbitrary functions of r and t with V 1 = 0. We now project the UFL along these two vectors: respectively.Adding Equations ( 57) and (58), we obtain (55), and subtracting (58) from (57), we obtain (56).These two vectors will be divergence free simultaneously if we take where φ 2 (r) and φ 2 (t) are arbitrary functions of r and t alone, respectively.As these vectors are divergence free, then we can find two conserved currents C µ = G µν U ν and Cµ = G µν Ūν associated with these vectors U ν and Ūν , respectively, i.e., ∇ µ C µ = 0 and ∇ µ Cµ = 0. Depending on the components, both of these vectors may be time-like, space-like or null.

Clausius Relation from the UFL
It has been shown that first law of black hole thermodynamics can be obtained by projecting the UFL along the trapping horizon [10][11][12][13]17], i.e., where z is a vector in the direction tangent to the trapping horizon.We shall now show the validity of the Clausius relation for the event horizon.
The normal vector to the hypersurface R = R E is given by n a = (−1, a, 0, 0).The vector n a is a null vector.This vector is also tangential to the event horizon.Then, we can show that the Clausius relation is not satisfied for the event horizon.
The flow of energy across the surface of the sphere of area radius R E during infinitesimal time dt is given by: Using the second Friedmann equation, the above relation can be written as: where we have used . Now, to derive the Clausius relation, we shall use the generalized Hawking temperature [26]: . Therefore, we have the right-hand side of the Clausius relation as: From ( 62) and (64), we can easily obtain the Clausius relation δQ = T (G) E dS E on the event horizon.The present approach to obtain the Clausius relation is a more generalized way, and this can be applied for any horizon.

A Redefinition of Surface Gravity
The surface gravity at the trapping horizon for the FLRW model is defined as: Here, we assume the above definition of surface gravity, and we shall examine the interrelationship between the UFL and Friedmann equations.For any horizon, the surface gravity is defined as: or explicitly, we can write it as: We know that the Hawking temperature on the horizon of area radius R h is given by T = 1 2πR h ; the amount of heat flow along this surface of area radius R h is given by dQ = 4πR 3 h (ρ + p)dt; and the horizon entropy is S h = πR 2 h G ; then, we can see that dQ = TdS.Therefore, we need to modify the Hawking temperature.In this section, we redefine surface gravity for a local Rindler observer.
The Unruh temperature T = h f 2πck B is proportional to the acceleration [27], and the form is very similar to the Hawking temperature of a black hole.If A is the acceleration of the FLRW model, so in analogy, we take the temperature to be proportional to the acceleration of the model (h = c = k B = 1).
Here, A is defined as: where we have introduced a 0 as a dimensionless constant of proportionality and R is introduced to match the dimension on both sides.Using the Friedmann equation, the temperature on the local Rindler causal horizon becomes: Here, we introduce a new vector V µ similar to the Kodama vector: Projecting the change in internal energy along the above-mentioned vector, we obtain: The right-hand side in the above expression is nothing but T RH dS, V µ , i.e., we obtain the Clausius relation: i.e., δQ = T RH dS.
where we have taken a 0 = 3 2 and S is the Bekenstein entropy.It is very interesting to note that the above-mentioned vector has properties similar to the Kodama vector.
(i) V µ lies in the 2D r − t plane.The vector is neither parallel, nor orthogonal to the Kodama vector, as well as U µ .(ii) Depending on the components, the norm of the vector may be positive, negative or zero, i.e., the vector may be space-like, time-like or null.(iii) The divergence of the vector vanishes.Therefore, we can find a conserved current C µ = G µν V ν associated with the vector, i.e., ∇ µ C µ = 0.

Conclusions
In this paper, we have described different aspects of the UFL for homogeneous and isotropic FLRW space-time.We have derived the Friedmann equation from the UFL in three different ways.Firstly, we have projected the UFL along the Kodama vector, and we have obtained the second Friedmann equation, then projected it along a divergenceless vector perpendicular to Kodama vector in the two-dimensional r − t plane and obtained the first Friedmann equation.In the next method, projecting the UFL along the vectors 2 HR 3 , 0 and 0, − 2 3R 2 a , we have obtained the Friedmann equation for pressure and density components separately.These two vectors are linearly independent, and they form a basis of the r − t plane.In the final method, we have projected the UFL along the vectors U µ = (V 0 , V 1 ) and Ūµ = (V 0 , −V 1 ), which may be divergenceless for certain choices and obtained from the Friedmann equations in a combined form.It is also interesting to note that these two vectors form a basis of the two-dimensional r − t plane for V 1 = 0.
We have also shown the validity of the Hayward-Kodama definition of surface gravity (65) for dynamical FLRW space-time.Assuming the generalized Hawking temperature, we can also obtain the Clausius relation from the UFL.Finally, we have redefined the surface gravity as a quantity proportional to the acceleration of FLRW space-time, and we can obtain the Clausius relation by projecting the UFL along the divergence-free vector V µ .In the future, we shall generalize the work to more generalized space-time models.