Surface gravity in spherically symmetric collapsing stars

Here we consider the generalized Oppenheimer-Snyder collapse of a star into a four-dimensional Einstein-Gauss-Bonnet black hole as well as a class of regular black holes labeled by the polytropic index of the stellar matter. We then analyze the nature of the horizon and the corresponding surface gravity outside and inside the star. The Hayward and Nielsen-Visser dynamical surface gravity are in agreement with the one resulting from the Killing vector of the outer static metric. However, these two definitions inside the star do not coincide with the Killing surface gravity outside the star when the star crosses the event horizon. This motivates us to study the surface gravity using Fodor's approach to have a unique surface gravity at the mentioned moment. Then the extremality condition and the first law of thermodynamics are discussed at the trapping horizon of the star.


Introduction
Many efforts have been made to analyze the theory of general relativity (GR) by studying its solutions and their physical interpretations.Black hole (BH) is one of the predictions of GR which is characterized by the existence of horizons, trapped surfaces, and singularities [1].
Although the naked singularities are covered by event horizons according to the cosmic censorship hypothesis [2], many efforts are made to introduce some BH geometries that are regular at all spacetime points.The first attempt was made by Einstein and Rosen by formulating the wormhole structure [3].In the next attempt, Sakharov proposed P = −ρ as the equation of state for the dense region of a star [4].One of the most important steps to remove the singularity was made by Bardeen by coupling the GR action to the nonlinear electrodynamics [5].Then, Hayward introduced a regular BH solution with a fundamental length scale that prevents the formation of a singularity [6].Both the Bardeen and Hayward metrics which are of interest here, reduce to the deSitter metric at small distances while they go to On the other hand, there are some extensions of GR which include functions of the curvature invariants in the action.One of the important generalizations of GR is the Gauss-Bonnet (GB) theory, generalized by Lanczos [9] and Lovelock [10].This theory prevents the Ostrogradsky instability [11] by giving second-order field equations.Although the consequences of GB gravity were supposed to be the same as those of GR in fourdimensional spacetime, recently it has been shown that a redefinition of the coupling constant of GB gravity can lead to new results even in four-dimensional spacetime [12].This theory is known as 4D Einstein-Gauss-Bonnet (4D-EGB) gravity.It should be noted that some authors have the opinion that the regularization method of [12] is questionable and propose new regularization methods [13], [14].Fortunately, even when other regularization methods are applied, the spherically symmetric BH solution of [12] is still valid.For this solution, the metric tensor is regular in the center.However, the scalar curvature diverges at this point.
To get a physical intuition of a regular BH, the Oppenheimer-Snyder-Datt (OSD) collapse scenario [15] for a general class of regular BHs is studied in [7] and for the 4D-EGB BH geometry in [16].In these works, the outer BH geometry is smoothly connected to the inner geometry of the star which is assumed to be spatially flat Friedmann-Robertson-Walker.It has been shown that the stellar matter could have uniform density and pressure to obtain a smooth transition at the star surface according to the Israel junction conditions [7].In the generalized OSD collapse, the stellar matter will violate the strong energy condition (SEC) in the final stages of the collapse.This prevents the star from reaching the center in the case of a regular BH while it reduces the collapse velocity for 4D-EGB collapse and the star reaches the center with zero velocity.As mentioned above, regular BHs are the solutions of non-vacuum Einstein equations but EGB BH is the so-lution of vacuum EGB gravity.Therefore, EGB-gravity boundary conditions should be used for the latter.In [17], the boundary conditions of charged anisotropic spherically symmetric stars in higher dimensional EGB gravity with a non-zero cosmological constant are obtained.Using the modified field equations, the authors of [17] have generated a complete stellar model in EGB gravity and show that the radial pressure at the stellar surface should be zero, similar to the case of GR.As expected, when the GB coupling constant tends to zero, the boundary conditions of N-dimensional GR are obtained.Here we have described 4D-EGB gravity as GR with an effective stress-energy tensor containing higherorder curvature terms.Thus, we have used the Israel junction conditions of GR for the collapse to 4D EGB BH [16].We will return to this point in Section 2.
In this paper, we want to investigate the nature of horizons in the generalized OSD collapse in terms of the ingoing and outgoing null geodesic congruences.The surface gravity (SG) is then calculated using the nonaffinely parameterized geodesic on the trapping horizon (TH) [18].
We also consider different approaches to SG, such as Hayward, Nielsen-Visser, and Fodor's approach.We see that it is possible to obtain a unique SG on the horizon according to Fodor's approach.To do this, we use null vectors with normalization coefficients that are arbitrary functions of spacetime, the affinely parameterized ingoing null vector, and the condition of the uniqueness of the SG at the crossing moment.
We study the extremality condition for the outer and inner geometry and see that outside the star, SG vanishing and the coordinate-invariant definition of extremality are the same.But this is not the case for the TH inside the star.We have also found the work term by writing the first law of thermodynamics using different definitions of SG.
The outline of this paper is as follows: In Section 2, we recall the OSD gravitational collapse for the 4D-EGB and the regular BHs.We consider different concepts of the TH and evaluate the SG in Section 3. In Sections 4 and 5, we study the expansion parameter and the Lie derivative of the expansion parameter and then we evaluate the properties of the horizon in the outer and inner geometry, respectively.We study the SG of the event horizon (outside the star) and TH (inside the star) using different approaches in Section 6.In Section 7, following Fodor's approach, we obtain a condition that gives a unique SG at the horizon crossing time by choosing appropriate null vectors.Then, in Section 8, we study the extremality condition for the external and internal geometries.Section 9 is devoted to obtaining the first law of thermodynamics for different approaches of Section 3. Section 10 provides a summary and conclusion.
Throughout this paper, the signature of the metric tensor is assumed to be (−, +, +, +).We use the natu-ral units in which G = c = 1.All variables are dimensionless with respect to the Schwarzschild radius.A dot denotes differentiation with respect to the proper time of a freely falling particle on the surface of the star and a prime means the derivative with respect to the radial coordinate.

OSD collapse into regular and 4D-EGB BHs
Consider a particular class of static, asymptotically flat spherically symmetric BH geometries as where the function ω(r) is defined by [7] for a regular BH with deSitter core and by [16] ω for a 4D-EGB BH.In the above expressions, β is a free parameter, n is the polytropic index of the stellar matter for the regular BH case, and α is the GB coupling constant.For n = 1 and n = 3/2, the metric (1) is reduced to the well-known Hayward [6] and Bardeen [5] metrics respectively while an arbitrary value of n gives a new family of static regular BHs with deSitter core.The location of the horizons of metric (1) can be obtained by setting rω(r) = 1.In the regular case, there are two horizons for β < 4/27, the BH becomes extremal for β = 4/27, and there is no horizon for β > 4/27.
In the case of the 4D-EGB, there are two horizons for α < 1/64π.The value of α = 1/64π corresponds to an extremal BH, while larger values of α lead to a horizonless geometry.
It is convenient to change the coordinates to the Painlevé-Gullstrand (PG) coordinates to construct a collapsing model for a star in the background geometry (1).In terms of the PG coordinates we can write the general form of the spacetime metric as follows where f (r, τ ) is given by and R is the radius of the star and the geometry inside the star is assumed to be described by the spatially flat FRW metric.Given (5), in order to have a smooth matching of the two geometries, we must have at the surface of the star.This means that the extrinsic curvature is the same on both sides of the star's surface and thus, there is no stress-energy layer on the surface of the star according to the Israel junction conditions.As mentioned before, here, we treat 4D EGB gravity as GR with an effective stress-energy tensor [16].In this way, the vacuum 4D EGB field equations can be written as where G µν = R µν − 1 2 g µν R is the Einstein tensor and includes the higher-order terms of curvature.This picture is usually used for any other modified gravity theory, and so one can express the field equations in the Einstein form G µν = 8πT eff µν where T eff µν = T EGB µν + T m µν is the effective stress-energy tensor.It includes the matter stress-energy tensor T m µν and the curvature terms that arise from the 4D EGB gravity, T EGB µν = −2αH µν .This interpretation has been used extensively in f(R) gravity [19], curvature-matter couplings [20], Weyl gravity [21], and in braneworlds [22].In this way, the effective stress energy tensor includes higher order curvature terms responsible for violating the energy conditions.Moreover, this picture allows us to use the junction conditions of GR.Otherwise, one should consider the EGB boundary conditions [17].Regarding the 4D EGB BH solution, if we insert (3) into ( 1) and then substitute it into Einstein's equations, we get ρ EGB and the components of anisotropic pressure p EGB .These must be added to the corresponding quantities of the star, ρ star and the components of p star in order to obtain the effective density and the different components of effective pressure of the star.This is done in [16] and it is shown that for a collapsing star, the radial component of the star's pressure would be zero at the surface of the star.This is similar to the result given by [17], but for a dynamically collapsing star.
According to (6), each free particle on the surface of the star begins its free fall from rest at infinity and moves along a time-like radial geodesic.Introducing the comoving radius as R(τ ) = a(τ )R com and substituting it into (6), this equation has the form of Friedman's equation, This gives the density of the star.It is given by for a regular BH and by for a 4D-EGB BH.Then, the pressure at the surface (and inside the star) can be read from the continuity equation, ρ + 3H(ρ + P ) = 0, which can be simplified to for a regular BH and for a 4D-EGB BH.In the interior region of the star, the event horizon and TH can be found using the outgoing radial null geodesics.The time evolution of the horizons and the surface of the star is studied in detail in [7] and [16].This is shown in the Penrose diagram of figure 1.
Assuming that the stellar matter is a perfect fluid, the equation of state can be derived from the geodesic equation of the stellar surface.In the first case, this leads to polytropic matter, while in the second, a different form is obtained.See equations ( 49) and ( 22) in [7] and [16] respectively.

Trapping horizon and the surface gravity
Consider a congruence of radial ingoing and outgoing null geodesics with tangent fields n a and ℓ a which are orthogonal to the 2-spheres of symmetry and satisfying the cross normalization n a ℓ a = −1.The corresponding expansion parameters are denoted by θ ℓ and θ n .The expansion parameter of ℓ a is defined as and n ↔ ℓ is substituted to obtain θ n , the expansion parameter of the ingoing null vector field.A TH [23] is a hypersurface that can be foliated by 2-surfaces such that θ n h < 0 , and past if θ n h > 0. A marginally trapped tube T (MTT) [24] is a hypersurface that is foliated by 2-surfaces S, called marginally trapped surfaces (MTS), such that θ ℓ T = 0 and θ n T < 0. If a MTT is spacelike (timelike) everywhere, it is called a dynamical horizon (time-like membrane) 1 .
For THs, the change in the area of the horizon can be related to the local value of T ab ℓ a ℓ b [23].From the Raychaudhuri equation with θ ℓ = 0 and assuming that l α is hypersurface orthogonal, one has where r a is a vector that is tangent to the TH and normal to the foliation by 2-surfaces, δA is an infinitesimal circle drawn on the spacelike 2-surface, σ ab is the shear tensor corresponding to the outgoing null vector and ζ is a scalar field at the TH.Assuming that ζ > 02 , since θ n < 0 and n b ∇ b θ ℓ < 0 for a future outer TH (FOTH) and n b ∇ b θ ℓ > 0 for a future inner TH (FITH), the sign of the area change is determined by the null energy condition 3 .The sign is negative for FITH and positive for FOTH if the null energy condition is satisfied.The existence of a decreasing area horizon is reported for a time-like membrane in some shell collapse models [26].
In stationary spacetime, the SG is a measure of the non-affinity of the Killing vector at the Killing horizon, where it becomes null.The SG of a non-Killing horizon in a spherically symmetric BH is usually defined as the inaffinity of the outgoing radial null geodesic ℓ a at the horizon.Several proposals have been made to fix the normalization of ℓ a .In this paper, we focus on three approaches 4 .
The first one is the proposal of Fodor et al [29] which is based on an affinely parametrized ingoing null geodesic n a whose asymptotic form is such that ξ a n a = −1 in an asymptotically flat geometry where ξ a is the asymptotic time-translational Killing vector.The crossnormalization condition n a ℓ a = −1 leads to Following Fodor's method, there will be no unknown degrees of freedom left in the ingoing null vector field and hence in the outgoing corresponding one, when the SG is read out in advanced Eddington-Finkelstein coordinates.So the SG is completely determined [29].However, this is not true in PG coordinates as we will see in Section 7.
The second is Hayward's proposal [30].It is based on the definition of the SG of spherically symmetric spacetime in terms of the Kodama vector K a [31].Here we want to explain the Kodama SG in detail.This is because there is a small error in this definition that has been repeated in the literature and needs to be clarified.Assume that the spacetime metric is given by where R is the areal radius and ǫ ab is the volume form of h ab .The Kodama vector is defined as To find the SG, it is necessary to calculate where in the last line the Hayward definition of SG and ( 19) are used.We see that ǫ bc appears on the right side of ( 20).The Kodama equation, which appears in [18] and [32] as well as in other literature, omits ǫ bc and uses K c on the right hand side.However, the presence of ǫ bc on the right side of ( 20) is expected since it is a consequence of the antisymmetrization of the left side of (20).It should be noted that the Kodama vector and the value of SG remain unchanged by this correction.Both Hayward and Fodor SGs are valid at all points in spacetime however we are usually interested in their values at a FOTH.The third definition is that of Nielsen and Visser [33] who use the PG coordinates.Consider the general metric form in PG coordinates The authors of [33] define the outgoing and ingoing radial null vectors as follows In the above, the normalization factors are chosen so that the null vectors satisfy two conditions.The first is the cross-normalization relation, n a ℓ a = −2 and the other is related to the SG of the TH.This function must be consistent with what is found in the first law of BH thermodynamics written by partial differentiation of the Misner-Sharp mass [34] m(τ, r) = rv 2 (τ, r)/2c 2 (τ, r) with respect to proper time.Doing this and then evaluating the result at the evolving horizon 2m(r, τ ) = r, we get This can be rewritten in terms of the area of the evolving horizon This is the first law of BH thermodynamics if The authors of [33] define the SG based on the inaffinity of the outgoing null vector field and as we have mentined, fix its normalization to obtain the SG (28).This gives ( 23) and ( 24).We will return to this point in Section 9.
The key concept to note is that the Hayward and Nielsen-Visser definitions of SG are not the same as Killing SG in static spacetime when g tt g rr = −1 in Schwarzschild-like coordinates [18], [35].Thus, these definitions are identical to the Killing SG for the static outer geometry of the star, given by (1).If the FRW geometry inside the star is reduced to the deSitter geometry where the Hubble parameter is independent of time, the above condition on the metric components holds.Otherwise, the Hayward and Nielsen-Visser SGs lead to different results, as we will see in Section 6.
The next two sections discuss the properties of the outer and inner geometry of the collapsing star.Then we will use the three definitions (17), (21), and (28) above to calculate the SG in Section 6.

External geometry of the star
Outside the star, the geometry is static.Therefore, the event horizon coincides with the OTH.The radius of the TH of metric (1) is given by r ± ω(r ± ) = 1 which gives Here we have considered the case with two horizons (i.e.β < 4/27 and α < 1/64π).Using metric (1) and crossnormalization n a ℓ a = −1, a specific set of radial null vectors is as follows By a simple calculation, the expansion parameter of n and ℓ using equation ( 15) can be obtained as where at the horizons which means that r ± are MTTs.Furthermore, using (1) it can be shown explicitly that The sign of the expression (36) depends on the sign of 1 + r 2 ± ω ′ (r ± ).For regular BH where . The expression (37) is positive in the outer horizon r + (so (36) is negative) and negative in the inner horizon r − for all allowed values of β and n.The same result is obtained for the 4D-EGB BH for all allowed values of α.This can be seen from and substituting r ±(EGB) from (30) in the above.Therefore, as expected, r + is the FOTH for both BHs.

Internal geometry of the star
Inside the star, spacetime is dynamical and so, the definition of the event horizon turns out to be substantially different from the quasi-local horizon definitions mentioned in Section 3. Using (4), and cross-normalization n a ℓ a = −1, a specific choice of the radial null vectors is given by It is then easy to show that in this region So, it is a TH.On this surface, the Lie derivative of the expansion parameter becomes Using Friedmann's equations, H 2 = 8πρ/3 and Ḣ = −4π(ρ + P ), we see that The pressure inside and on the surface of the star is negative, according to ( 13) and ( 14), for regular and 4D-EGB BHs.Therefore, expression (44) is positive and the TH is a FITH.We have already shown that the total energy-momentum tensor inside the star satisfies the NEC [7], [16].Therefore, the area of the horizon inside the star can be reduced even if the null energy condition is satisfied, according to (16).This is consistent with the result of the horizon evolution equation obtained in [7] and other collapse models [36].

Surface gravity
In this Section, we will calculate the SG for the inside and outside of the star.As mentioned in Section 3, the definition of Hayward's SG in terms of Kodama's vector and the definition of Nielsen-Visser's SG are not necessarily the same.To clarify this point, let us consider the following general metric in PG coordinate [35] where σ(t, r) is equal to one if g tt g rr = −1.In the static case, the Killing SG is: For the dynamical case, according to the definition of SG by Hayward and Nielsen-Visser, we obtain Thus, the two definitions give the Killing SG for static spacetime, such as the outer region of the star.However, they are not the same in the dynamical case due to the existence of ṁ, such as the inner region of the star, see (5).Below, we will see that these two definitions are identical for the TH inside the star if Ḣ = 0, the de Sitter geometry.
Here, both the inner and outer regions of the star are represented by the PG coordinates.Therefore, we use the general form of the spacetime metric as (4).First, we use Hayward's method to calculate the SG.Substituting the two-dimensional part of metric (4) into Hayward formula (21), we get Therefore, where it can be evaluated at each constant time slice of the TH1 .According to Fig 1, the TH is formed when the surface of the star crosses r + .At this moment, τ c , i.e. the crossing time, the radius of the star is R c where from ( 5) and a natural question is, are the SGs of (50) and (51) equal at this moment?To answer this question, we must express the inner SG (51) in terms of the density and pressure of the star by using Friedmann's equations and then use equations ( 9) and ( 12) to write (51) in terms of the ω function and its derivative at the star's surface.Performing this calculation, we find that at crossing time.This is different from (50) where its value at this moment is Now, it is necessary to compare the two metrics ( 4) and ( 22) in order to find the SG at crossing time using the Nielsen-Visser method.Remembering the definition of the Misner-Sharp mass (25), we then obtain from (28) on the horizon where f = 1, These can be simplified at the crossing time to The result of (57) shows that the Nielsen-Visser SG of the outer horizon is the same as (50), as expected.So it is obvious that at the crossing time, when the surface of the star reaches the outer event horizon, it is reduced to (54).This is because the geometry outside the star is static and has a time-like Killing vector.According to the standard definition of the Killing SG for Killing horizons, a simple calculation gives (50) as the Killing SG.
In general, however, the SG evaluated on the r + and TH of the star will not necessarily coincide at the crossing time.This is because the metric inside the star is inherently dynamic for which there is no unique definition of SG.This point can be easily seen from ( 53) and (58).

A unique surface gravity at crossing time
Now we want to study SG using Fodor's approach.In this approach, the ingoing radial null rays satisfy the affinely parametrized geodesic equation2 .To obtain the same inner and outer surface gravities at the crossing time, we use some arbitrary normalizing functions in the null vectors and then specify these functions.The outgoing and ingoing radial null vectors for metric ( 4) are Assume that the geodesic equation for n a is affinely parameterized.Then we obtain the following equation for D(r, τ ): where (4).Substituting these symbols into (61) yields with solution where E, E ′ , γ and γ ′ are integration constants.As mentioned in Section 3 for an asymptotically flat geometry with the time-translational Killing vector ξ, we require that ξ a n a = −1 at spatial infinity, therefore This means that (66) From the zeroth component of the outgoing radial null geodesic, ℓ a ∇ a ℓ b = κ F ℓ b , we have Now, if we use the cross-normalization ℓ a n a = −1, we have 2CD = 1 and equation (68) becomes This gives the SG at r + and TH using ( 4), ( 66) and ( 67) where in the second line of (70), the definition of the Hubble parameter H = Ṙ(τ )/R(τ ) is used.Also, we choose E = 1 to match (71) with the SG given by ( 50) and (57) for the outer metric.Evaluating (70) and (71) at the crossing time gives: Now, we find the coefficient E ′ with the condition that (72) and ( 73) are equal at the crossing time.At this moment, using ( 9) and ( 12), the right-hand side of (73) can be rewritten in terms of the density and pressure of the star and thus from Friedmann's equations in terms of the Hubble parameter and its derivative, i.e.
Substituting this expression into (73) and equating it with (72), yields where in the second equality we have used ( 9), ( 12), (52) and Friedman's equations.Now, the Fodor's SG (70) is determined by the value of (75) as showing that Fodor's SG is proportional to the stellar contraction velocity.Our results are summarized in Table 1.For two typical examples of regular BHs, Hayward and Bardeen, as well as for 4D-EGB, the external SG is shown in Table 2. Substituting (75) into (67), the null vectors (60) for r < R become and for r > R , 0, 0) (79) From (6), it is clear that the SG (76) is, in general, a function of proper time It is therefore of interest to study the evolution of the SG as the star contracts.The time dependence of the SG: evaluation at the event horizon outside the star SG: evaluation at the TH inside the star Rc τc<τ <τe -SG: inside the star at crossing time Bardeen BH (n = 3/2 in (2)) radius of a star collapsing into the regular BH, (2), and 4D-EGB BH, (3), has been calculated in [7] and [16] as in which where 2 F 1 is the hypergeometric function and the integration constant R 0 is chosen such that τ = 0 at R = R 0 .
Combining the above expressions with (2), ( 3), ( 29) and ( 30), we thus find the time dependence of the SG.It is worth noting that according to (82)-( 84), in the ex-tremal case (i.e.β = 4/27 and α = 1/64π) the radius of the star reaches the extremal horizon in a finite time, and the outer SG (73) becomes zero in this case.The evolution of the SG is shown in figure 2 in the time interval [τ c , τ e ].The TH appears at τ c and then vanishes at τ e .The maximum of the SG occurs when the surface acceleration of the star is zero.It should be remembered that in static spacetime, SG is actually the force that an observer at infinity would have to exert to keep an object at the horizon.This description does not hold in the dynamical case.Here we see that, according to (76), the SG of a spherically symmetric collapsing star is maximized when its radial acceleration is zero.This happens on a surface between r + and the TH of the star.In other words, inside the star, according to Friedman's equation the SG is a function of the density and pressure of the star.A comparison of ( 76) and (85) shows that the surface gravity is at a maximum when the SEC is marginally valid.This is somewhere between r + and TH of the star.Moreover, in the 4D EGB and regular BHs, the SG of the TH increases when the SEC is satisfied, and vice versa.This is in contrast to the case of gravitational collapse into the Schwarzschild BH, where the SG of the TH always increases.This is because the SEC is not violated at any radius of the star.For a star collapsing to a near extremal BH, the time interval in figure 2 becomes very small.In the limit of an extremal BH, the SG exists only at a certain value of the proper time and its value is zero.This can be well deduced from (76) where the numerator is zero for an extremal BH.

Extremality condition
The extremality condition for stationary BHs is that the Killing SG must be zero.This condition is equivalent to the fact that in (4) f ′ (r) = f (r) − 1 = 0 has a double root r + = r − = r ex and there is no trapped surface.A coordinate-invariant definition of extremality for dynamical horizon, given in [35] as By comparing (4) and ( 22), the outgoing radial null vector field of Nielsen-Visser ( 23) can be written as Then performing a simple calculation leads to where f = 1 at the horizon.From (5), the extremality condition (86) then reduces to for the outer BH metric.This is exactly the same as (57).Recall that from (4), f ′ (r) = 0 for an extremal static BH, so the extremality condition is trivially satisfied and this is equivalent to setting the SG (57) of the outer horizon to zero.For the dynamical metric inside the star Therefore the vanishing of this value is not equivalent to the vanishing of the SG of the TH (56).
Performing a similar calculations in Fodor's approach using (60), we find that Thus, for the r ex and TH of the star, we have where the last equality of ( 94) is a result of (68) and gives κ (out) F = 0.However, a comparison of (95) with (76) shows that the inner SG is not zero if the extremality condition is satisfied inside the star.

The first law of thermodynamics
In this section, we will consider the first law of thermodynamics for the evolving TH.Using the time derivative of the Misner-Sharp mass, the first law of thermodynamics takes the form where −w 0 dV TH dτ is the work term and the function w 0 must be determined.For each slice of constant time, the Misner-Sharp mass within the TH is with derivative Moreover As mentioned before, in the Nielsen-Visser [33] approach, the SG is obtained by assuming that the work term in the first law of thermodynamics does not exist for the TH.But we should be careful that ( 27) is the first law of thermodynamics with a partial time derivative of mass and surface area which gives a different result compared to the situation where we use the total derivative.In the latter case, substituting (98) and (99) into the first law (96) gives which is the SG of the dynamical horizon, as defined in [37] and differs from the Nielsen-Visser SG (56), as expected.
A similar calculation for Hayward's SG (51), leads to and for the SG (70) derived from Fodor's approach where in the second equality we have used the Friedmann equation and the third equality is written in terms of the star's radius.From our discussion above, it is easy to see that different definitions of SG lead to different versions of the first law of thermodynamics, as expected.In [30], Hayward writes the laws of thermodynamics for THs and uses the Kodama vector to define SG, resulting in a non-zero work term in the first law of thermodynamics.Ignoring the work term, Nielsen and Visser in [33] found a relationship similar to the first law by partial differentiation of the Misner-Sharp mass and defining the surface change coefficient as SG.On the other hand, the first law of thermodynamics is expressed in terms of the total derivative of the thermodynamic quantities.Therefore, the result of [33] may not be applicable without considering the work term in the first law of thermodynamics with total derivatives.For example in [29] using Einstein's equations, it is shown that Fodor's definition of SG in the first law of thermodynamics is not applicable without considering the work term.Therefore, in the above, we have considered the work term in the first law of thermodynamics for Fodor's definition of SG.As mentioned above, if we take the total time derivative of the horizon equation or the Misner-Sharp mass, regardless of the work term, the coefficient of variation of the area of the horizon will be equal to 1/2r T H , which is the SG of the dynamical horizon and defined in [24].

Conclusion
This paper considers the collapse of a star into a special class of regular and 4D-EGB BHs.The interior of the star is described by the spatially flat FRW geometry.To find the SG for evolving THs, we have considered several approaches.These include those of Hayward, Nielsen-Visser and Fodor.We have studied the SG for the outer and inner geometries using the above approaches.Since the outer geometry is static, the SG at the outer event horizon is the same in all approaches.However, due to the dynamical nature of the inner geometry, different approaches lead to different SGs.Tables 1 and 2 summarize the results for regular and 4D-EGB BHs.
We have also seen that different SGs are not necessarily equal to the outer SG at the moment when the surface of the star crosses the outer horizon.Following Fodor's approach, we can use null vectors with arbitrary normalization coefficients which are then determined by cross-normalization and SG uniqueness conditions at the crossing time.The crossing time is defined as the moment at which the star's surface, the event horizon and the TH coincide with each other.In this way, we obtain the time evolution of the inner SG for regular and 4D-EGB BHs.It has a maximum when the collapsing acceleration of the star becomes zero.This is completely different from the behavior of the SG in the OSD collapse to a Schwarzschild BH where the SG evolves uniformly and then diverges asymptotically as the surface of the star approaches the singularity.
Applying the extremality condition, we find that it does not necessarily lead to zero SG for the interior geometry.Here, we saw that using the PG coordinate time, a finite time is required for a star to collapse into an extremal BH.We then obtained the first law of BH thermodynamics for the evolving TH for each of the above mentioned approaches.This gives a different first law from that of Nielsen and Visser.

Figure 1 :
Figure 1: Penrose diagram of a star collapsing into a BH with two horizons.The horizons are denoted by r+ and r−.The blue line is the time-like TH.It appears when the surface reaches r+ and disappears when it reaches r−.A straight line and a wavy line indicate the final state of the collapse for the singularity-free case in the regular metric and with a singularity in the 4D-EGB case.

Figure 2 :
Figure2: Evolution of the SG as a function of proper time for Schwarzchild (ω = 1), Hayward (n = 1), Bardeen (n = 3/2) and 4D-EGB BHs.We set β = 0.09 for the regular BH and α = 0.0005 for the 4D-EGB BH.For the Schwarzschild BH, the SG grows uniformly and diverges at the singularity, where the TH disappears.For the regular and 4D-EGB BHs, however, the SG is finite.Its values at r− and r+ are the same.In between, it has a maximum value.The variables are dimensionless with respect to the Schwarzschild radius.

Table 1 :
SG of horizons of a collapsing star using different approaches.Outside the star, all of these approaches lead to the same result at the horizon.The TH appears at τc and then vanishes at τe.The variables are dimensionless with respect to the Schwarzschild radius.

Table 2 :
Different SGs for regular and 4D-EGB BHs.The variables are dimensionless with respect to the Schwarzschild radius.